Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,630,912	G→A	K60K (AAG→AAA)	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,912	0	G	A	91.4%	178.8 / 4.2	70	K60K (AAG→AAA)	slt	soluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base A (33/31);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGG  >  NC_000913/4630669‑4631160                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                              ttAGTATAAATCTAATCCGCGGCATGAAGCTGCATGGATGTATTTTCACTTAGAGGATGCGCTGGTGGAAAAAGCCAAAAAAGTTACCGGTCGGCTGTTGGCTGCAGGTGTCTGGCTGCGGACGGACAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGAGGTCGACCAATTGATGCCTGGACTGA‑GAGATTAt                                                                                                                                                                                                                                                         <  2:426374/250‑1 (MQ=255)   aGCATAAAGCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATcc                                                                                                                                                                                                                                                       <  2:114568/250‑1 (MQ=255)         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‑AAAAGCCCGGAACTACCGAAGCGCAATGTAAtt                           <  1:294110/248‑1 (MQ=255)                                                                                                                                                                                                                          gTCGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGACCTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AACAGACCGGAACTACCGAAGCGCACAGTAAtt                           >  1:51851/1‑250 (MQ=255)                                                                                                                                                                                                                            cGAACAAATGATGCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTTTTCT‑CCTTTATCACCGG‑AAAGACCCGGA‑CAACCCAAGCGCCATGTAAAtact                        >  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‑AAAAGCCCGGAACTACCGAAGCGCCATGTAATTACTACTATGca                >  1:358856/1‑249 (MQ=255)                                                                                                                                                                                                                                      atgCCTGGACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTG‑CCTTTTGC‑CCGGGAAAACCCCGGA‑CTAACGGAGCGCCATGGAATTActactaagcaa               >  1:2578/1‑245 (MQ=255)                                                                                                                                                                                                                                            ggACTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCCAAATggg         >  1:603647/1‑249 (MQ=255)                                                                                                                                                                                                                                               cTGA‑AAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGA‑CTGT‑CACTCACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGCACCGGCGCGGCGTGAAGACTGGCGTTGCTTGTTTA‑CCTTTAGC‑CCGG‑GCAAGCCCGGGACCACCGCAGAGCAATGTACTTACTAATATTCAAAACGGAacac     >  2:651609/1‑251 (MQ=255)                                                                                                                                                                                                                                                  a‑aGGATTATCCCCCTTATCCCTGCCTGGAATTCCGTCAGATCCCCGATGATCTGATGAAGCAACCGGCGGTGA‑CTGC‑CACTAACTTTGTTCGCTCTAACCCTACGCTTCCTCCCGCTCGCACGCTGCAATCTCGCTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTG‑TCTTTAGC‑CCGG‑AAAAGCCAGGAACTGCCGAAGCGCAATGTAATTACTACTATGCCAAATggtccaccgt  <  2:1980/251‑9 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                              TTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGA‑AGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGA‑CGGT‑CACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAG‑CCTTTAGC‑CCGG‑AAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGG  >  NC_000913/4630669‑4631160

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