Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,331,5570CT79.1% 67.0 / 19.0 45G352G (GGC→GGT) proPosmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (19/15);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-02
Rejected as consensus: Frequency below/above cutoff threshold.

CACCAACGTGACTTACTACATGTTGCTGACCTATATGCCGAGTTATTTGTCGCATAACCTGCATTACTCCGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATC  >  NC_000913/4331311‑4331800
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
caccaATGTGACTTACTACATGTTGCTGACCTATATGCCGAGTTATTTGTCGCATAACCTGCATTACTCCGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTg                                                                                                                                                                                                                                                   >  1:110526/1‑250 (MQ=255)
           cTTACTACATGTTGCTGACCTATATGCCGAGTTATTTGTCGCATAACCTGCATTACTCCGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCggg                                                                                                                                                                                                                                       >  1:351127/1‑251 (MQ=255)
                 acATGTTGCTGACCTATATGCCGAGTTATTTGTCGCATAACCTGCATTACTCCGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGGTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTg                                                                                                                                                                                                                                 <  2:30186/251‑1 (MQ=255)
                            aCCTATATGTCGAGTTATGTGTTGCAGAACCTGCATTACTCCTAAGACCACGGGGTGCTGATTATTAGCGCCATTATGATTGGTATGCGGTTTGTCCAGCCGGGG‑ATGTGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCGACTTGGTAGTGTTGGCCTGGTTTTGGTTGAGATCGCGGCTTTTATTCTGATTATGAGTGTCGTCATCGGCGTGATGGTTGCAGGGTGACTGATGCTGGCGGt                                                                                                                                                                                                                      <  2:365513/251‑1 (MQ=255)
                            aCCTATATGCGGAGTTATTTGTCGCATAACCTGTATTGCTCCGACGACCTCGGGCTGCTGATTATTCTCGCCATTTTGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTTCTGAGTGACCGATTTGGCCTTCGTCCGTTTGTGCTACTTGTTAGTGTTGCCCTGCTTGGGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTAGCCGGGTTACTGATGCTGGCggg                                                                                                                                                                                                                      >  2:211782/1‑250 (MQ=255)
                                 tttGCCGAGTTATTTGTCGCATAACCTGCATTATTCCGAAGCCCCCTGGTTGTTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTT‑ATGTGCTTGCTGATTTACCTGTTTGGCCGTCGTCCGTTTGTGCTACTTGGTATTGTTGCCTAGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCATATGTTTGCCGGTTTGCTGAttcttgttgtgatc                                                                                                                                                                                                                  <  2:342869/248‑14 (MQ=255)
                                      cGAGTTATTTGGCGCATAACCTGCATTACTCCGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACGGATGCTGGCGGTGATCCTTAAc                                                                                                                                                                                                            <  1:154219/251‑1 (MQ=255)
                                           tATTTGTCGCATAACCTGCATTACTCCGAAGACCACGGGGTGCTGCTTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGCTTAACAGTAACGTCATCGGTCTGATTTTTGCCAGGTTACTGATGCTGGCGGTGATCCTTAACTGCtt                                                                                                                                                                                                       >  2:24921/1‑251 (MQ=255)
                                            aTTTGACGCATTACCTGCATTACTCCGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGCAATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGTGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGCTACTGATGCTGGCGGttgtatttacctgtttt                                                                                                                                                                                                      >  2:320514/1‑235 (MQ=255)
                                             tttGTCGCATAACCTGCATTACTCCGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGACTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTTACTGCtttt                                                                                                                                                                                                     >  1:75585/1‑250 (MQ=255)
                                                        aCCTGCATTACTCCGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATg                                                                                                                                                                                          >  1:166761/1‑251 (MQ=255)
                                                                tACTCCGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGGTCCTTAACTGCTTTACGGGCGTTATGGCTTCTAc                                                                                                                                                                                  >  1:6509/1‑251 (MQ=255)
                                                                     cGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGc                                                                                                                                                                             <  1:213575/251‑1 (MQ=255)
                                                                       aaGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGccc                                                                                                                                                                           >  1:272492/1‑250 (MQ=255)
                                                                                  gTGCTTATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTTCGGGCGTTATGGCTTCTACCTTGCCAGCGATGTGccc                                                                                                                                                                <  2:110596/251‑1 (MQ=255)
                                                                                     cTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGGTCCCGAc                                                                                                                                                             >  1:318973/1‑251 (MQ=255)
                                                                                          tattatCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCat                                                                                                                                                         >  1:200344/1‑250 (MQ=255)
                                                                                            ttatCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTTCCGACGCATATc                                                                                                                                                      >  1:248793/1‑251 (MQ=255)
                                                                                                   cATTATTATCGGTATTCTGTTTGTCCAGCCGGTG‑ATGGGCTGGCTGAGTGACCGTTTGGGCCTTCGTCCGTTTGTGCTACTTGTTAGTGTTGCCCTGTTTGTGTTTGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTATCTGCTTTACGGGCGGTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCAGATCCGATACa                                                                                                                                               <  2:379164/251‑1 (MQ=255)
                                                                                                      ttggATCGGGATGCTTTTTGTCCAGCCT‑TGTATGGGCTTGCTTAGTGTCCGTTTTGGCCGTCGGCCGTTTGTGCTACTTGGTAGTGTTGCCCTTGTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATGTTTGCCGGGTTACTGATGCTGGCGGTGATCATTTACTGCTTTACGGGCGTTATTGCTTCTACCTTGCCAGCGATGTTCCCGACGCatagctgatacagc                                                                                                                                             <  2:422242/247‑12 (MQ=255)
                                                                                                      tATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTTTCCGTATTGGCCGTCGTCCGTTTTTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCGTGCCAGCGATGTTCCCGACGCATATCCGATACAgc                                                                                                                                             <  1:383517/250‑1 (MQ=255)
                                                                                                            cGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGGGACCGTTTTGGCCGTCGGCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGGTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTGGCCGGGTTACTGATGCTGGCGGTGATCCTTAACGGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCTATGTTCCCGACGCATATCCGTTACAGCGCGCTg                                                                                                                                       <  2:398251/250‑1 (MQ=255)
                                                                                                                aTGCTGTTTGTCCAGCCGGTT‑ATTGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGGCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTAGGGCGTCTACCTTGCCAGGGATGTTCCCGACGCATATCCGATACAGCGCGCTggcggc                                                                                                                                  <  2:200929/251‑1 (MQ=255)
                                                                                                                aTGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTACCAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGGTCCTTAAACGCTTTACGGGCGTTGTGGCTTCTAACTTGCCAGCCATGTGTTCGACGCAAATCCgataacacgcgcgggcggg                                                                                                                                  >  2:309031/1‑233 (MQ=255)
                                                                                                                       ttggccAGCCGTTG‑ATGGGCTTGCTGAGGGGCCGTTTTGTCCGTCGGCCGTTTGTGCTACTTGGTGGTGTGGCTCTGTTTGTGTGGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGAAGCTGGCGGTTCTCCTTATCTGCTTTACGGGCGTTATGGCTTCTCCCTTGCCAGCGCTGTGCCCGACGCATATCCGATTCAGCGCGCTGGCGGCGGCAtt                                                                                                                            <  2:200652/246‑1 (MQ=255)
                                                                                                                       ttGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCAttt                                                                                                                           >  1:319473/1‑251 (MQ=255)
                                                                                                                         gTCCAGCCGGGG‑AGGGGCTTGCTGAGTGACCGTTTTGGCCTTCGTCCGTGTGTGCTGCGTGGTAGTGTTGCTCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGAAGCATATCCGATACAGCGCGTTGGCGGCGGCATTTa                                                                                                                          <  2:303893/250‑1 (MQ=255)
                                                                                                                               ccGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTTCCCGCGATGTTCCCGACGCATATCCGATACAGCGCGCCTGCGGCGGCCTTTAATAtttt                                                                                                                   >  2:312142/1‑250 (MQ=255)
                                                                                                                               ccGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTc                                                                                                                   <  2:6509/251‑1 (MQ=255)
                                                                                                                                ggggg‑ATGGGCTGGCTGGGTGACCGTTTTGGCCGGCGTCCGTTTGTGCTCCTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCGTCTACCTTGGCAGCGATGTTCCCGACGCATAGCCGATACAGCGCGTTGGCGGCGGCATGTAATATTTCg                                                                                                                  <  2:351127/247‑1 (MQ=255)
                                                                                                                                  gTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTGGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTGGCCGGGTTGCTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCGTGCCATCGATGTTCCCGACTCAGATCCGAGACAGCGCGCGGGCGGCGGCATTTAATAGTTCGGt                                                                                                                <  2:183357/251‑1 (MQ=255)
                                                                                                                                      aTGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCt                                                                                                             >  1:348280/1‑251 (MQ=255)
                                                                                                                                           cTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATGAACAGTAACGTCATCGGGCTGATTTTTGCCGGGGTACGAATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCtt                                                                                                                                                                               <  1:17551/180‑1 (MQ=255)
                                                                                                                                            ttGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTTCCGACGCATATCCGATACAGCGCGCTGGCGGCGGGATTTAATATTTCGGGGCTGGTTGc                                                                                                       >  1:290812/1‑251 (MQ=255)
                                                                                                                                             tGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCCACGCATATCCGATCCAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGcc                                                                                                      >  1:57249/1‑251 (MQ=255)
                                                                                                                                                                 cgtcggcagttggggttACTTGGTAGTGGTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGTCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTggc                                                                                 <  2:360875/235‑1 (MQ=255)
                                                                                                                                                                           ttGTGCTCCTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCctggct                                                                        <  2:319473/250‑1 (MQ=255)
                                                                                                                                                                                    ttGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATTGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCCTTTAATATTTCGGTGCTGGGTGCCGGGTTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGc                                                              >  1:392847/1‑251 (MQ=255)
                                                                                                                                                                                        tAGTGTTACCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTTTGGCTTCTACATTTCCCGGGATGTTCCCGACGCATACTCGATAAAACGCGCTGGCGGCGGGCTTTTATATTTCGGTGCTGGGTGCCGGTCTGCCGCCCAAC‑CTggcgggcagggcgggcggaagcccgc                                                          >  2:350320/1‑230 (MQ=255)
                                                                                                                                                                                                                        cGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCAGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATg                          <  1:24921/251‑1 (MQ=255)
                                                                                                                                                                                                                          tttATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCCGCGATGTTCCCGCCGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCgggg                                                                                                               >  1:200888/1‑163 (MQ=255)
                                                                                                                                                                                                                          tttATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTg                                                                                                               <  2:200888/165‑1 (MQ=255)
                                                                                                                                                                                                                                  gATTAACATGAACGTCATCGGTCTGATTTTGGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTggagg                <  2:104737/251‑4 (MQ=255)
                                                                                                                                                                                                                                               gTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATTCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTTGGTGCTGGTTGCCGGTCTGACGCC‑CACGCTGGCGGCCTGGCTGGGCGGAAGCTCGCAGAATCTGATGATGCCTGCCTATTTCCTGATGGTAGTGGCgggtggtggggtttt   >  1:17265/1‑238 (MQ=255)
                                                                                                                                                                                                                                                 cATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGGTGCCTGCCTATTACCTGATGGTAGTGGCgggggggggtgttttt  >  1:417536/1‑236 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
CACCAACGTGACTTACTACATGTTGCTGACCTATATGCCGAGTTATTTGTCGCATAACCTGCATTACTCCGAAGACCACGGGGTGCTGATTATTATCGCCATTATGATCGGTATGCTGTTTGTCCAGCCGGTG‑ATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATC  >  NC_000913/4331311‑4331800

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