Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,100,3570TC31.8% 49.2 / 32.6 45V21A (GTG→GCG) ycdZinner membrane protein YcdZ
Reads supporting (aligned to +/- strand):  ref base T (11/19);  new base C (6/8);  total (17/28)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGC  >  NC_000913/1100111‑1100590
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
ccTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGCATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCt                                                                                                                                                                                                                                       >  1:373192/1‑251 (MQ=255)
catttCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCt                                                                                                                                                                                                                                       >  1:373190/3‑251 (MQ=255)
 cTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTg                                                                                                                                                                                                                                      >  1:343498/1‑251 (MQ=255)
       ggCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTc                                                                                                                                                                                                                                  <  2:165339/249‑1 (MQ=255)
            aCCATGGCACCGCTAACCCGCGATGCCATTAGTGCCATGTGGGACGAGATGGCAGCAGAATCTGAAGAATAAGTGTGATCTACGTCACTAATAACTGCAACGGATAATTTGTTGTTGAATAAAATGTGTGCTCGATCTCCTTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATCAACATACTTCTCTCTATTGCAATCACAACAGGAATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGc                                                                                                                                                                                                                           >  1:218665/1‑251 (MQ=255)
            aCCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGc                                                                                                                                                                                                                           <  2:89233/251‑1 (MQ=255)
             ccATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACATCAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGCt                                                                                                                                                                                                                          <  2:325205/251‑1 (MQ=255)
                                    gCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCt                                                                                                                                                                                                    <  1:364425/250‑1 (MQ=255)
                                      cATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTggg                                                                                                                                                                                                 >  1:173079/1‑251 (MQ=255)
                                                          aGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCATCGGATAATTTGTTGTTGCATAAAAGGTGTGCTCGATCTCAGTCATGGCCTCGTTTTCTGCTTTCAGGCGCGGCATGAACATACTTCTCTCTATTGCAATCCCAATAGGCATTCTCTCCGGGATCTGGTGATGGGTGGCTGTGTCCCTGGGCTTACTAAGCTGGGCGGTCTTCCGGGGCTGTACGGCCTACTTtgcct                                                                                                                                                                             <  2:191529/251‑1 (MQ=255)
                                                             tGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTtgcctgcc                                                                                                                                                                          >  1:147468/1‑251 (MQ=255)
                                                                   aagaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACtt                                                                                                                                                                                  <  1:120449/237‑1 (MQ=255)
                                                                   aagaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACtt                                                                                                                                                                                  >  2:120449/1‑237 (MQ=255)
                                                                    agaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAgg                                                                                                                                                                   <  2:49509/251‑1 (MQ=255)
                                                                      aTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTg                                                                                                                                                                 <  2:228491/251‑1 (MQ=255)
                                                                                    gtgtGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTGCCTGGGCTGTACGGCCTACGTTGCCTGCCCGCGATGTGGTCTGAAAGGGCTg                                                                                                                                                   <  2:216587/251‑1 (MQ=255)
                                                                                         aTCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGCACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGAt                                                                                                                                              >  2:141634/1‑251 (MQ=255)
                                                                                                 cACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTTGGCTGAAAGGGCTGGCGATCTCCGCTg                                                                                                                                      <  2:281092/251‑1 (MQ=255)
                                                                                                            gCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTGCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGTTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGc                                                                                                                             <  2:186355/249‑1 (MQ=255)
                                                                                                               ccGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGt                                                                                                                        <  2:110507/250‑1 (MQ=255)
                                                                                                                     aaTTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCgtg                                                                                                                  >  1:234924/1‑251 (MQ=255)
                                                                                                                          gttgttgCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGgtgt                                                                                                              <  1:188233/250‑1 (MQ=255)
                                                                                                                          gttgttgCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGgtgt                                                                                                              <  1:321424/250‑1 (MQ=255)
                                                                                                                                  aTAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATgg                                                                                                     >  1:199512/1‑251 (MQ=255)
                                                                                                                                             gCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTAc                                                                                          >  1:210010/1‑251 (MQ=255)
                                                                                                                                                             ggCCGCGTTTTCTGCTATCATGCGCGGCATGACCATCCTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGCTCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGG‑CGGTGGTCAGTAttttccggtgtgcaccgggccc                                                                         >  2:16120/1‑232 (MQ=255)
                                                                                                                                                                          gCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGaa                                                               <  2:33240/249‑1 (MQ=255)
                                                                                                                                                                                aTGCGCGGCATGAAAATACTTCTCTCTATTGCAATCTCAACAGGCATTCTCACCGGGAACTGGGAATGGGGTGCTGTTTCCCGGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCTCTGCAACCCTGCTTTGTGGCGTGGTGTGGGCGATGGTCATTAGTTACGGTAGTGCACTGGCTCCACATCTGGAAACTCTCgg                                                       <  2:282402/251‑1 (MQ=255)
                                                                                                                                                                                  gcgcGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGt                                                      <  1:353816/250‑1 (MQ=255)
                                                                                                                                                                                   cgcgGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTa                                                    <  1:228148/251‑1 (MQ=255)
                                                                                                                                                                                        cATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCa                                               <  1:173344/251‑1 (MQ=255)
                                                                                                                                                                                          tGAACATACTTCTCTCGATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGTTCATTATTTACGGTAGTGCACTGGCAACACATCTGGAAATGCTCAGTTATGTCAt                                              <  2:47793/250‑1 (MQ=255)
                                                                                                                                                                                             aCATACTTCTCTCTATTGCACTCAAAACAGGCATTCCCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGATGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGAGTGGGGTGGGCGATGGTCATTCTTTACGGTAGTGCACTGGCACCACATCTGGAAATTCCCGGTTATGTCATAAcc                                          >  2:235100/1‑251 (MQ=255)
                                                                                                                                                                                                  cTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCCCCACATCTGGAAATTCTCGGTTATGTCATAACCGGTAt                                     >  1:261125/1‑251 (MQ=255)
                                                                                                                                                                                                     cgctctATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGt                                  <  2:261125/249‑1 (MQ=255)
                                                                                                                                                                                                         ctATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTAATTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTc                                 <  2:213868/248‑1 (MQ=255)
                                                                                                                                                                                                         ctATTGAAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTcg                                <  1:284286/249‑1 (MQ=255)
                                                                                                                                                                                                             tGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGttt                           <  1:62019/250‑1 (MQ=255)
                                                                                                                                                                                                                aaTCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTg                        >  1:202882/1‑250 (MQ=255)
                                                                                                                                                                                                                       ccaGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGTCGGTCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTAt                <  2:62685/250‑1 (MQ=255)
                                                                                                                                                                                                                         aGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTAtt               <  2:110592/250‑1 (MQ=255)
                                                                                                                                                                                                                         aGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTAtt               >  1:110592/1‑250 (MQ=255)
                                                                                                                                                                                                                             aTTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGc          <  2:140099/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tcCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTCTTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAAc     >  1:336604/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cGGTATCTGGGGATGGGCGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAgc  <  1:141634/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         
CCTTTCTGGCGCACCATGGCACCGCTAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAGGCATTCTCTCCGGTATCTGGGGATGGGTGGCTGTTTCCCTGGGCTTACTAAGCTGGGCGGGCTTCCTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGC  >  NC_000913/1100111‑1100590

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