Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 497,862 C→T intergenic (+309/‑201481) unknown → / → unknown tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1497,8620CT92.5% 184.5 / 3.0 80intergenic (+309/‑201481)unknown/unknowntRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (5/1);  new base T (39/35);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

GCCGACTATCAGCTTCGGTTGCAGCGTTGCCACCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCC  >  1/497618‑498054
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                 
gCCGACTATCAGCTTCGGTTGCAGCGTTGTCACCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCaa                                                                                                                                                                                             >  1:698424/1‑251 (MQ=255)
                 gTTGCAGCGTTGTCACCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAAc                                                                                                                                                                                           <  1:72046/236‑1 (MQ=255)
                         gTTGTCACCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGt                                                                                                                                                                    >  1:442439/1‑251 (MQ=255)
                               aCCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGt                                                                                                                                                              >  1:490786/1‑251 (MQ=255)
                                cccGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTg                                                                                                                                                             >  1:689712/1‑251 (MQ=255)
                                cccGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTg                                                                                                                                                             >  1:656487/1‑251 (MQ=255)
                                cccGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTg                                                                                                                                                             >  1:845976/1‑251 (MQ=255)
                                  cGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGa                                                                                                                                                                                 <  1:639863/229‑1 (MQ=255)
                                  cGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAg                                                                                                                                                           >  1:199425/1‑251 (MQ=255)
                                   gATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAgg                                                                                                                                                          >  1:783682/1‑251 (MQ=255)
                                    aTAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGa                                                                                                                                                                                 <  1:304500/227‑1 (MQ=255)
                                    aTAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATc                                                                                                                                                                 <  1:711753/243‑1 (MQ=255)
                                                gCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAg                                                                                                                                                                                         <  1:333025/207‑1 (MQ=255)
                                                   tGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAgg                                                                                                                                                                                        >  1:41959/1‑205 (MQ=255)
                                                      ttACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGa                                                                                                                                        <  1:649568/250‑1 (MQ=255)
                                                         cATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTa                                                                                                                                                                   >  1:682684/1‑220 (MQ=255)
                                                                 aTTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTtct                                                                                                                             >  1:79872/1‑250 (MQ=255)
                                                                           aTTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTa                                                                                                                  <  1:253078/251‑1 (MQ=255)
                                                                                   aTTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTgg                                                                                                                                                                     >  1:49601/1‑192 (MQ=255)
                                                                                   aTTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAgg                                                                                                                                                          <  1:723290/203‑1 (MQ=255)
                                                                                        gttgATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAAc                                                                                                                                                                                           <  1:457569/165‑1 (MQ=255)
                                                                                                cAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGaa                                                                                              <  1:619853/250‑1 (MQ=255)
                                                                                                 aaCGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAag                                                                                            <  1:54247/251‑1 (MQ=255)
                                                                                                   cGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAagag                                                                                          <  1:342170/251‑1 (MQ=255)
                                                                                                     ggCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAggg                                                                                        >  1:244119/1‑251 (MQ=255)
                                                                                                          cAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGa                                                                                                                                                                                 <  1:499887/157‑1 (MQ=255)
                                                                                                            ggTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCg                                                                                                                                 <  1:598584/203‑1 (MQ=255)
                                                                                                              tAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg                                                                                <  1:359912/250‑1 (MQ=255)
                                                                                                                ggCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAg                                                                                                                         >  1:177415/1‑207 (MQ=255)
                                                                                                                ggCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGt                                                                             >  1:349809/1‑251 (MQ=255)
                                                                                                                ggCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGt                                                                             >  1:40265/1‑251 (MQ=255)
                                                                                                                ggCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGt                                                                             >  1:34491/1‑251 (MQ=255)
                                                                                                                    aTTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑aaaa                                                                        <  1:381233/251‑1 (MQ=255)
                                                                                                                            aCTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGaa                                                                                                                                                   <  1:675338/169‑1 (MQ=255)
                                                                                                                                  aaTTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATaa                                                           <  1:369192/250‑1 (MQ=255)
                                                                                                                                    ttAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAcc                                                        <  1:529715/251‑1 (MQ=255)
                                                                                                                                     tAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAccc                                                       <  1:8185/251‑1 (MQ=255)
                                                                                                                                                    aTTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt                                        <  1:330172/251‑1 (MQ=255)
                                                                                                                                                    aTTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt                                        <  1:764886/251‑1 (MQ=255)
                                                                                                                                                      tttCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg                                      <  1:367703/251‑1 (MQ=255)
                                                                                                                                                        tCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg                                    <  1:173318/251‑1 (MQ=255)
                                                                                                                                                        tCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg                                    <  1:126220/251‑1 (MQ=255)
                                                                                                                                                         ccTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAAt                                                                                  <  1:809433/205‑1 (MQ=255)
                                                                                                                                                         ccTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAcc                                                        >  1:65219/1‑230 (MQ=255)
                                                                                                                                                         ccTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAc                                                         >  1:700199/1‑229 (MQ=255)
                                                                                                                                                           ttATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAg                                 <  1:735380/251‑1 (MQ=255)
                                                                                                                                                            tataAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAccc                                                       <  1:148805/228‑1 (MQ=255)
                                                                                                                                                            tataAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAcc                                                        >  1:723118/1‑227 (MQ=255)
                                                                                                                                                            tataAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAg                                 >  1:563030/1‑250 (MQ=255)
                                                                                                                                                             ataAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTgg                                                                                                                                                                     <  1:824448/118‑1 (MQ=255)
                                                                                                                                                                  cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAgg                                                                                                                                                          >  1:755469/1‑124 (MQ=255)
                                                                                                                                                                  cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGACAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAccc                                                       >  1:453298/1‑222 (MQ=255)
                                                                                                                                                                  cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTtatat                                                             >  1:243785/1‑216 (MQ=39)
                                                                                                                                                                  cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTtatat                                                             >  1:169288/1‑216 (MQ=39)
                                                                                                                                                                  cacCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTACGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGa                          >  1:89271/1‑251 (MQ=37)
                                                                                                                                                                   gatctgcaccccaaaacctggacacGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGAt                         >  1:342208/25‑251 (MQ=11)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCAtt                                                                                                                    >  1:210003/1‑160 (MQ=255)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGcc                                                                                                       <  1:506908/173‑1 (MQ=255)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCt                                                                                                      >  1:528298/1‑174 (MQ=255)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAag                                                                                            >  1:90337/1‑184 (MQ=255)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg                                                                                >  1:329445/1‑196 (MQ=38)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg                                                                                <  1:768899/196‑1 (MQ=38)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAAt                                                                       >  1:635466/1‑204 (MQ=38)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTtatat                                                             >  1:413223/1‑214 (MQ=38)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAc                                                         <  1:586638/218‑1 (MQ=38)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg                        >  1:249194/1‑251 (MQ=25)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg                        >  1:292542/1‑251 (MQ=25)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg                        >  1:551074/1‑251 (MQ=25)
                                                                                                                                                                    ccTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATg                        >  1:645188/1‑251 (MQ=25)
                                                                                                                                                                        gcaccccaaaacctggacacGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAcc                                                        <  1:411520/196‑1 (MQ=11)
                                                                                                                                                                           ccccaaaacctggacacGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAAt                                                                                  >  1:783098/17‑187 (MQ=11)
                                                                                                                                                                           ccccaaaacctggacacGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaag                  >  1:744744/17‑251 (MQ=11)
                                                                                                                                                                           ccccaaaacctggacacGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg                                      >  1:508810/17‑230 (MQ=2)
                                                                                                                                                                              gCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaagaaag              >  1:883966/1‑251 (MQ=17)
                                                                                                                                                                                 cctcAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGCGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCACCTTATATAAAccc                                                       >  1:612535/4‑207 (MQ=9)
                                                                                                                                                                                   ccAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAag                                                                                            <  1:187019/169‑1 (MQ=17)
                                                                                                                                                                                      aaaCCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGaa      <  1:539103/251‑1 (MQ=14)
                                                                                                                                                                                       gacacGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGa                          >  1:629087/5‑230 (MQ=11)
                                                                                                                                                                                         cccGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt                                        >  1:444670/1‑214 (MQ=14)
                                                                                                                                                                                          ccGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATcc  <  1:540026/251‑1 (MQ=11)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                 
GCCGACTATCAGCTTCGGTTGCAGCGTTGCCACCCGATAGAAGCCATCGCATGGTTACATTTGGCATTTGAATCCATTCACCCATTTTGTTGATGGCAACGGGCGACAGGTAGGCTATTGCTAAACTTTGAATTAATGAAAAACGGTTATTTTCCTTATAAACACCTATGATCTGCCCCCCAAAACCCGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCC  >  1/497618‑498054

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑