Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,630,921	G→A	P63P (CCG→CCA)	slt →	soluble lytic murein transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,630,921	0	G	A	92.8%	177.5 / ‑0.5	71	P63P (CCG→CCA)	slt	soluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (35/29);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑ATGTAATTACTACTATGCGAAATGGAACACCGGGC  >  NC_000913/4630677‑4631162                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        aaTCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCATCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGTAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTAt                                                                                                                                                                                                                                                   <  1:405626/251‑1 (MQ=255) aaTCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTAt                                                                                                                                                                                                                                                   <  1:270830/251‑1 (MQ=255)            cgGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTgg                                                                                                                                                                                                                                         >  1:286030/1‑250 (MQ=255)            cgGCATGACGCTGCATGGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTgg                                                                                                                                                                                                                                         >  2:489045/1‑250 (MQ=255)                aTGACGATTCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTCCCGGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGAAGCGTGGCGCGAGCCGACTCACGTGATGAGCAGAGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTTGAATAc                                                                                                                                                                                                                                    <  2:295135/251‑1 (MQ=255)                aTGACACTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATAc                                                                                                                                                                                                                                    >  1:102375/1‑251 (MQ=255)                     gCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCGGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCgtc                                                                                                                                                                                                                                >  2:443889/1‑248 (MQ=255)                          aTTGATTTTTTTCCACTTTGAGTATGCGCTGGTGGAAAAAGCTAAAAAAGTTACCTGGCGGCTGGTGGCTGCCGGTGTCTGTCTGTTTACGTGCAGCATCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATAAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTAGCCCTACCTTGAAGACCGTCAGATCa                                                                                                                                                                                                                          <  2:477726/251‑1 (MQ=255)                          aTTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTATGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAgctc                                                                                                                                                                                                                           >  1:627631/1‑247 (MQ=255)                          aTTGATGTATTTACACTTAGAGGAGGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTTCCTGGAATACCGTCAGATCa                                                                                                                                                                                                                          <  1:464466/251‑1 (MQ=255)                             gATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACACACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCg                                                                                                                                                                                                                       >  2:71199/1‑251 (MQ=255)                                           tAGAGGATGCGCTTGGGGAAAAAGCCAAACAAGTTACCTGGCGGGTGTTGGCTGTCGGTGTCTGTCTGCTGACGGTCAGGAGCGTGGCGCGAGCCGACTCCCTGGATGAGCAGCGTAGCCGTTACGCGAAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATCCCGTCAGATCACCGATGACCtgatga                                                                                                                                                                                                           <  1:126402/249‑1 (MQ=255)                                            agagGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATc                                                                                                                                                                                                        >  1:559521/1‑251 (MQ=255)                                                 aTGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACcg                                                                                                                                                                                                   <  1:496349/251‑1 (MQ=255)                                                        tgtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGACTGCTGACGGTCAGCAGCGTGGCGCGAGCCGCCTCACTGGATGAGCAGCGTAGCAGTTACGCGCAAATCATGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAAAGATGCCTGGACTGA‑AAGATTATCCACATTATCCCTACCTGGCATACCGTCAGATCACAGATGACCTGGTGAAACAAAcggcag                                                                                                                                                                                               >  2:366454/1‑246 (MQ=255)                                                         gtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGACGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGGTCACCGATGACCTGATGAATCAACCGGCGGTg                                                                                                                                                                                             >  2:82426/1‑249 (MQ=255)                                                         ggggaaaaccccaacaaaGTTACCTGGCGGGTGTGGGCGTCCGGTGTCGGTCTTCTGACGGTCAGCAGCTTGGAGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCACATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AATATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑c                                                                                                                                                                                          <  2:589543/235‑1 (MQ=255)                                                           ggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGTCTGCCGGTGCCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑ctg                                                                                                                                                                                        <  1:352323/251‑3 (MQ=255)                                                           ggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑ctg                                                                                                                                                                                        <  1:461966/251‑3 (MQ=255)                                                                        cAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGA‑CTGT‑CCCTAAATTTGt                                                                                                                                                                          >  2:20449/1‑251 (MQ=255)                                                                         aaGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAACTCAAGCAGGCCTGGGATAATCGACAAATGTATGTGGTCGAACAAATGCTACCTGGACTGC‑AAGATTAACCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGTTGAATCAACCGGCGGTAAATGGT‑‑ACTAACTTTGtt                                                                                                                                                                         >  2:353660/1‑251 (MQ=255)                                                                          aGTTACCGGGCGGCTGTTGGCTGCCGGTGTCTGTTGGCTGACGGTCAGGAGCGTGTCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCTGCGGTTA‑CTGT‑CACTAAGTTTGTTc                                                                                                                                                                        <  2:598542/251‑1 (MQ=255)                                                                              aCCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGCCGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTg                                                                                                                                                                                                                                          >  2:180581/1‑183 (MQ=255)                                                                              aCCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTg                                                                                                                                                                                                                                          <  1:180581/183‑1 (MQ=255)                                                                                 tGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CCCTAACTTTGTTCGCGCTa                                                                                                                                                                  >  1:353216/1‑250 (MQ=255)                                                                                   gcggcTGTTGGCTGCCGGTTTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgatga                                                                                                                                                                                                           <  2:364474/209‑1 (MQ=255)                                                                                   gcggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgatga                                                                                                                                                                                                           >  1:364474/1‑209 (MQ=255)                                                                                      gcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCa                                                                                                                                                            <  1:20449/251‑1 (MQ=255)                                                                                        tGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGA‑ACGT‑CACTAACTTTGTTGGCGCTTACCCCAcc                                                                                                                                                          >  2:373474/1‑250 (MQ=255)                                                                                           tGGCTGCGGGTGTCTGTATGCTGACGGTCCGCAGCGTGGCGCGAGACGCCTACCTGGATGAGCAGCGTAGCCGTTAAGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACACATGATGCCTGGACTGACAGG‑TTATCAAATTTATCCCTACCTGGAATACCGTCAGATCACCGATGCCCTGCTGAATACACCGGCGGTGT‑CTGTGCA‑TAACTTTGTTCGCTCTAAACCCACGcat                                                                                                                                                       >  2:335787/1‑249 (MQ=255)                                                                                                gccgggggcggcctgctgACGGTCAGCAGCGTGGCGCGAGCCGACTCACTTGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGTGGA‑CTGT‑CTCTAACTTTGCTCGCGCTAACCCCACGCTtccg                                                                                                                                                    <  2:248744/237‑2 (MQ=255)                                                                                                      gtctgtctGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACTCTTCCTCCCGCTCg                                                                                                                                            <  1:560023/251‑1 (MQ=255)                                                                                                       tctgtctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGc                                                                                                                                           >  1:618920/1‑251 (MQ=255)                                                                                                               ctgACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGAGGTGA‑CTGT‑CACTAACTTTGTTATCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGcc                                                                                                                                   >  2:627777/1‑250 (MQ=255)                                                                                                                       cagcagCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGCAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGCTGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGtt                                                                                                                           >  1:279981/1‑251 (MQ=255)                                                                                                                           agcGTGGCGCGAACCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACATGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTACCGCTCGCACGCTGCCATCTCGTTTcat                                                                                                                       >  1:318062/1‑249 (MQ=255)                                                                                                                                  cgcgAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGCCTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGG‑CACTAACTTTGTTCACGCTAACCCCACGATTCCTCCCGCTCGAAAGCTGCAATCTCGTTTCGTCAATGa                                                                                                                 >  1:1980/1‑250 (MQ=255)                                                                                                                                    cgAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTa                                                                                                                                                                                                                                                    >  1:106730/1‑119 (MQ=255)                                                                                                                                    cgAGCCGACTCACGGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGTGATAATCGACAAATGGATGTGGTCGAACAAAGGATGCCTGGACTGA‑AAGATTATCCACTTTa                                                                                                                                                                                                                                                    <  2:106730/119‑1 (MQ=255)                                                                                                                                      aGCCGACTCACTGGATGAGCCGCGTAGCCGTTACGCGCAAATCAAGCAGTCCTTGGATAAACGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTTA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCGCGTTTCGTCAATGAACTgg                                                                                                            <  2:87369/251‑1 (MQ=255)                                                                                                                                         cGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGTACTGA‑AAGATTAGCCACTTTa                                                                                                                                                                                                                                                    <  2:671252/114‑1 (MQ=255)                                                                                                                                         cGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTa                                                                                                                                                                                                                                                    >  1:671252/1‑114 (MQ=255)                                                                                                                                         cGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTAACCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGgcgc                                                                                                         <  1:353660/251‑1 (MQ=255)                                                                                                                                               aCTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTg                                                                                                   <  1:510282/251‑1 (MQ=255)                                                                                                                                                      gAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCATGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTg                                                                                            >  1:223674/1‑251 (MQ=255)                                                                                                                                                         cagcGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCg                                                                                         <  1:373474/251‑1 (MQ=255)                                                                                                                                                                gccGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTg                                                                                  >  1:412412/3‑251 (MQ=255)                                                                                                                                                                 ccGTTACGCTCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTTGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCTCCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGt                                                                                 <  2:321099/250‑1 (MQ=255)                                                                                                                                                                    ttACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGGGGCtttttct                                                                              >  2:533931/1‑246 (MQ=255)                                                                                                                                                                             aaTCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGc‑c                                                                   <  1:660107/251‑1 (MQ=255)                                                                                                                                                                               ccAAGCAGGCCTGGGATAATCGACAAATGGATGTGGGCGAACAAATGATGCCTGGACTGA‑AAGATTATCCCCTTTATCCCTACCTGGAATACCGTCAGATCACCGATTACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGc‑cag                                                                 <  2:538987/250‑3 (MQ=255)                                                                                                                                                                                  aGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑aa                                                             <  1:533939/251‑1 (MQ=255)                                                                                                                                                                                        ccTGGGATAAGCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATCCGCTTTATCCCGACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCGCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGccc                                                       <  2:193017/251‑1 (MQ=255)                                                                                                                                                                                                  acgtcAAATGGATGTGGTCGAACAAATGATGCCTGGACTGG‑AAGATTATCCGCTTTATCCCTACCTGGAATACCGTCAGACCGCCGACGACCTGATGATTCAGCCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCGCGCTTCCGCCCGCTCGCACGCTGCAATCTCGGTTCGTGAATGCACCGGCGTTTCGTGAGCAGTGGCGTGTATTGTTAG‑CGTTTAGC‑CCGG‑AAAGGCCGGGatctgccg                                             <  2:318266/247‑8 (MQ=255)                                                                                                                                                                                                   cGACAAATGGATGTGGACGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACCAACTTTTGTCCCGCTAACCCCAAGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGACAATGAACTGGGGCGGCGTGAGGGAGGGCGTGGCTTTTTAGGCCTTTTAC‑CCGG‑aaaaaccgggacaccatg                                             >  2:178438/1‑237 (MQ=255)                                                                                                                                                                                                     aCAAAGGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAAGCCCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAg                                          <  2:353216/251‑1 (MQ=255)                                                                                                                                                                                                      cAAATTGATGTGGTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAgc                                         <  2:286030/251‑1 (MQ=255)                                                                                                                                                                                                                 gggcGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑ATGTAAtt                             <  1:294110/248‑1 (MQ=255)                                                                                                                                                                                                                  gTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGACCTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AACAGACCGGAACTACCGAAGCGCA‑CAGTAAtt                             >  1:51851/1‑250 (MQ=255)                                                                                                                                                                                                                    cGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTTTTCT‑CCTTTATCACCGG‑AAAGACCCGGA‑CAACCCAAGCGCC‑ATGTAAAtact                          >  2:287449/1‑251 (MQ=255)                                                                                                                                                                                                                    cGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTTTTAG‑CCTTTTGC‑CCGG‑AAAAACCCGGGACCACCAAAGCGCattgaacttact                           >  2:447467/1‑240 (MQ=255)                                                                                                                                                                                                                           atgatgCCTGGACTGA‑AGGATTATCCGCTTTATCCCTCCCTGGAATACCGCCAGATCACCTATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATg                                                                                                                  <  2:240342/160‑1 (MQ=255)                                                                                                                                                                                                                           atgatgCCTGGACTGA‑AGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATg                                                                                                                  >  1:240342/1‑160 (MQ=255)                                                                                                                                                                                                                             gatgCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCC‑ATGTAATTACTACTATGca                  >  1:358856/1‑249 (MQ=255)                                                                                                                                                                                                                              atgCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTG‑CCTTTTGC‑CCGGGAAAACCCCGGA‑CTAACGGAGCGCC‑ATGGAATTActactaagcaa                 >  1:2578/1‑245 (MQ=255)                                                                                                                                                                                                                                    ggACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑ATGTAATTACTACTATGCCAAATggg           >  1:603647/1‑249 (MQ=255)                                                                                                                                                                                                                                       cTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTCACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGCACCGGCGCGGCGTGAAGACTGGCGTTGCTTGTTTA‑CCTTTAGC‑CCGG‑GCAAGCCCGGGACCACCGCAGAGCA‑ATGTACTTACTAATATTCAAAACGGAacac       >  2:651609/1‑251 (MQ=255)                                                                                                                                                                                                                                          gaaagATTATCCACTTTATCCCTACCTGGAATCCCGTCAGATCACCGATGACCTTATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CTGG‑AAAAGCCCGGAACTACCGAAGCGCA‑ATGTAATTACTACTATGCAAAATGGAACAcc      <  1:471997/246‑1 (MQ=255)                                                                                                                                                                                                                                          a‑aGGATTATCCCCCTTATCCCTGCCTGGAATTCCGTCAGATCCCCGATGATCTGATGAAGCAACCGGCGGTGA‑CTGC‑CACTAACTTTGTTCGCTCTAACCCTACGCTTCCTCCCGCTCGCACGCTGCAATCTCGCTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTG‑TCTTTAGC‑CCGG‑AAAAGCCAGGAACTGCCGAAGCGCA‑ATGTAATTACTACTATGCCAAATggtccaccgt    <  2:1980/251‑9 (MQ=255)                                                                                                                                                                                                                                             agagtctCCGCTGTATCTCTACCTGGAATACCGACATATCACCGATGACCTGATGAATCATTCGGCGGTGA‑CTGG‑CACTAACCTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCATTCTCGTTTCGTCAAGGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGTCCGGAACTACCGAAGCGCA‑ATGTAATTACTACTATGCAAAATGGAACACCGGGc  <  2:181659/245‑1 (MQ=255)                                                                                                                                                                                                                                              gATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTA‑CCTTTAGA‑CCGG‑AAAAGACCGGGACTAACGAAGAGCACTTGTA‑TTAATACTATGCAAAAAGGAaccacgggc  >  1:410801/1‑243 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        AATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCA‑ATGTAATTACTACTATGCGAAATGGAACACCGGGC  >  NC_000913/4630677‑4631162

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