Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,308,9950TA25.8% 126.1 / 12.1 62intergenic (‑71/+38)nlpI/pnplipoprotein NlpI/polynucleotide phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (31/15);  new base A (0/16);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 2.20e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCATACGTGCCAGAATCACTTCCTGCTGTAAAGTCGGTTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCAGCGATTTGAGAGATGTGGACCAGACCTTCTTT  >  NC_000913/3308750‑3309236
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
cATCGGTTAAAGCCCGACTGGCAAGGAGTTGTTCCATACGGGCCAGAATCACGTCCTGCTGTAAAGTCGGTTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGGCGCAACGAAACGCCAGCGCAAAAAAGGCTTCATGTCCCACTCCCGAAGACCACGGTTGAAAGAACTGCCTGTTCCCGGTTGCGAACAAGGCGTCCTGCCCGCTTaaaa                                                                                                                                                                                                                                                <  2:227519/250‑1 (MQ=255)
 aTCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCATACGTGCCAGAATCACTTCCTGCTGTAAAGTCGGTTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGc                                                                                                                                                                                                                                             <  1:127198/251‑1 (MQ=255)
      ttAAAGCCCGACTAGCCAGGGTTTGTTCCATACGTGCCAGAATCACTTCCTGCTGTAAAGTCGGTTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGAATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCCACGAAACACCAGCGCACAAAAGGCTTCATTTCCCACTCCCGCAGAACACGGTTCAATGAACGTCCTGTTCCCGCTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGcccccc                                                                                                                                                                                                                                        >  2:46314/1‑251 (MQ=255)
        aaaGCCCGACTGGCAAGGATTTGTTCCATACGTGCCAGAATCACTTCCTGCTGTAAAGTCGGTTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCccgc                                                                                                                                                                                                                                      >  1:561797/1‑251 (MQ=255)
         aaGCCCGACTGGCAAGGATTTGTTCCATACGTGCCAGAATCACTTCCTGCTGTAAAGTCGGTTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAAGGTCCTGTTCCCGGTTGCTAACAAGGCGGCCTGCCCGGTTAAAATCCCgccgcc                                                                                                                                                                                                                                     <  1:409542/251‑1 (MQ=255)
                  tGGCAAGGATTTGTTCCATACGTGCCAGAATCACTTCCTGCTGTAAAGTCGGTTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCgg                                                                                                                                                                                                                              >  1:487800/1‑249 (MQ=255)
                                      cgtgccaaaagctctgcctgtggTAAAGGCGGTTGCAATGGTCCCGGGAGGATTTCACTTTTCCGCCAGGAAGGATTACGGAATCCGGCAAGGGTAAGTGTGGTCACAACGAACCCCCGGCGAAAAAAAGGTTTCTTTCCCCACCCCCGAAGACCACGGTTAAATGAACGTCGTGGCCCCGGTGGCAAACAAGGCGGCCTGCACGTTAAAAACCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTc                                                                                                                                                                                                        <  2:494788/229‑1 (MQ=255)
                                       gTGCCAGAATCACTTCCTGCTGTAAAGTCGGCTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTCCATCCTGCCAGCGTAAGTGCTGACGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTACCCAGGCGTCCTGCCCGGTTAAAAGCACCCCGCCGCAGCGGCGGGCAAAGGGCAACCTTACTct                                                                                                                                                                                                       >  1:34859/1‑250 (MQ=255)
                                       gTGCCAGAATCAATTCCTGCTGTAAAGTCGGTTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCaa                                                                                                                                                                                                                       >  1:152101/1‑235 (MQ=255)
                                                            gTAAAGTCGGTTGAAATGGTACCGCGAGTACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGAAAGGTTAAGTGCTGTCCCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACCCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGTTTGCTAACAAGGCGTCCTGCCCGTTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAg                                                                                                                                                                                    <  1:159787/249‑1 (MQ=255)
                                                                     gTTGCAATGGTACCGCGAGGACTTCGCTTTGACGCCAGGAAGTATTACTGCACCCGGCAAGCGTAAGTGCTGTCGCAACAAAACACCAGCGCAAAAAAGGCTTAATTTCCCCCCCCCGAAAACCACGGTTAAATGACGGTCCTGTTCCGGGTTGCAAACAAGGCGTCCTGCCCGGTAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTc                                                                                                                                                                                               <  1:487935/229‑1 (MQ=255)
                                                                     gTTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACCCCAGCGCAAAACAGGCTTCATTTCCCACTCCCGAAGACAACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCACCCTTACTCGCCCTGTTc                                                                                                                                                                                               >  2:487935/1‑229 (MQ=255)
                                                                       tGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtgctgc                                                                                                                                                                       >  1:493520/1‑251 (MQ=255)
                                                                        gCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGCTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtgctgc                                                                                                                                                                       >  1:474222/1‑250 (MQ=255)
                                                                             ggTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCAACCTGCAAGCGTAAGTGCTGTCGCAACGAAACCCCAGCGAAAAAAAGGCTTCATTTCCCCCCCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCAAACAAGGCGTCCTGCCCGTTAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAgg                                                                                                                                                                 <  1:322340/251‑1 (MQ=255)
                                                                                      aGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTg                                                                                                                                                        >  1:439924/1‑251 (MQ=255)
                                                                                        gACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCCAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGGGCTTCCGGTGCTGCAGCAGGTTGAGACTGCt                                                                                                                                                      >  2:310834/1‑251 (MQ=255)
                                                                                             cacgtttACGCAAGGAAGTATTCCTGCACCCTCCAAGCGTAAGTGCTGTCCCAACGAAACACCAGCACAAAAAAGCCTCCATTTCCCCCTCCCGAAGACCACGTTTGAATGACCGTCCTTTCCCCCGTTCCTAACAGGGCCTCCTGCCCGGTTAAAGGCCCCCCGCCGCAGCGGAGGGCAAATGGCATCCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTc                                                                                                                                                 <  1:288015/247‑1 (MQ=255)
                                                                                              aCTTGTACGCCAGGAAGTAGTACTGCATCCTGCAAACGTAAGTTCTGTCGCAACGAAACACCAGCGCAAAAAAGTCGTGGTTTCCCACGCCCGAAGTCCACGGTTGAATGACGGGCGTGTTCCCGGTTGCTAGCAAGGCATCCTGCCCGGTTAAAAGccccc                                                                                                                                                                                                                                         <  1:513382/162‑1 (MQ=255)
                                                                                                  ttACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAAAACCGTCGACACGTAACACCAGCGCAAAACAAGCTTACTTTCCCACTCCCGAAGAACACGGTTGAATGAACGACCTGATCCCGGTTGCTCACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGAAGCGGAAGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTcc                                                                                                                                                                               >  1:505346/1‑216 (MQ=255)
                                                                                                           gAAGTATTACTGCATCCTGCACGCGAACGCGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCt                                                                                                                                   >  1:454038/1‑251 (MQ=255)
                                                                                                           aaaGTATTACTGCATCCTGCAAGCGTAATTGCTGACGCAACGAAACCCCAGCGCAAAAAAGGTTTCATTTCCCCCGCCCGAAGACCAGTGTGAAAGGACCGCCCTGTTCTGGGGTGCTAACAGGGCGTCCTGCCCGGTAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCGTACTCGCCCTGTTCAGCAGCCGGAGCGTGCGGTGCTGCAGCAGGTTTAGACGGCTCAGTCGCTTCTTTAATGCt                                                                                                                                   <  2:394088/250‑1 (MQ=255)
                                                                                                             aGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCACAAAAGGCTTCATTTCCCCCTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCCGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTcc                                                                                                                                 >  2:476956/1‑250 (MQ=255)
                                                                                                                         cccctgaacggaatgTGCGGTCGCAACGAAACACCGGCCAAAAAAGGGCTCCATTCCCCACCCCCGAGAACCACGGTGGAATGAACTCCCTGTCCCCGGTTGCTAACAAGGCGTCCTGCCCGGTAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCATGTTGAGACTGCTCAGTCGCTGCTTTAATGCTCAGACGGATACGGc                                                                                                                     <  2:291179/237‑1 (MQ=255)
                                                                                                                               aagcgtagttccggccGCAACGAAACACCAGCGCAAAAAAGGCTTCATGTCCCACCCCCGAAGCCCCGGTTGTAATGAACGTCCTTTTCCCGGTTGCTAACAAGGCGTCCTCCCCGGTTAAAACCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGc                                                                                                               <  2:493520/236‑1 (MQ=255)
                                                                                                                               aaGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACTGCCCTGGc                                                                                                               >  1:321563/1‑251 (MQ=255)
                                                                                                                                aGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCg                                                                                                              >  1:156875/1‑251 (MQ=255)
                                                                                                                                  ggTAAGTGCTGTCGCACCGAACCCCCAGCGAAAAAAAGGCTACATTCCCCACCCCCAAAGACCACGGTGGAATGAACGCCCTGTTCCCGGTTGCAAACAGGGCGTCCTGCCCGGTTAAAACCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGAt                                                                                                            <  1:9969/250‑1 (MQ=255)
                                                                                                                                     aaGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGATTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCa                                                                                                          >  1:345856/1‑250 (MQ=255)
                                                                                                                                            gTCGCAACGAAACACCAGCGCCAAAAAGGCTTCATTTCCCACTCCCGCAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTCACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCACATGGCAACCTTACTCGCCCTGTTCAGCCGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCACCTTcc                                                                                                    >  2:398009/1‑249 (MQ=255)
                                                                                                                                              cGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGa                                                                                                 >  2:266628/1‑250 (MQ=255)
                                                                                                                                                caacaaacccccAGCACAAAAAAGGCTTAATTCCCCCCCCCCAAATCCCTCGTTTGAAGAACCGTCCTGTCCCCGGTTGCTAAAAAGGCGCCCTCCCCGGTAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTGGAGACTGCTCAGGCGCTTCTTTAATGCTCAGACGGATACGTCCCTGGCGATCAACTTCCAGAACt                                                                                              <  2:345856/241‑1 (MQ=255)
                                                                                                                                                        cacCAGCGCAAAAAAGGCTTCATTTCCAACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCgg                                                                                      >  1:489100/1‑251 (MQ=255)
                                                                                                                                                           cAGCGCAAAAAAGGCTTCATTCCCCACCCCCGAAGACCACGGTTGAATGAACGTCCTTTTCCCGGTTGCAAACAAGGCGTCCGGCCCGGTAAAAACCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTAc                                                                                   <  1:121724/251‑1 (MQ=255)
                                                                                                                                                               caaaaaaaGGCTCCATTTCCCCCTCCCTAACCCCACGTTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGCCTTCCTCCCTGTTTAAAACCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTcc                                                                               <  1:25426/249‑1 (MQ=255)
                                                                                                                                                               aaaaaaaaGGCTTCTTTCCCCACTCCCAAAACCCAGGGTTGAGTGAACTCCCTGTCCCCGGTTCTTAACAAGGCGTCCTCCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCACCCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATGCGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTcc                                                                               <  1:398009/249‑1 (MQ=255)
                                                                                                                                                                    aaaGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGAcc                                                                          >  1:494204/1‑251 (MQ=255)
                                                                                                                                                                                ttCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAggc                                                               >  2:280242/1‑250 (MQ=255)
                                                                                                                                                                                 ccccACCCCAGAAGACCAGGGTTGAAGGACCCTCCTGTCCCCGGTCGCTACAAAGACGTCCGGCCCGGTAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACGGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTGCACCGGTACTTCCTGACCCATCTGCAggga                                                              <  2:451265/250‑3 (MQ=255)
                                                                                                                                                                                  cccccccccaaaGACCAGGGTTGAATAACCGTCTTGGTCCCGGTTGATAACAGGGCGTCCTGCCCGGTAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACGTCCGGACCCATCTGCAGGTa                                                              <  1:476956/240‑1 (MQ=255)
                                                                                                                                                                                     ccccccgaaacccgcGGTTGAATGAACTTCCTTTTCCCGTTTGTTAACAAGCCCCCCTCCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCg                                                          <  1:306957/237‑1 (MQ=255)
                                                                                                                                                                                     ccccccGAAGCCCCCGGTGGAATGAACGCCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTAAAACGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCg                                                          <  2:439924/248‑1 (MQ=255)
                                                                                                                                                                                          cGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCAc                                                     >  1:511375/1‑251 (MQ=255)
                                                                                                                                                                                           aaagcccACGGTTGAATGACCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTAAAAACCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGGACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACt                                                    <  1:385331/246‑1 (MQ=255)
                                                                                                                                                                                            aaGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCAAACAAGGCGTGCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTAATCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGcc                                                                                                                    >  1:56393/1‑186 (MQ=255)
                                                                                                                                                                                             aGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACttt                                                  >  1:445291/1‑251 (MQ=255)
                                                                                                                                                                                                             gaaCTTCCTGTTCCCGTTTGCAACCAAGGCGTCCGGCCCGGTAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCTGAGCTTTCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCGTTAATGCTCGGACGGATACGGGCTTTTCGATCAACGTCCAGTACTTTCACGGGTACTTTTTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGt                                    <  1:285704/250‑1 (MQ=255)
                                                                                                                                                                                                                   tCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCAGCg                               >  1:79209/1‑250 (MQ=255)
                                                                                                                                                                                                                     cggTTCCCGGTTACTAACAAGGCGTCCTCCCCGGAAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAAGGCGCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCAGCGAtt                            <  2:106212/249‑1 (MQ=255)
                                                                                                                                                                                                                       tgcccctggtGCTAACAGGGCGTCCCCCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGGCAGCGAtt                            <  2:487800/240‑1 (MQ=255)
                                                                                                                                                                                                                       gTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCAGCGATTTg                          >  1:15849/1‑251 (MQ=255)
                                                                                                                                                                                                                       gTTCCCGGTTGATAACAAGGCGTCCTCCCCGGTAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCAGCGATTTg                          <  1:50277/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGAAATCGGTCACTTTCTCAACGCGTTTGTCAGAGATTTgag                        >  1:27181/1‑248 (MQ=255)
                                                                                                                                                                                                                             gtttGCTAACAAGGCGCCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCGCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTGCACCGGTACTTCCTGACGCAGCTGCAGGTAAGCGGTCACTTTCTCAACGCGTTTGTCAGCGATTTGAGAGAt                    <  1:67382/249‑1 (MQ=255)
                                                                                                                                                                                                                             ggTTGATAACAAGGCGTCCTGCCCGGTTAACCGCCCCACGCCGCAGCGGCGGGCACATGGCCACCTTACTCGCACTGTTCAGCAGCCGGAGCTTCAGGTGCTGCAGCAGGTTGAGCCTGCTCAGTCGCTTCTTTAATGCTCAGCCGGATACGGCCCTGGCGATCAACTTCCCGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGACACTTTCTCACAGCGTTTGTCAGCGATTTgtg                        >  2:230636/1‑245 (MQ=255)
                                                                                                                                                                                                                                    aacagggcccccggCCCGGTAAAAACCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTGCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCAGCGATTTGAGAGATGTGGAcc             <  2:564958/238‑1 (MQ=255)
                                                                                                                                                                                                                                       aaGGCGTCCTGCCCGGTTAAACGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCAGCGATTTGAGAGATGTGGAcc             >  2:501034/1‑248 (MQ=255)
                                                                                                                                                                                                                                       aaGGCCTCCTGCCCGGTTAAAAGCCCCCCGCCTCAGCGGAGGGCAAATGGCAACCTTACTCGTCCTGTTCAGCAGCGGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGAGCATCTTCGAGATCTTGCACCGGTGCGTCCGGACGCATCTGCAGGTAAGCTGTCTCTTTCTCAACGCGTTTGTCAGCGATGTGAGAGATGTGGACCAGa          <  1:367726/251‑1 (MQ=255)
                                                                                                                                                                                                                                        agggcacccgCCCGGTTAAAAAGCCCCCCGCCGCATCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGTATACTGCCCTGTCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAAGCGGTCACTTTCGCAACGCGGTTGTCATCTATTTGAGAGGTGTTCACCAGa          <  2:247206/242‑1 (MQ=255)
                                                                                                                                                                                                                                          gCGTCCCGCCGGGTAAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGTAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTGTCACCGTTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCATCGATTTGAGATATGTGGACCAGACct       <  2:213499/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCCAATGGCACCCTTACTCGCCCTGTTCAGAAGCCGGAGATTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCAGCGATTTGAGAGAtg                   >  1:348437/1‑237 (MQ=255)
                                                                                                                                                                                                                                                tGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCAGCGATTTGAGAGATGTGGACCAGACCTTCttt  >  1:142814/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
CATCGGTTAAAGCCCGACTGGCAAGGATTTGTTCCATACGTGCCAGAATCACTTCCTGCTGTAAAGTCGGTTGCAATGGTACCGCGAGGACTTCACTTTTACGCCAGGAAGTATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTTCCTGACCCATCTGCAGGTAATCGGTCACTTTCTCAACGCGTTTGTCAGCGATTTGAGAGATGTGGACCAGACCTTCTTT  >  NC_000913/3308750‑3309236

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: