Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,930,7180TG28.2% 186.5 / 14.0 79T126P (ACC→CCC) ravAhexameric AAA+ MoxR family ATPase, putative molecular chaperone
Reads supporting (aligned to +/- strand):  ref base T (9/47);  new base G (22/0);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 1.42e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGATTGTCGTTTTCATCCTGTTGACTGGTCAGCATGGAGCGGAAATTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCG  >  NC_000913/3930495‑3930954
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
ggATTGTCGTTTTCATCCTGTTGACTGGTCAGCATGGAGCGGAAATTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCgg                                                                                                                                                                                                                                        >  1:1435870/1‑239 (MQ=255)
ggATTGTCGTTTTCATCCTGTTGACTGGTCAGCATGGAGCGGAAATTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATCGCCgg                                                                                                                                                                                                                                        <  2:1711808/239‑1 (MQ=255)
     gTCGTTTTCATCCTGTTGACTGGTCAGCATGGAGCGGAAATTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCg                                                                                                                                                                                                                                  <  1:1629489/239‑1 (MQ=255)
               tCCTGTTGACCTATCAGCATGAAGCGGAAATTCGCTTTAACCTCCCCTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACGGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGAt                                                                                                                                                                                                                        <  2:1057757/239‑1 (MQ=255)
                ccTGTTGACTGGTCAGCATGGAGCGGAAATTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGGATTAAGAATTTCCGGGCCCGCTTTCCCGa                                                                                                                                                                                                                         >  1:212294/1‑237 (MQ=255)
                              aGCATGGAGCGGAAATTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCACGGGATTAAGAATTTCCGGGCCCGCTTTCCCGCTCTCATCCACAAAGAc                                                                                                                                                                                                        >  1:1537031/1‑239 (MQ=255)
                                     aGCGGAAATTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAAAATTGCCGGGGCCGCTTTTCAGAACTCCTTCCGAAAAACCAAttt                                                                                                                                                                                                   >  1:604577/1‑238 (MQ=255)
                                      gCGGAAATTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGGTAGCAAGGGATTAAGAATTGTCCGGGCCGCCTTCCCGCTCTCATCCCAAAAGACAAta                                                                                                                                                                                                    >  2:1461594/1‑234 (MQ=255)
                                          aaaTTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGGATTAAGAAATTCCGGGCCCGCTCTCTCGATCTCTTCCAGAGAAACAAattttcttt                                                                                                                                                                                             >  2:2141065/1‑232 (MQ=255)
                                            aTTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGGGTGCAAGGGATTAAGAAATTTCGGGGCCCCTTTCCAGATCTCATCCAGAAAGACAAATTTTGCTTccc                                                                                                                                                                                         >  1:2017091/1‑238 (MQ=255)
                                               cGCTGTATCCTGGCCTTCATCTCCACACAGACGAACCCGACGGGGTTCATATAACGTCTCAAGCCTCTGGTCGGCTGCCGGCGGCTCTTTGGAGGCCGCCCCCCGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGGGGTGAGCAAGGTATTAAGAATTCCCGGGCCCGCTTTCCAGATCTCACCCAGAAAGACAATTTCTGCTTCCgg                                                                                                                                                                                        <  2:2555237/237‑1 (MQ=255)
                                                      accTGCACTTTACCTAACCCCAGCCGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCg                                                                                                                                                                                         <  2:2521582/228‑1 (MQ=255)
                                                       ccTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCg                                                                                                                                                                            <  1:135589/239‑1 (MQ=255)
                                                       ccTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCg                                                                                                                                                                            <  2:143381/239‑1 (MQ=255)
                                                           cACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGGATTTATAAATTCCCGGCCCGCTTTCCAGATCTCCCCCAGAAAGACAAATTCTTCCTCCCGCACGGAACCg                                                                                                                                                                            >  2:804879/1‑235 (MQ=255)
                                                              tttATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTTTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTa                                                                                                                                                                    <  1:2368158/239‑1 (MQ=255)
                                                              tttATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTa                                                                                                                                                                    <  1:1062567/239‑1 (MQ=255)
                                                              tttATCTAACCACAGACGAAGCAGCCTGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGTAGCGCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAGCTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTa                                                                                                                                                                    <  2:617192/239‑1 (MQ=255)
                                                               ttATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTaa                                                                                                                                                                   <  1:2620860/239‑1 (MQ=255)
                                                               ttATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTaa                                                                                                                                                                   <  2:1898916/239‑1 (MQ=255)
                                                                tATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTTGGCAATGGATTAAAAATTGACCGGCCCCCTCTCCAAATGTCATCCCGAAAGACACATTCTGCCTCCCGCCGGGTAACCACGGGtcta                                                                                                                                                                   >  2:1062567/1‑236 (MQ=255)
                                                                tATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGccc                                                                                                                                                                                                                                    >  1:571042/1‑179 (MQ=255)
                                                                tATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGccc                                                                                                                                                                                                                                    <  2:571042/179‑1 (MQ=255)
                                                                tATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTa                                                                                                                                                                    <  2:442905/237‑1 (MQ=255)
                                                                tATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTa                                                                                                                                                                    >  1:442905/1‑237 (MQ=255)
                                                                tATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGGATTAAGAATTTCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGGAAcc                                                                                                                                                                             >  1:2521582/1‑229 (MQ=255)
                                                                tATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGGATTAAGAATTGCCGGGGCCGCTTTCCAGATCTCATCCCGAAAGACAATTTTTGCTTCCGGCCGGTAACCGCTGGGTaa                                                                                                                                                                   >  1:1717666/1‑238 (MQ=255)
                                                                tATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTCCGGAACTGGCGCTCGTTAATGGCGGTCAGCAAGCTATTAAGAAATTCCGGGGCCGCGTTCCAGGTGTCCTCCCATAACGCAATTTTTTCTGCTCGCAGGTAACCGcgtgcttc                                                                                                                                                                    >  2:2553719/1‑231 (MQ=255)
                                                                      aCCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGGGCAAGGTATTAAGAATTGCCGGGCCCCCTTTCCAGATCTCATCCCGAAAAACAATTTTTGCTTCCGGCAGGGAACCGCTGGGTAAACGTTCa                                                                                                                                                            >  1:689715/1‑239 (MQ=255)
                                                                          caGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGGGCAAGGGATTAAGAATTGTCGGGCCCGCCTTCCAGATCTCATCCAAGAAGAAAAATTCTTCCTCCTGCCGGGAACCCCCGGCTAAAAGTTCATat                                                                                                                                                        >  2:1443656/1‑237 (MQ=255)
                                                                                 aTCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACTGCTGGTTAAACGTTCATAGCGCCCTTc                                                                                                                                                <  1:334746/239‑1 (MQ=255)
                                                                                       aTGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTTCAGATCTCATCCAGAAAGACAATTTTTGCTTCCGGCCGGTAACTGCCGGGTTAACGTTCCTAGCGCCCCTCATCttt                                                                                                                                          >  2:334746/1‑239 (MQ=255)
                                                                                        tGCGGTCATATAACGCTTCCAGGCTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCtttt                                                                                                                                         <  1:2375054/239‑1 (MQ=255)
                                                                                        tGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCtttt                                                                                                                                         <  1:837985/239‑1 (MQ=255)
                                                                                           ggTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGGTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAg                                                                                                                                      <  1:1686294/238‑1 (MQ=255)
                                                                                             tCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACg                                                                                                                                                                <  1:1362429/212‑1 (MQ=255)
                                                                                              cATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAgcgc                                                                                                                                   <  1:878549/238‑1 (MQ=255)
                                                                                              cATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGcc                                                                                                                                  <  1:2141065/239‑1 (MQ=255)
                                                                                               atatAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGt                                                                                                                                                                      <  1:2756277/204‑1 (MQ=255)
                                                                                               atatAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCt                                                                                                                                 <  1:1872217/239‑1 (MQ=255)
                                                                                               atatAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCt                                                                                                                                 <  1:495407/239‑1 (MQ=255)
                                                                                                tataACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTcc                                                                                                                                                                                                                            <  1:52413/154‑1 (MQ=255)
                                                                                                      gattCCAGGCTGGTGTATGCTTACGTCAGCTCGTTGGGAGCCGGCACCAGCAGACGCATCGGGATTGTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCACGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGa                                                                                                                          <  2:1441892/237‑1 (MQ=255)
                                                                                                       cTTCCAGATTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTTTAATAAGTATTTAAAACTGACCGGCCCGCGTCCCCGATCTCATTCAGAAAGACAAATTTTTTTCTCGGCAGGGTACCGCTCGGTTAACGTTCCTAAGGCCCCTTCACTTTTCGCGCCTGacatgaa                                                                                                                         >  2:497533/1‑233 (MQ=255)
                                                                                                        ttCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGc                                                                                                                                                                                                                                           <  1:2535321/132‑1 (MQ=255)
                                                                                                         tCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGTAAGGGATTAAGAATTGCCGGGCCCGCCTTCCAGATCTCATCCAGAAAGACAATTTTTGTTTCCGCCCGGAGATCGCCGGTTCTAAATTCATAACGTCCCTCCTCTTATCTGGCCTCAATAgcata                                                                                                                       >  2:541969/1‑235 (MQ=255)
                                                                                                          ccAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGGATTAAGAATTTCCGGGCCCCCTTTTCAGATCCCATc                                                                                                                                                                                                                  >  2:52413/1‑154 (MQ=255)
                                                                                                          ccAGACTACTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCAACCAGCAGACGCATCGGGATTTTTTCCACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGCAAGCGAAGGATTTAAAAATGAGGAGCCCGGTTTTCCAGTCTCCACCCGAAAGACAAATTCTTACTCCGGACAGGTACAGCTGGGTAAAAGTTCACAACGCCCTTCATCTTTTTGAGCCCGAAATGaaaca                                                                                                                     >  2:1354248/1‑237 (MQ=255)
                                                                                                                  ctgTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGGAAGGGATTTAGAAATGCCGGGCCCGcc                                                                                                                                                                                                                                >  2:2535321/1‑131 (MQ=255)
                                                                                                                     tCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGAATTAAAAATTACCGGGCCCGCTTTTCAGAACTCATACCGAGAAACAATTTCTGCTTCCGGGAGGTAACCCCTCGGTACACGATCATAGCGGCCTTCATCTTTTTGCGCCTCAATAGATATGGGGACCaaaaa                                                                                                           >  2:231175/1‑239 (MQ=255)
                                                                                                                          ttCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAActtctt                                                                                                     <  1:804879/239‑1 (MQ=255)
                                                                                                                                   cTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTTAAAAATTCCGGGACCGCCTTCCAGATCTCCTCCAGGAAGACAATTTTTTCTTCCCGCCGGGAACCGCTGGGTTAACCGTTATAAGCCCCTTCATCGTTTATGGGCTTGATAGAAAAGGGACCAAAAACTTCCTCCTGCGCgg                                                                                            >  2:432652/1‑239 (MQ=255)
                                                                                                                                    tCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGt                                                                                              <  1:231175/236‑1 (MQ=255)
                                                                                                                                    tCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGGGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTTGGCAAGGGAATAAAAATTTCCGCGCCCGCTTTTCAGGTGTCATCCAGGAAGAAAATATCTTCTTGTCGGAGGGGACCGCGGATGACACGTTTATAACGGCCTTCATCTTTTCGCGCCGGAATAGAAAAAGGCGCAAAAACAACGTCCGGCGCGga                                                                                           >  2:1933735/1‑239 (MQ=255)
                                                                                                                                        tGGAAGCCGACACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACAGTTGCGGAACTGGCGCTCGTTAAAGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAg                                                                                       <  2:2687812/239‑1 (MQ=255)
                                                                                                                                              ccgccACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCggg                                                                                   <  1:428986/237‑1 (MQ=255)
                                                                                                                                              ccgccACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCggg                                                                                   <  1:432652/237‑1 (MQ=255)
                                                                                                                                                     cagcagACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGcc                                                                                                                                  <  1:68030/184‑1 (MQ=255)
                                                                                                                                                        cagACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGGGTTAAAAATTGACGGGCCCCCCTTCCAAGTCTCATTCCGAAAGAAAATTTTTTCCTCCTGCAGGTAAACCCCTGGTAAAAGGTTAAAACGACCATCGTTTTTTAGGCGCCTAAATGAAAAGGGGACAAAAAACTCTTCCCGCGGCGTGAAGGGGGtcgttcagaatt                                                                       >  2:393074/1‑229 (MQ=255)
                                                                                                                                                            cGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCa                                                                     <  1:1443656/237‑1 (MQ=255)
                                                                                                                                                               aTCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGGAAGGGATTAAGAATTTCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTTTGCTTCCGGCCGGTAACCGCCGGGTAAACGATCATAGCGCCCTTCATCTTTTTGCGCCCGAATATaaa                                                                                                                       >  2:68030/1‑184 (MQ=255)
                                                                                                                                                                   ggATTGTTTCAGCGTGTGGAGCGTTGCTGAACTGGCGCTGTTTAATGGCGGTGAGCAAGGTATTAAGAAGTGCCGGGCCCGCGTGCCAGATCTCATCCAGAAAGACAATTTCTGCGTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAAcgagcg                                                             <  2:848401/238‑1 (MQ=255)
                                                                                                                                                                   ggATGTTATCTACGTGTGCATCATTGCGGAACTGGCTCGCGATAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGACCCTGTTTCCAGAACTCAACCCGAAAGACACTTTCTACTTCCGGCAGGTAACCGCGGGTTAAACGTTCATAGCGCCCTTCAACTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAACGGGTCATCAGATTTTCAAgcgggcgg                                                            <  2:604577/239‑1 (MQ=255)
                                                                                                                                                                    gATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCgg                                                            <  1:2519559/238‑1 (MQ=255)
                                                                                                                                                                                gtgtGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGTCGGGc                                                                                                                                                                                                                                      <  2:1704670/65‑1 (MQ=255)
                                                                                                                                                                                gtgtGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTaa                                                                                                                                                                   <  2:1834634/126‑1 (MQ=255)
                                                                                                                                                                                   tgCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGGATTAAGGATTGCCGGGCCCGCCTTTCAAATGGCATCCAAAAAGACAAATTTTGCTTTCGGCAGGGAACCGCCGGTGTAACCTTCATAGCGCCCTTTATCTTTTTCGCCCCGAAAAGAAATGGGAGCAAAAACTTCTTTCCGCGCGGAGAACCGGGGCACAAGACATTCACAAGCGCGCGCCATCTGAGAGGGGa                                            >  2:2529733/1‑239 (MQ=255)
                                                                                                                                                                                          ttGCGGAACTGGCGCTCGTTAATGGCGGTGAGCCAGGGAATAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTTTGCTTCCGGCAGGTAACCGCTGGTTAAAACGTCATAGCGCCCTTCCTCTTTTATCGCCTTCATAGAAAAGGGACAAAAAAATTCTTCCACCGTGGGTCAGTGGGGCATCAAAAATTTAAAACCGCGGGCATTCTCAAAggtggatattta                                     >  2:2084997/1‑229 (MQ=255)
                                                                                                                                                                                          ttGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGGATTTAGAATTTCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTTTGCTTCCGGCAGGTAACCGCTGGGTAAACGTTCATAg                                                                                                                                                        >  1:1834634/1‑126 (MQ=255)
                                                                                                                                                                                          ttGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGGATTTAGAATTGCCGGGCCCGCTTTccc                                                                                                                                                                                                                           >  1:1704670/1‑64 (MQ=255)
                                                                                                                                                                                           tGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCgg                                                                                    <  1:2363082/191‑1 (MQ=255)
                                                                                                                                                                                                  cTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGcc                                                                                                                                                    <  1:2199705/122‑1 (MQ=255)
                                                                                                                                                                                                           gTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCt                                                                                                                                 <  1:183307/131‑1 (MQ=255)
                                                                                                                                                                                                             tAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACtt                  <  1:38803/239‑1 (MQ=255)
                                                                                                                                                                                                             tAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACtt                  <  1:2298083/239‑1 (MQ=255)
                                                                                                                                                                                                                    ggTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCa            <  1:497533/238‑1 (MQ=255)
                                                                                                                                                                                                                     gTGAGGAAGGGATTAAGAAATGCCGGGCCCGCTTTCCAGATCTCATCCCGAAAAAAAAATTTTGCTTCCGGCAGGGAAACGCTGGGTAAACGTTCATAGCGCCCTTCCTCTTTTAGCGCCTGAATAGAAAg                                                                                                                      >  2:183307/1‑131 (MQ=255)
                                                                                                                                                                                                                           aaGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTgg    <  1:1613936/239‑1 (MQ=255)
                                                                                                                                                                                                                           aaGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTg     <  1:1461594/238‑1 (MQ=255)
                                                                                                                                                                                                                               tATTAAGAATTGCCGGGCCCGCTTGCCAGATCTCATCCAGAAAGTGCATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTGTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCTATCAAACTTTTGGCAATACCTGGCg  <  2:1435870/237‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
GGATTGTCGTTTTCATCCTGTTGACTGGTCAGCATGGAGCGGAAATTCGCTTTATCCTGCACTTTATCTAACCACAGACGAATCAGCATGCGGTCATATAACGCTTCCAGACTGCTGTCTGCTTCCGGCAGCTCGTTGGAGGCCGCCACCAGCAGACGCATCGGGATTTTTTCTACGTGTGCACCGTTGCGGAACTGGCGCTCGTTAATGGCGGTGAGCAAGGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCG  >  NC_000913/3930495‑3930954

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