Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	452,427	A→C	16.1%	G373G (GGT→GGG)	ampG ←	muropeptide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	452,427	0	A	C	16.1%	191.0 / 2.0	69	G373G (GGT→GGG)	ampG	muropeptide transporter
Reads supporting (aligned to +/- strand):  ref base A (42/14);  new base C (2/11);  total (44/25)
Fisher's exact test for biased strand distribution p-value = 1.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCGTCAATAGCAGTAACCACACGGCCAACA‑GGCTGA‑CGCCCGCCGCCAGTGTCCACATGGCAAAGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAAT  >  NC_000913/452202‑452657                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       ctacGTAAATAGAAGT‑ACCCCCCCGCCAACC‑AGCGGACCGCCCTCCGCCAGTGTCCACAGGGCAAGGGCATAACCTGCCGGATAGGGGGTGGGGGAGATAAAGTTGTCATTTCCTCGTGTATATACAAGCGTCTGGCGGCAACCCAGCAGCAAAAAAAGCCCTGGTAC‑GGCAGCGGCGACGGAAAATAGATAGAAGGTCGACCAGCGGT‑GTGCTAAAACAAACCA‑ACCCGCCACGGGGc                                                                                                                                                                                                                           <  2:1081830/236‑1 (MQ=255) catcaaaataaacataacccccacGG‑CAACAGGGCTGA‑CGCCCGCCGCCAGTGTCCACATGGCAAAGGCATACCCTGCCGAATAGGCGGTAGGGGAGATAAAGTGGTAATTTACTCGTGTATATTAAAGCGTCGGGCGGAAAACCACAAGAAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAAGGTCGCCCAGCCGT‑GTGCTTCAAAAAACCA‑CCCCGCCACGGGGc                                                                                                                                                                                                                           <  2:1083660/218‑1 (MQ=255)    cGTCAATAGCAGTAACCACACGGCCAACA‑GGCTGA‑CGCCCGCCGCCAGTGTCCACATGGCAAAGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCg                                                                                                                                                                                                                         >  1:1225992/1‑238 (MQ=255)              attAACCACACGGCCAACA‑GGCTGA‑CGCCCGCCGCCAGTGTCCACATGGCAAAGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTAAAGCGTCTGGCGGCAAACCAGCAGAAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑CCCCGCCACGGGGCCGACATAAACt                                                                                                                                                                                                                <  1:849522/235‑1 (MQ=255)                tAACCACACGGCCAACA‑GGCTGA‑CGCCCGCCGCCAGTGTCCACATGGCAAAGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGc                                                                                                                                                                                                             >  1:208591/1‑238 (MQ=255)                 acccccccGGCCAACA‑GGCGGA‑GCCCGGCCGCCAGTACCAACAGGGAAAAGGCATAACCTGCGGGATAGGCGGTAGGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCGGGCGGAAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCGGCGGCGACGGAGAATAGATAGAATGCCGACCGGCCGT‑GTCCTAAAACAAACCA‑ACCCGCCACGGGGCCGACATAAACt                                                                                                                                                                                                                <  2:111423/227‑1 (MQ=255)                   aacaccagggccacagggtggaggcccgccGCCAGTGTCCACATGGCAAAGGCATAACCTGCCGGTTATGCGGTACGGGAAAAAAAGTTGTCATTTACTGGTGTATATTCAAGCGTCGGGGGGAAACCCAGCAGAAAAATAAGCCCGGTTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGCCCAGCCGT‑GTGCTTCAACAACCCA‑ACCCGCCACGGGGCCGACATAAACTCg                                                                                                                                                                                                              <  2:665869/213‑1 (MQ=255)                                    ctgctggtttgccgccAGTGTCCACATGGCAAAGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAATCGCTg                                                                                                                                                                                          >  1:79723/11‑239 (MQ=255)                                    cTGA‑CGCCCGCCGCCAGTGTCCACATGGCAAAGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAAAAAGCGCTgc                                                                                                                                                                                         >  2:365208/1‑238 (MQ=255)                                      gA‑CGCCCGCCGCCAGTGTCCACATGGCAAAGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑CCCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTgagc                                                                                                                                                                                       >  2:269947/1‑238 (MQ=255)                                              gccgccAGTGTCCACATGGCAAAGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGAAAACCAGCAGCAAAAAAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAAGGCCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGc                                                                                                                                                                                <  2:1022242/239‑1 (MQ=255)                                                         ccACATGGAAAAGGCATAACCTGCCGGATATGCGGTACGGGAAATAAATTTTTCATTTCCCGGTTTATTTACAGGCGTCTGGCGGCAAACCAGCAGAAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGAAAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAg                                                                                                                                                                     <  2:208591/239‑1 (MQ=255)                                                                gCAAAGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGaaaa                                                                                                                                                              >  1:1135381/1‑239 (MQ=255)                                                                  aagatgtgtataagagcaagaggCGGTACGGGGAATAAAGTTGTCATTTACTCGTGTATATCAAAGCGTCGGGCGACAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAACCCA‑CCCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATg                                                                                                                                                            <  2:577973/217‑1 (MQ=255)                                                                  aaaggaccaaCCTGCCGGATATGCGTTCCGGGAGATAAATTTGTCATTTACTCGTGTATACTCAAGCTGCGGACGGAAAACCATCAGCAAAATAAGCCCGGGTAC‑GGCAGCGCCAACGGAAAAGGGATAGAATCGCGGCCAGCCGT‑GTGCTTCAACAACCCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATg                                                                                                                                                            <  2:682547/231‑1 (MQ=255)                                                                  aaaGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAACCCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATg                                                                                                                                                            <  1:160262/239‑1 (MQ=255)                                                                  aaaGGAAAAACCTGCCGGAAAGGCGGTACGGGAGAAAAAGTTGTAATTTCATCGGGTATATACAAGCGGCTGGGGACAAACCAGAAGAAAAATAAGCCCTGGACC‑GGAAGCGGCGACGGAGAATAGATAGAATGTCGCCCAGCCGT‑GTGCTAAAACAAACAA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATg                                                                                                                                                            <  2:1152858/239‑1 (MQ=255)                                                                    aGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACAA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATg                                                                                                                                                            >  1:447139/1‑237 (MQ=255)                                                                        aaaaCCTGCCGGTTTGGGGGTCCGGGAGATAAAGTTGTCTTTTCCTCGTTTATATCAAGGCGTCGGGCGGCAAACCAGCAGAAAAACAAGCCCGGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCCACCAGCCGT‑GTGCTTCAACAACCCA‑CCCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACtta                                                                                                                                                       <  2:1240236/236‑1 (MQ=255)                                                                        aTAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGCCCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACtta                                                                                                                                                       <  1:531124/238‑1 (MQ=255)                                                                         tAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACttatt                                                                                                                                                     >  1:966561/1‑239 (MQ=255)                                                                         tAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACttatt                                                                                                                                                     >  1:797524/1‑239 (MQ=255)                                                                            ccTGCCGGAAATGGGGTACGGAAGAAAAATTTGTCATTTACCGGTGTATATTCAAGCGTCTGGGGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGAGGGAGAATAGATAGAATGTCGACAAGCCGT‑GTGCTTCAACAAACCA‑CCCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACa                                                                                                                                                  <  2:21303/239‑1 (MQ=255)                                                                              tccagtatagGCGGTACGGGAGATAAATTTGTTATTTACTGGTGTATATACAAGCGCCTGGCGGAAAACCGGAGACAAAATAGGCCCGGGTAC‑GGAAGCGGCGACGGAAAAGAAAAAAAATGCCAACCAGCCGT‑GTGCTTAAACAAACCA‑CCCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATa                                                                                                                                                <  2:875215/233‑1 (MQ=255)                                                                              ggcggaatagggggTAGGGAAAATAAATGTGTCATTTACTCGGGTATATCCAAGCGTCTGGCGGAAAACCAGCAGCAAAATGAGCCCTGGTAC‑GGCAGCGGCGACGGAAAATAGATAGAATGGCGCCCGGCCGT‑GTGCTTCAACAAACCA‑CCCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATa                                                                                                                                                <  2:79723/227‑1 (MQ=255)                                                                                ccGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATg                                                                                                                                                            >  2:1045431/1‑225 (MQ=255)                                                                                ccGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATg                                                                                                                                                            <  1:1045431/225‑1 (MQ=255)                                                                                ccGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGc                                                                                                                                              >  1:501268/1‑239 (MQ=255)                                                                                     tatGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCAt                                                                                                                                         >  1:373198/1‑239 (MQ=255)                                                                                      atGCGGTACGGGAGATAAAGTTGTCTTTTACTCGTGTATATTCAAGCGTCTAGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGATAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCAtt                                                                                                                                        >  2:69699/1‑239 (MQ=255)                                                                                            taggggAAATAAAGTTGTAATTTCTCCGTGTATATTCAATCGTCTGGCGGCAAACCGGCAGAAAAAAAAGCCCTGGTAG‑GGCAGGGCCGACGGAAAAAAAATAAAAGCCCGACCACCCGT‑TTGTTTCAACAAACCA‑CCCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAg                                                                                                                                   <  2:844394/235‑1 (MQ=255)                                                                                              cGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCAGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAgcgc                                                                                                                                >  1:844394/1‑239 (MQ=255)                                                                                                 gagaTAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACa                                                                                                                                                  >  1:577973/1‑218 (MQ=255)                                                                                                 gagaTAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACATATTACATAGCGTCATTAACAGCGCGAc                                                                                                                             >  2:64493/1‑239 (MQ=255)                                                                                                   gaTAAAGTTCTCTTTTACCCTTTTATATTCAAGCGTCTGGCGCCAAACCAGCAGCAAAATAACCCCTGGTTC‑GGCAGCGGCGTCGGAGAATAAAAAGAATGTCGTCCAGCCGT‑GTTCTTCAAAAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGAc                                                                                                                             <  2:960855/237‑1 (MQ=255)                                                                                                     tAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAg                                                                                                                         >  1:877221/1‑239 (MQ=255)                                                                                                      aaaGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAg                                                                                                                         >  1:111950/1‑238 (MQ=255)                                                                                                        aatatccATTTGCTCGTTTATATCCAAGCGTCAGGCGGAAAACAGGCAGAAAAAAAAGCCC‑GGAACGGGCAGCGGCGACGAAGAATAGATAGAATGCGCACAGGCCGT‑TTGCTTCAACAAACCA‑CCCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACtta                                                                                                                                                       <  2:251477/200‑1 (MQ=255)                                                                                                         tttGTCATTAACTCGTTTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAAAAAGCCCGGGTCC‑GGCAGGGGCGACGGAGAATAGATAGAATGTCGCCCAGCCGT‑GTGCTTAAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTg                                                                                                                     <  1:112915/238‑1 (MQ=255)                                                                                                             tCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACCTAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGAtgg                                                                                                                 >  2:1102838/1‑238 (MQ=255)                                                                                                                  tACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCg                                                                                                                                             >  1:251477/1‑206 (MQ=255)                                                                                                                         gtATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGcc                                                                                                     >  1:964349/1‑239 (MQ=255)                                                                                                                                 aaGCGGCTGGTGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTCCAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGCAAATGACTTATTACATAGCGACAATCACAGCGCGAGAAAGGCACACGTTCACATCCCGCCACAGAGGt                                                                                             >  2:680805/1‑239 (MQ=255)                                                                                                                                      tCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCTGACACAGAGGTTTTCg                                                                                        >  2:622752/1‑239 (MQ=255)                                                                                                                                      tCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCg                                                                                        >  2:449408/1‑239 (MQ=255)                                                                                                                                        tGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGaa                                                                                      >  2:814564/1‑239 (MQ=255)                                                                                                                                                           aaaaaaagcccgggtaacgggcgggcgACGGAAAAAAAAATAAATGTCCGCCCGGCGTGGTGCTAAAAAAACCCA‑CCCCGCCACGGGGCCGACATAAACCCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCa                                                                     <  2:966561/216‑1 (MQ=255)                                                                                                                                                                      ggTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGAt                                                         >  1:200276/1‑238 (MQ=255)                                                                                                                                                                      ggTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GCGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGAAAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGAt                                                         >  1:25921/1‑238 (MQ=255)                                                                                                                                                                       gTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGATTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGAt                                                         >  2:23062/1‑237 (MQ=255)                                                                                                                                                                                gcggcgACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAg                                                                                                                                                                     >  1:1019976/1‑121 (MQ=255)                                                                                                                                                                                gcggcgACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAg                                                                                                                                                                     <  2:1019976/121‑1 (MQ=255)                                                                                                                                                                                 cggcgACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTa                                              >  1:531577/1‑239 (MQ=255)                                                                                                                                                                                        ggcgaATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑CCCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATAGAGTAGCGGCAAAGGACTTACTACATAGCGGCATTAACAGCGCGACATGGGCTGAGGTACCCATCCCGCCACAGAGGTTTCCGAAAAAGACGGCTGCGCCCACGCTGGAGAGTTGCCTATCATTACTCGAca                                       >  2:354269/4‑239 (MQ=255)                                                                                                                                                                                                   agaATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTa                             >  1:547936/1‑238 (MQ=255)                                                                                                                                                                                                    gaATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTaa                            >  1:947433/1‑238 (MQ=255)                                                                                                                                                                                                      aTGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCg                         >  2:334141/1‑239 (MQ=255)                                                                                                                                                                                                      aTGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCg                         >  1:460817/1‑239 (MQ=255)                                                                                                                                                                                                         tCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAGGGCTGATGTGCCCATCCCGCCACAGAGGTATTCGAAAAAGACGCCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCg                      >  1:239176/1‑239 (MQ=255)                                                                                                                                                                                                          cGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGt                     >  1:539551/1‑239 (MQ=255)                                                                                                                                                                                                           gACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGtt                    >  1:562754/1‑239 (MQ=255)                                                                                                                                                                                                            aCCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTa                   >  1:880803/1‑239 (MQ=255)                                                                                                                                                                                                            aCCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTa                   >  2:453731/1‑239 (MQ=255)                                                                                                                                                                                                               aGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGAc                >  1:364714/1‑239 (MQ=255)                                                                                                                                                                                                                 ccgggttcctaaaaaaaaCCA‑CCCCGCCACGGGGCCGAAATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGCCTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACg               <  2:539551/224‑1 (MQ=255)                                                                                                                                                                                                                       tttctcaaaaaaacacaccCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGAAGGGCAAATTGAGTAGCGGAAAATGACTTATTACAGAGCGTAATTAACAGCGCGAAAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCAcc           <  1:622752/222‑1 (MQ=255)                                                                                                                                                                                                                       tttctaaaaaaaCCACACCCGCCACGGGCCGGACATAAACCGGCCATACAGCAGAAGGGGCTGAGAGAAGGGCAAATTGAGTAGCGGAAAAGGACTTATTACATCGCGTCATTAACAGCGCGAAAAAGGTTGAGGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCAcc           <  2:200276/229‑1 (MQ=255)                                                                                                                                                                                                                       tgCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCtt         >  1:1168779/1‑239 (MQ=255)                                                                                                                                                                                                                            caacaaACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTCTTCGAAAAAGACGGCAGCGCCCATGCTGTAGAGATGCTTATCTGTAATCGTCAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTaaa    >  1:864468/1‑239 (MQ=255)                                                                                                                                                                                                                               caaACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAAt  >  2:436945/1‑238 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                       CATCGTCAATAGCAGTAACCACACGGCCAACA‑GGCTGA‑CGCCCGCCGCCAGTGTCCACATGGCAAAGGCATAACCTGCCGGATATGCGGTACGGGAGATAAAGTTGTCATTTACTCGTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTAC‑GGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGT‑GTGCTTCAACAAACCA‑ACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAAT  >  NC_000913/452202‑452657

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: ‑