Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,185,1990TG22.2% 132.7 / 10.6 72T35P (ACC→CCC) thiMhydroxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base T (12/44);  new base G (16/0);  total (28/44)
Fisher's exact test for biased strand distribution p-value = 7.39e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAAGTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCG  >  NC_000913/2184953‑2185448
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
aaaGTTCATGACAAAAAGGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCGTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTGTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTTGCGGTAAAGGTTTGc                                                                                                                                                                                                                                                           <  2:438270/251‑1 (MQ=255)
   gTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacc                                                                                                                                                                                                                                                       <  1:403869/251‑1 (MQ=255)
             aaaaTGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCAtttg                                                                                                                                                                                                                                             >  1:239868/1‑249 (MQ=255)
                        gATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGGGTTGccccccctccattggcatggagggg                                                                                                                                                                                                                                   >  1:612/1‑227 (MQ=255)
                                   ccaccccccGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGa                                                                                                                                                                                                                      <  2:25724/250‑1 (MQ=255)
                                         cccGCTACTGGATAAAGCGTCCAGGGTGTTTGAGAGCTTTGTGCTTGCTCAACGGCAGCACGCATCGCCGGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTGGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATGGGTCATGCAGTGCACAAGAGGGGAATGTTg                                                                                                                                                                                                                <  2:359431/250‑1 (MQ=255)
                                                             ccAGGTTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCAGTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGt                                                                                                                                                                                           <  1:475467/251‑1 (MQ=255)
                                                               aGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTgc                                                                                                                                                                                         <  1:461976/251‑1 (MQ=255)
                                                                 ggTGTTTGAGCGCTTTTTGCTTGCTCAAAGGCAGCACGCATCGCCTGAGCGCGAGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCATCTTCGGTTTCGATAACCATCGCTGGCGATGCCCCGAGCGCCAGCAAGGTATTGGCGGTCAAGGATTTGTCCCACTCATTGGTCAGTCAGGGCACACGGGGGGAATGTGGGTGAGAAAggcggaacacttgagc                                                                                                                                                                                       >  2:486761/1‑237 (MQ=255)
                                                                  gtgtTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCATGCAGTGCACCAGAGGGGGATGTTGGGGGAAAAGGTGTAACGCGTgcgcc                                                                                                                                                                                      >  2:15064/1‑250 (MQ=255)
                                                                    gtTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGCAAAGGGTTGCACCCCCTCATTTGTCATGCAGTGCACAAGAGGGGGATGTTGGTGGAAAAGGTGGAAAGCGTgcgc                                                                                                                                                                                       >  1:54003/1‑248 (MQ=255)
                                                                       tGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGAtt                                                                                                                                                                                  <  2:239868/250‑1 (MQ=255)
                                                                        gagaGCTTTTTGCTTGCTCAACGGCCGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGGGTTGTACCCCATCATTTGTCATGCCGTGGACAAGAGGGGGATGTTTGTGGAAAAAGGGTCACCCGCGCGCAGATTgc                                                                                                                                                                                >  1:427198/1‑251 (MQ=255)
                                                                              tttttGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGCAAAGGGTTGCACCACATCATTGGGCATGCAATACACCAGGGGGGGAGGTTGGTGGAAAAGGGGTAAAGCGGGCGGAGATTTCGCTGa                                                                                                                                                                           >  1:478650/1‑250 (MQ=255)
                                                                                   gcttgctCAACGGCAGAACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTCCCCACCTCATTGGTCATGCCGTGCCCAAAAGGGGGATGTTTGTGAAAAAGGGGTAACGCGCGCGCAGGTTGCCCTGAAccc                                                                                                                                                                       >  1:371439/1‑249 (MQ=255)
                                                                                           aaCGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTGACCACATCATTTGGCATGCAGTGCACAAGAGGGGGATGTTGGTGGAAAAGGGGTAAAGCGTGGGGCGAATGCGCTGAACCCAGCAgggcg                                                                                                                                                              >  1:265783/1‑247 (MQ=255)
                                                                                              ggCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACtt                                                                                                                                                          <  1:144746/251‑1 (MQ=255)
                                                                                              ggCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACtt                                                                                                                                                          <  2:127387/251‑1 (MQ=255)
                                                                                              ggCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACtt                                                                                                                                                          <  1:268001/251‑1 (MQ=255)
                                                                                              ggCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGCACCACATCATTGGTTATGCCGTGCACAAGAGGGGGATGGTGGGGAAAAAAGGGTAACGCGGGCGCCGATTTGGGTGAACCCAGCAgggcggctt                                                                                                                                                          >  1:376311/1‑244 (MQ=255)
                                                                                              ggCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGAAAAGGGTTTGACCCCAACAATGGGCATGCAGTGCACAAGGGGGGGATGTTGGTGGAAAAGGGGGAACGCGTGCGCGGAGTGCGGTGaaacccgacgggtggatt                                                                                                                                                          >  1:187945/1‑235 (MQ=255)
                                                                                                 agcaCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCa                                                                                                                                                       <  2:155894/251‑1 (MQ=255)
                                                                                                        aTCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAAc                                                                                                                                                                                               <  1:346189/204‑1 (MQ=255)
                                                                                                        aTCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTg                                                                                                                                                <  1:234982/251‑1 (MQ=255)
                                                                                                          cGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAAGCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTATCGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCt                                                                                                                                              <  2:52786/251‑1 (MQ=255)
                                                                                                            ccggAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTcc                                                                                                                                            <  2:380511/248‑1 (MQ=255)
                                                                                                             cTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGGTTTGCACCACATCATTGGTCATGCAGTGCCCAAGAGGGGGATGTTGGTGGAAAAGGGGTAACGCGGGCGCAGATTGCGCTGAACCCCGCCGGGCGACTTGCATAGTTTGCTccc                                                                                                                                           >  1:32488/1‑250 (MQ=255)
                                                                                                                      tGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGCAAAAGTTTGTCCCCCCTCATTGGGAACGCGGGGCACAAGAGGGGGATGTTGGTGGAAAAAGGGTAAAGCGTGCGCCGATTGCGCTGAACCCCGCACGTCGACTTGGATTGttttcttctcccaatacc                                                                                                                                  >  1:168796/1‑236 (MQ=255)
                                                                                                                       ggCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGt                                                                                                                                 <  1:524626/251‑1 (MQ=255)
                                                                                                                       ggCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGt                                                                                                                                 <  1:263721/251‑1 (MQ=255)
                                                                                                                                  ggtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATGGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAAt                                                                                                                      <  2:612/250‑1 (MQ=255)
                                                                                                                                   gtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAAt                                                                                                                      <  1:9363/250‑1 (MQ=255)
                                                                                                                                   gtgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACAGCATTGGTCATGCAGTGCACAAGAGGGTAATTTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATg                                                                                                                     <  1:40718/251‑1 (MQ=255)
                                                                                                                                       cAACGTTAATCAACAAGGCGCTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTTCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCt                                                                                                                 <  1:519225/251‑1 (MQ=255)
                                                                                                                                           gTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTTCACCACCTCCTTGGGCATGCCGTGCACAAGAGGGGGAAGTTGGGGGAAAAAGGGTTACGCGTGCGCAGATTTCGCTGAAACCCGCAGGGCGACTTGCCTAGTTTGCTCCTGCCATCAAGTGACGAAGCagatcgggaggg                                                                                                             >  2:156388/1‑240 (MQ=255)
                                                                                                                                           gTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGc                                                                                                                         <  1:156388/239‑1 (MQ=255)
                                                                                                                                              aaTCAACAAGGCACTGGCGTTAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGCACCACATCATTTGTCCTTCCGTTCACAAAAGGGGGATGGTGGGGGAAAAGGTGTACAACGTTGGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATCGTTTGCTCCTGCCCTCACGTGAAGACGCCATGGCCtggtgttc                                                                                                         >  2:309921/1‑248 (MQ=255)
                                                                                                                                              aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTTGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAGGGTTTGCACCACATCATTGTTCATGCAGTGCACAAGAGGGGAAGGTTGGTGAAAAAGGTGGAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGGTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTggtggt                                                                                                          <  1:315120/250‑1 (MQ=255)
                                                                                                                                               aTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACCTCCTTGGTCAAGCCGTGCACAAGAGGGGGATGTTGGGGGAAAAAGTGGAACACGTGCGCCGAATGCGCTGAACCCAGCAGGGCGACTTGCATAGTTTTCTCCTGCCCTAAACTGAAGAAGCAAAGACCTtgtggttc                                                                                                        >  1:494346/1‑249 (MQ=255)
                                                                                                                                               aTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCCGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGTTTGCCACCCATCATTCGTCCTGTCGTGCCAACGAGGGGGATGTTGGTGAAAAAAGGGTAACGCGGGCCCAGATTTCGCTGAACCCAGCACGGCGCCTTGCCCAGGTTGTTCCTGCCCTCACATGAAGAAGCAATGACCttgtgggtc                                                                                                        >  1:80041/1‑248 (MQ=255)
                                                                                                                                                 caacaaGGCCCTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGCCCCCCACCCTTTGTCATGCAGTGCCCAAGAGGGGGATGTTGGAGAAAAAAGAGGAACGCCTGCGCCGATTGCGCTGGACCCCGCAGGGCGACCTGCATAGTTTGCTCCTGCCATTACGTGAAGAAGCAAAGACCTTggtgggcgg                                                                                                       >  1:421436/1‑247 (MQ=255)
                                                                                                                                                 caacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCg                                                                                                       <  1:447917/250‑1 (MQ=255)
                                                                                                                                                 caacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGCACCACCTCCTTGGTCATGCAGTGCACAAGAGGGGGATGTTGGGGAAAAAGGTGTAAAGCGTGCGCAGGTTGCGCTGGACCCAACAGGGCGACTTGCCTAGTTTGCTCCTTCCATAACGTGAAGAAGCAATGACCTTGTGGTCCgg                                                                                                      >  1:127154/1‑250 (MQ=255)
                                                                                                                                                    caaGGCACTTGCGATAGCCGCAATCTGACTGGCCTCGTCGGTTGCGATAACCATCGCTTGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAg                                                                                                                                                    <  2:145263/203‑1 (MQ=255)
                                                                                                                                                     aaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTTCACCACATCATTGGGCACGCAGTGCACAAGAGGGGGATGTTGGGGGAAAAAGGGTTACGCGGGCGCCGATTGCGCTGAACCCCGCCGGGCGACTTGCCTAGGTTGCTCCTGCCCTTACGTGAAGAAGCAAACGCCTTGTGGTCCGGGACt                                                                                                  >  1:194319/1‑251 (MQ=255)
                                                                                                                                                              gCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGc                                                                                         <  1:289541/251‑1 (MQ=255)
                                                                                                                                                               cGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGCTGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTTTACAACACCATTCGGCATTCAACGCACCCGCGGGGGGTGTTGGTGGAAAAAGGGTAACGCGTGCGCACATTGCGCTGAACCCCACAGGGCGACTTTCATAGTTTGCTCCTGCCATAAAGTGAAGAAGCAATGACCTGGGGGGTCGGGACTTCCCTAcctc                                                                                        >  1:224855/1‑248 (MQ=255)
                                                                                                                                                               cGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCt                                                                                        <  1:15064/251‑1 (MQ=255)
                                                                                                                                                                gATAGCCGCAACCTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTGGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTg                                                                                       <  2:421436/251‑1 (MQ=255)
                                                                                                                                                                         aaaCTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATcc                                                                              <  1:475090/251‑1 (MQ=255)
                                                                                                                                                                              gactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGTTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGTCCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATc                                                                         <  1:433005/251‑1 (MQ=255)
                                                                                                                                                                              gactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTGTGCTCCTTCCATAACGTGAAGAAGCAAGGGCCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATc                                                                         <  1:486761/251‑1 (MQ=255)
                                                                                                                                                                                    ccTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATATTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGa                                                                   <  2:494346/251‑1 (MQ=255)
                                                                                                                                                                                                            cTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTATGCAACACATCATTGGTCATGCAGTGCACAAGAGGGGGAAGTTGGGGAAAAAGGGGTAACGCGGGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCCTAAAGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACCCTGGCATTATCCAGATCAGGTGATACGGGGATTTCTCAGCCTTCACg                                           >  1:42028/1‑251 (MQ=255)
                                                                                                                                                                                                            cTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCCCATCATTGGTCATGCCGTGCCCAAGAGGGGAATGTTTGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCCTAGTTTGCTCCTGCCCTAACGTGAAGAAACAAAGACCTGGTGGTCCGGGACTTCCCTACGCTGGCCTTATTCAGATCAGGGGGTACGGGTATTTTTCAGCCTTcccc                                           >  1:183619/1‑248 (MQ=255)
                                                                                                                                                                                                             tGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCCCATCCTTGGTCCTGCCGTGCACAAGAGGGGGATGGTGGTGGAAAAAGTGTAACGCGTGCGCAGATTTCGCTGAAACCCGCCGGTCGACTTGCATAGGTTGCTCCTGCCCTAACGTGAAAAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTCTCCCGATCAGGGGCTACGGGTATTTCTCACCCTTCACGc                                          >  2:402398/1‑251 (MQ=255)
                                                                                                                                                                                                                 gATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGGCCCACATCATTTGTCATGCCGTGGCCAACAGGGGGATGTTTGGGGAAACGGTGTGAAGCGTGCGCCGATTTCGCCGGAACCAGCCGGGCGCCTTGCATAATTTTATCCAGCCATAACGTGGAGAAGCAATGAACTGGTGGTCCCCGGCTTCCCTAAGCTTGCCTTATACAGGACAGGGGATACGGGGATTTTTCAGCCATCACGCAgc                                       >  1:494561/1‑249 (MQ=255)
                                                                                                                                                                                                                 gATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTTGTCATGCAGTGCACAAAGGGGGGATGTTTGGGGAAAAGGGGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGGCGACTTGCATAGTTTGCTCCTGCCATAAAGTGAAGAAGCAATGACCCGGGGGGCCGTGACTTCCCCTCGCTGGCATTATCCAGATCAGGGGATACGGGGATTTCTCAGCCTTCCCGCAGa                                       >  1:510460/1‑250 (MQ=255)
                                                                                                                                                                                                                    gccccGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTGGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAGCCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAgg                                    <  2:168796/247‑1 (MQ=255)
                                                                                                                                                                                                                    gCACCGAGCGCCAGCAAGGTATTGGCTGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAgg                                    <  2:205714/250‑1 (MQ=255)
                                                                                                                                                                                                                    gCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAgg                                    <  1:30557/250‑1 (MQ=255)
                                                                                                                                                                                                                      aCCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCAGAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCa                                 <  1:253624/251‑1 (MQ=255)
                                                                                                                                                                                                                          aGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGCACCACATCATTGGTCATGCAATGCACCAGAGGGGGCTGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGGACCCCGCCGGTCGACTTGCATAGTTTGCTCCTGCCCTAACGTGAAGAAGCAATGAACTGGTGGGCCGTGACTTCCCTACCCTGGCATTATCCAGATAAGGGGGTCAGGGGATTTTTTAGCCTTCACGCCGAAAGGCAcccc                             >  2:388098/1‑251 (MQ=255)
                                                                                                                                                                                                                             ggcAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCATTCACGCAGAAGGGCACCCCGAg                          <  2:92219/249‑1 (MQ=255)
                                                                                                                                                                                                                              ccAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGt                         <  1:338279/251‑1 (MQ=255)
                                                                                                                                                                                                                               cggCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTc                        <  2:54003/249‑1 (MQ=255)
                                                                                                                                                                                                                                        aTTGGCGGTAAATGTTTGCACCACATCATTGGCCATGCAGGGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACTCGTGCGCAGATTGCGCGGAACCCAGCAGGTCTTCGTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTTGTCCGTGACTTCTTTATGCTGGCATTATCCAGATCAGGTGATACGGGTATTGCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTGTGGtt                <  2:167060/250‑1 (MQ=255)
                                                                                                                                                                                                                                             cGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACTGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCgatg          <  1:44496/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ggTAAAGGTTTGCACCACATCATTCGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCCTTATCCAGATCAGTTGATACGGGTATTTCTCAGCCTTCACGCATAAGGGCACCCCGAGTCGTTTGGTTGCgatgt         >  2:257502/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   aGGTTGGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCATGCGCAGATTGCGCGGAACCCAGCAGGGCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGTCTTCCCGACGCTGGCATTAGCCAGATCAGGTGATACGGGTATTTCTCAACCTTCACGCAGAAGGGCACACCGAGTCCTTTGGTTGCGATGATAAAt    <  2:262402/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     ggTTGGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATc   <  2:265783/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      gTTTGCACCAAAGCATGGGTCATGCAGGGCACAAGAGGGGAAGGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCg  <  2:371439/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
AAAGTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATGATAAATCG  >  NC_000913/2184953‑2185448

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: