Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	452,547	G→C	14.9%	G333G (GGC→GGG)	ampG ←	muropeptide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	452,547	0	G	C	14.9%	170.0 / 2.8	70	G333G (GGC→GGG)	ampG	muropeptide transporter
Reads supporting (aligned to +/- strand):  ref base G (47/11);  new base C (2/10);  total (49/21)
Fisher's exact test for biased strand distribution p-value = 4.06e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTGTTAACCACGCCTACTTCACCCGCATCA  >  NC_000913/452318‑452781                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             gtgtATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTACGGCAGCGGCGACGGAGAATAGATAGAATATCGACCAGCCGGGGGCTTCAACAAACCAACCCGCCACGGGGCCGACCTAAACTCGCCCTACAGCAGAAAACGCTGAGAACAGGGCCAATTGAGTAGCGGCAAATGACTAATTACCTAACGCCATTAACAACCCGACAAAAACCAATGTG‑CCCA‑Tccc                                                                                                                                                                                                                                    >  1:161107/1‑238 (MQ=255)        ttCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACaga                                                                                                                                                                                                                            >  2:1161049/1‑239 (MQ=255)         tCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAACACAACCCGCCACGGGGCCGACATAAACCCCCCCTCCAGCCGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTAATACAAAGCGTCATTAACAACGCAGCAAAGGCTCAAGTG‑CCCC‑ACCCACCACaaag                                                                                                                                                                                                                           >  1:283459/1‑239 (MQ=255)             gCGTCTGGCGGCCAACCAGCAGCAAAATAAGCCCTGGTACGGCAGCGGCGACGGAGACTAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCAAGGGGCGCGCATAAACTCGCCCTACAGCAGAAAGCCCTGGGAGCAGGGCAACTTAAGTAGAGGAAAATGACTACTTCAATAGCATCATTAACACAGCGGCCAAGGCTGAAGTG‑CCCC‑TCCCCCAACAAAGGttt                                                                                                                                                                                                                       >  1:225497/1‑239 (MQ=255)             gCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGGG‑CCCA‑TCCCGCCACAGAGGttt                                                                                                                                                                                                                       >  2:640494/1‑239 (MQ=255)                  tGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGaa                                                                                                                                                                                                                  >  2:599019/1‑239 (MQ=255)                                     aaaTAAGCCCTGGTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGAGTCATTACCAGCGCGACAACGGCAGATGGG‑CCCA‑TCCCCCCCCAGAGGGTTTCGAAAAAAACGGAGGAGCCCac                                                                                                                                                                                                >  1:15245/1‑237 (MQ=255)                                           gCCCTGGTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGCAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGGAAATGACTTATTACATAGCGCCCTTAACACCGCCACACAAGCTGCAGGG‑CCCA‑GCCCCCCACACAGGTTTTCCCAAAAGACGCAGGCGCCCAtccgttag                                                                                                                                                                                         >  1:168733/1‑232 (MQ=255)                                                ggTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGCCTTATTACATAGCGTCATTAACCGCGCGACAAAAGCTGATATG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAAACCGCTGCGCCCATTCTGTAaagaag                                                                                                                                                                                    >  1:628009/1‑237 (MQ=255)                                                ggTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACCTAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCCCCACAGAGATTTTCGAAAAAAACGGCCGCGCCCATGCTGTAgagaa                                                                                                                                                                                     >  1:358574/1‑237 (MQ=255)                                                ggTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGAGGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGAt                                                                                                                                                                                     >  2:1260570/1‑238 (MQ=255)                                                       cAGCCGCGACCGAGCACAGATAGACTATCGACCAGCCCTGTGCTTCAACAAACCAACCCGCCACGGGGCCGCCATAAACCCGCCACACACCAGAACGCCCTGAGACCGGGGCCAATTGAGTATCGGAAAAGTACATATTACAACACTTCATAAACAGCCGGACAAATGCCGATGTC‑CACA‑TACCTGCACAGCGGTTCTAGAAAAAAACGGCCGCGCCCagcctctaggaaatcacaacc                                                                                                                                                                             >  1:401803/1‑219 (MQ=255)                                                                             gaATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGACTGATGTG‑CCCA‑TCCCGCCACAGAGGCTTTCGAAAAAGACGGCGACGCCCCTGCTCTAGAGATGCTTATCATTAATCGACAGCAGCCAGTaa                                                                                                                                                        >  2:322524/1‑238 (MQ=255)                                                                               aTGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCg                                                                                                                                                     >  2:502153/1‑239 (MQ=255)                                                                               aTGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAACGGAAAATGACTTATTACATAGCGTCATTACCAGCGCGACAAAGGCTGATGAG‑CCCA‑TCCCGCCACAGAGGGTTTCGTACAAGACGGCTGCGACCACACTGCAGCGCTGCATAGCAGAAATTGAAAGAAGCCAGGAACCg                                                                                                                                                     >  1:733326/1‑239 (MQ=255)                                                                                  tCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTTTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCg                                                                                                                                                  >  2:995495/1‑239 (MQ=255)                                                                                          ccggcccctacaaaaaCCCACCCCCCCCGGGGGCCAAAAAAAACCCCCCCAAAAGCAGAAAGCGCGGAGAGGAGGCCAAATTGAGTAGGGAAAAATGACTTATTACAATGCGTCAATAACAGCGCGAAAAAGGCTATTGTC‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGc                                                                                                                                          <  1:1260570/229‑1 (MQ=255)                                                                                          ccGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACg                                                                                                                                           >  2:24161/1‑238 (MQ=255)                                                                                           ggtggctaaaaaaaacccacccCGCCAGGGGGCCGACATAAACTCGCCCTACAGAAGAAAGCGCTGAGAGCAGGGAAAATGGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGTTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCa                                                                                                                                         <  2:358574/223‑1 (MQ=255)                                                                                            gtgtGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGAAAGCACCCAGCAACCGGGGTTAGACgaa                                                                                                                                         >  1:864383/1‑236 (MQ=255)                                                                                                    caacaaACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAACGACGGTCGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCCGCCAGGAAACCGCGTgagaagacccgtggaaa                                                                                                                                >  1:1209282/1‑226 (MQ=255)                                                                                                    caacaaAACAACCCGCCCCGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGCATTATTACATAGCGTCATTACCAGCGCGACAAAGGCTAATGGC‑CCCC‑ACCCGCCACAGAGGCTTTCGAAAAAGAAGACCGCGCCCATACTGTAGCGATGCTCATCAGAAATCCAAAGCAACCAGGAACCAGCGTTACACGAACCtaaaaat                                                                                                                                >  1:714415/1‑233 (MQ=255)                                                                                                          cccccccccccAGGGGCCCGAAAAAACCCCCCCCTACAGGAAAAAGCGCTAAGAGCAGGAAAAGTGGAGGAGGGGAAAATGACTAATACCATGGCGTAATAAACGGCGAGCAAAGGGGTGATTTC‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCg                                                                                                                          <  1:640494/234‑1 (MQ=255)                                                                                                          accaaccCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAACCGGCTGCGCCCATGCTGCAGAGATGCGTATCAGTAATCGACCACAGCCCGTAACCGGCGTTAAACGCCCCTTGTAAAATAcca                                                                                                                          >  1:1022241/1‑238 (MQ=255)                                                                                                               cccGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGCGGTTTTCGAAAAAGACGGCTGAGCCCATGCTGTAGGGATGCTTACCAGTAATCGAAAGCAGCCAGGAACCGGCGTTAGACGCACCTTGTAAAATGCCGaaaac                                                                                                                     >  2:517486/1‑238 (MQ=255)                                                                                                                          gccaaaaTAAACCCGCCCTACAGCAGAAAGCGCTAAGAGCAGGCAAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCAACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGcc                                                                                                           <  2:367121/232‑1 (MQ=255)                                                                                                                             caaaaaaaCCCCCCCTACAGCAAAAAGCGCTGAGAGGAGGGCAAATTGAAGAGGGGAAAATGATTTATTACAGAGGTTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGa                                                                                                       <  2:994821/234‑1 (MQ=255)                                                                                                                              aCATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGCCACCAGCCAGTAACCGGCGTTAGACGAACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGaa                                                                                                      >  2:573912/1‑239 (MQ=255)                                                                                                                              aCATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATAAGCAGTGCCCGGaa                                                                                                      >  2:697543/1‑239 (MQ=255)                                                                                                                                    aCTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTg                                                                                                 >  2:1057646/1‑238 (MQ=255)                                                                                                                                         cccTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGgc                                                                                           >  2:1099658/1‑239 (MQ=255)                                                                                                                                          ccTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCCCCTTGTAAAATGCCGAAAATCAGCAGGGCCCGGAACAGTGACAGgc                                                                                           >  2:1046412/1‑238 (MQ=255)                                                                                                                                                    aaaGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATATCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGCTGCTTATCAGTAATCTACAGCAGCCAGTAACCTGCGTTAGACGCACCTAGAAAAATGCCGAAAATCAGCAGAGCCCGGAAAAGGGACAGGCGCCACATCaaa                                                                                >  1:785416/1‑239 (MQ=255)                                                                                                                                                    aaaGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGTTGCTTATCAGGAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCATGGCCCGGAACAGTGACAGGCGCCGCATCaaa                                                                                >  2:802823/1‑239 (MQ=255)                                                                                                                                                    aaaGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCaaa                                                                                >  2:1025913/1‑239 (MQ=255)                                                                                                                                                     aagccgtgaaggCAGGGCAAATTGGGTAGCGGAAAATGACTTATTACATGGCGTAATTAACAGCGCGAAAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCaaaa                                                                               <  2:1209282/228‑1 (MQ=255)                                                                                                                                                       taaaagagacgggggAAAATTGAGTGGGGAAAAATCTCTTATTAAAGGGCTTAATTAACAGCGCAAAAAGGGCTGAGCGG‑CCAA‑CCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATc                                                                             <  1:437045/227‑1 (MQ=255)                                                                                                                                                         gcTGAAAGAAGGGCAAATTGAGTAGCGGAAAAGGACTTATTACATAGCGTCATTAACAGCGCGAAAAAGGCTGAGGTG‑CCAA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATc                                                                             <  2:628009/237‑1 (MQ=255)                                                                                                                                                                  aGGGAAAATGGAGTAGGGGAAAATGACTTATAACAAGGCGTCATAAACCGGGCAAAAAAGGCTGATGTC‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACa                                                                   <  2:1194014/238‑1 (MQ=255)                                                                                                                                                                       aaaTTGAGTAGCGAAAAATTACTTATTTAATAGCGTCATTAACAGCGCGAAAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATgcg                                                              <  2:1167119/238‑1 (MQ=255)                                                                                                                                                                         aTTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCca                                                           >  2:215231/1‑239 (MQ=255)                                                                                                                                                                              gTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACa                                                                   >  2:437045/1‑226 (MQ=255)                                                                                                                                                                                   ggAAAAGGACTTATTACATAGCTTCATTAACAGCGCGACAAAGGCTGATGTC‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGt                                                   <  2:15245/237‑1 (MQ=255)                                                                                                                                                                                   gaaaatttctttttcaCATAGCGTCATTAACAGCGAGAAAAAGGGTTGGGTC‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGt                                                   <  1:215231/222‑1 (MQ=255)                                                                                                                                                                                    gAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCa                                                                          >  2:58682/1‑213 (MQ=255)                                                                                                                                                                                            cTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCTTACAATGCGCCACCAATGGTCGCTAGTAAGc                                        >  2:1071741/1‑239 (MQ=255)                                                                                                                                                                                              ttttccATGGCGTTAATAACCGCCCGAAAAGGGCTG‑TGTGCCCCC‑CCCCGACCAAGAGGTTTTCGAAAAAGACGGCTGCGCCCCTGCTGTAGGGATGCTTATCAGTAACCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCa                                      <  1:1071741/234‑1 (MQ=255)                                                                                                                                                                                               attaCATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGGCGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCTCACAAAGCGCCAACAAAggcccccagaaaagaaa                                     >  1:494100/1‑224 (MQ=255)                                                                                                                                                                                                 taCATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGCACAGTGACAGGCGCCACATCCAAATCCACCCCTACAATCCGCCAACAAAGGTCGCCAGTAAGCCaacc                                   >  1:862145/1‑237 (MQ=255)                                                                                                                                                                                                     tAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGttt                               >  2:1306652/1‑239 (MQ=255)                                                                                                                                                                                                      aGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGtttt                              >  2:813565/1‑239 (MQ=255)                                                                                                                                                                                                        cgttattaacgggccaaaaaaGGCTGATGTC‑CCCA‑CCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGtt                                <  2:168733/217‑1 (MQ=255)                                                                                                                                                                                                        cGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATACCACCGTACAACGCGCCAACACAAGCCGCCAGTAAACCAAGCGTTTTGt                            >  1:600021/1‑239 (MQ=255)                                                                                                                                                                                                          tCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTg                             >  2:229707/1‑236 (MQ=255)                                                                                                                                                                                                            ataaaCAGGGCGAAAAAGGCTGATGTC‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTGTTAAc                        <  2:1022241/236‑1 (MQ=255)                                                                                                                                                                                                            aTTAACAGCGCGACAAGGGCTGATGGC‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTGTTAAc                        <  2:347777/239‑1 (MQ=255)                                                                                                                                                                                                            aTTAACAGCGCGACAAAGGTTGAGGGC‑CCCA‑CCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTGTTAAc                        <  2:972785/239‑1 (MQ=255)                                                                                                                                                                                                             ttAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACg                                                                                                                                           >  1:584680/1‑123 (MQ=255)                                                                                                                                                                                                             ttAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACg                                                                                                                                           <  2:584680/123‑1 (MQ=255)                                                                                                                                                                                                             ttAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAACCAGCAGGGCCCGCAACGGTGACAGGCGCTGCATTCAAATCCACCCGTACACTGAGCAAAAAATGGAGGCTAAGAAACAAAACGATTTTTTaaac                       >  1:619492/1‑237 (MQ=255)                                                                                                                                                                                                                aCAGCGGGAAAAGGGCTGATTTG‑CCCA‑CCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTGTTAACCACg                    <  2:1094629/239‑1 (MQ=255)                                                                                                                                                                                                                  aGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACg                                                                                                                                           <  2:973072/118‑1 (MQ=255)                                                                                                                                                                                                                  aGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACg                                                                                                                                           >  1:973072/1‑118 (MQ=255)                                                                                                                                                                                                                    cgcgACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAACTCAGCAGTGCCCGGAACAGTGCAAGGCGCTGCATCAAAATCCCACCGAACAATGCGCCACAAATGGCCGCCAGTAAACCAAGGGTTTTGTTaacaaccgac                 >  1:560185/1‑234 (MQ=255)                                                                                                                                                                                                                    cgcgACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGGACAGTGACAGGCGCTGCATCCAAATCCCACCGTAAAAACGGCCAACAAGGGTAGCCAGGAAGCCAAACGATTTTTTAAACACGcca                 >  1:1094629/1‑237 (MQ=255)                                                                                                                                                                                                                    cgcgACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTGTTAACCACGCCt                 >  2:1217083/1‑238 (MQ=255)                                                                                                                                                                                                                    cgcgACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTGTTAACCACGCCt                 >  2:242468/1‑238 (MQ=255)                                                                                                                                                                                                                    cgcgACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGAGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTCGAAAGCCAAGCGTTTTGTTAACCACGcca                 >  1:247145/1‑237 (MQ=255)                                                                                                                                                                                                                       gaaaaggggtGATGTC‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTGTTAACCACGCCTACt              <  2:864383/229‑1 (MQ=255)                                                                                                                                                                                                                            aggggttgttg‑CCCACTCCCGCCAAAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCCCCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTGTTAACCACGCCTACTTCAcc         <  2:1101291/230‑1 (MQ=255)                                                                                                                                                                                                                                  aTGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGGGCCCGGACCAGTGACAGGCGCTGCATCAAAATCCCACCGGACAATGCGCCAACAATGGTCGCCAGTAACCCAAGCGCTTTGTTAAcacccccaaacccaccgcagcaa  >  1:1019620/1‑217 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             GTGTATATTCAAGCGTCTGGCGGCAAACCAGCAGCAAAATAAGCCCTGGTACGGCAGCGGCGACGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTG‑CCCA‑TCCCGCCACAGAGGTTTTCGAAAAAGACGGCTGCGCCCATGCTGTAGAGATGCTTATCAGTAATCGACAGCAGCCAGTAACCGGCGTTAGACGCACCTTGTAAAATGCCGAAAATCAGCAGTGCCCGGAACAGTGACAGGCGCTGCATCAAAATCCCACCGTACAATGCGCCAACAATGGTCGCTAGTAAGCCAAGCGTTTTGTTAACCACGCCTACTTCACCCGCATCA  >  NC_000913/452318‑452781

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