Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,911,877:1	+A	7.3%	coding (1773/2235 nt)	relA ←	(p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,911,876	1	.	A	7.3%	201.1 / 11.4	55	coding (1774/2235 nt)	relA	(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base . (30/21);  new base A (2/2);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

AATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCTACCACGCGGACCACCAGCGAATATCCGGCGGAGTAGCTCTCACCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCACGGTCCTGTTTACGGAA  >  NC_000913/2911592‑2912164                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                 aagggTCATGAAATCACGTACCAACCCATGTCGTTCATTAGCTACCACGCGGGCCACCAGCGAATATCCGGCGGAGTAGCTCTCACCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTg                                                                                                                                                                                                                                                                                              <  1:17514/286‑1 (MQ=255)                         cccACTACGATCATTAGCTACCACGCGGACCACCAGCGAATATCCGGCGGAGTAGCTCTCACCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCCTCTCCAGGAACCGGCTGGCAGCAGCCCGCGATGTGGTGCATCAGGTTGCCACCACCTTCGCCTACACCGCGCCCGTTCTCTTT‑ACTGCGGTTTCGCGGCgttgctcttttt                                                                                                                                                                                                                                                                      >  1:322146/1‑279 (MQ=255)                            aCTACGATCATTAGCTACCACGCGGACCACCAGCGAATATCCGGCGGAGTAGCTCTCACCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTAACTGCGGTTTTGCGGCGTGTAGCtttttt                                                                                                                                                                                                                                                                     >  1:107927/1‑288 (MQ=255)                            aCTACGATCATTAGCTACCACGCGGACCACCAGCGAATATCCGGCGGAGTAGCTCTCACCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTAACTGCGGTTTTGCGGCGTGTAGCTTTTTTg                                                                                                                                                                                                                                                                    >  1:481605/1‑289 (MQ=255)                                   tCATTAGCTACCACGCGGACCACCAGCGAATATCCGGCGGAGTAGCTCTCACCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGtt                                                                                                                                                                                                                                                            >  1:369795/1‑289 (MQ=255)                                            aCCACGCGGACCACCAGCGAATATCCGGCGGAGTAGCTCTCACCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAgc                                                                                                                                                                                                                                                    >  2:240196/1‑288 (MQ=255)                                                         ccggcgaagtgccGGCGGGGTAGCTGTCACCCCATCCCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAAGCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTg                                                                                                                                                                                                                                      <  1:240196/278‑1 (MQ=255)                                                                     cggcggAGTAGCTCTCACCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACt                                                                                                                                                                                                                          <  2:369795/289‑1 (MQ=255)                                                                                     aCCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATGAAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAAtt                                                                                                                                                                                                            >  2:281647/1‑288 (MQ=255)                                                                                     aCCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCCACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAATTTGCTTCACCGCGGCGGCGTCCTGCTCTTCGTCATTCGGCTTTTTAAAttcc                                                                                                                                                                                                          >  1:175825/1‑287 (MQ=255)                                                                                       cccATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGa                                                                                                                                                                                                        <  2:307433/289‑1 (MQ=255)                                                                                           tACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCGTCAGCGCGGCGGCGTCCTGCGCTTCGGCCCTCGGCTTATTAaattgcaattgc                                                                                                                                                                                                    >  2:309083/1‑289 (MQ=255)                                                                                            aCCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTAACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGc                                                                                                                                                                                                    <  2:107927/289‑1 (MQ=255)                                                                                            aCCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTTAACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGc                                                                                                                                                                                                    <  2:481605/289‑1 (MQ=255)                                                                                                cgTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCa                                                                                                                                                                                                   >  2:349657/1‑285 (MQ=255)                                                                                                cgTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCa                                                                                                                                                                                                   <  1:349657/285‑1 (MQ=255)                                                                                                 gTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTTGCTTTTTTGCTTAAGTTGCTTCAGCGCGGCGGCTTGCTGCTCTTCGGCACTCGGATTATTAAATTGCGATTGCAGGaat                                                                                                                                                                                              >  1:205666/1‑288 (MQ=255)                                                                                                              tCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGAt                                                                                                                                                                                 <  2:405680/289‑1 (MQ=255)                                                                                                                     cATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCa                                                                                                                                                                                                                            >  2:479654/1‑239 (MQ=255)                                                                                                                     cATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCa                                                                                                                                                                                                                            <  1:479654/239‑1 (MQ=255)                                                                                                                       tGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGCCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAACCGGCTAGCAGCAGCGCGCGATGTGGTTCATCAGGTTGCCAACACCTTCGACTACCACGCGTCCGTTCTCTTT‑ATTGCGGTTGTGCGGCGTGTACCTTTTTTGACGACGTTGCCTACATGCGGACGTGGCCCGCTTTCAGGCACGCGGCTGATTAACTCGCGCTTGCTGGCTGTTCACCATATCATTGACACg                                                                                                                                                                          >  1:232232/1‑287 (MQ=255)                                                                                                                        gggAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTg                                                                                                                                                                               <  2:495062/281‑1 (MQ=255)                                                                                                                        gggAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTg                                                                                                                                                                               >  1:495062/1‑281 (MQ=255)                                                                                                                                    ccGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGccacc                                                                                                                                                           >  2:377123/1‑289 (MQ=255)                                                                                                                                          gttgttCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCCCCAATcg                                                                                                                                                      >  1:351348/1‑288 (MQ=255)                                                                                                                                                            ggTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGa                                                                                                                                    <  2:351348/288‑1 (MQ=255)                                                                                                                                                              tgtACTGAAATACCGCGCCCCTGGGTAATGAAAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGAc                                                                                                                                   <  1:281647/288‑1 (MQ=255)                                                                                                                                                                 aCTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGAAAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGAcatcat                                                                                                                              <  2:175825/289‑1 (MQ=255)                                                                                                                                                                        tACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGtt                                                                                                                       >  1:399873/1‑289 (MQ=255)                                                                                                                                                                             cgcCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAAc                                                                                                                  >  1:200166/1‑289 (MQ=255)                                                                                                                                                                             cgcCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCCCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCACCGCCACCAATCGCCCCCAGCaagctgcctcactccttgaacatttaa                                                                                                                   >  1:347615/1‑263 (MQ=255)                                                                                                                                                                               ccccTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAAATCGTCGACATCATTGAAGTTGTAAcgc                                                                                                                >  2:53029/1‑289 (MQ=255)                                                                                                                                                                                          tG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTgag                                                                                                                                                                             <  1:344878/217‑1 (MQ=255)                                                                                                                                                                                          tG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTgag                                                                                                                                                                             >  2:344878/1‑217 (MQ=255)                                                                                                                                                                                                ccGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGCTATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCCTTGAAGTTGTAACGCGGCAGCAGATGtttt                                                                                                 >  1:187266/1‑288 (MQ=255)                                                                                                                                                                                                ccGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCCTTGAAGCTGTAACCAGGCAGCAGATGttttt                                                                                                >  1:88830/1‑289 (MQ=255)                                                                                                                                                                                                ccGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACCCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGtttt                                                                                                 >  2:145870/1‑288 (MQ=255)                                                                                                                                                                                                    cAATCTCATGCCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCAGCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGc                                                                                            <  1:377123/289‑1 (MQ=255)                                                                                                                                                                                                                    agtCGGCTGGCAGCAGCGCGCGATGGTGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTg                                                                               <  2:205666/284‑1 (MQ=255)                                                                                                                                                                                                                        ggCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCCGCAACTCGTCGAACTCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCCTTCAGGCTTGTcccc                                                                         >  1:96466/1‑288 (MQ=255)                                                                                                                                                                                                                        ggCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCa                                                                        >  1:114247/1‑289 (MQ=255)                                                                                                                                                                                                                                           tgtgGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGCTGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCgtcgtc                                                     >  1:239043/1‑289 (MQ=255)                                                                                                                                                                                                                                                            gccacacCTTCAACTCCCTCGCGACCTTTATCTTT‑ACTGCGGTTTGGGTGCTTGTAGCGTGTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCccag                                    <  1:309083/286‑1 (MQ=255)                                                                                                                                                                                                                                                               acacCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCccagcca                                 <  2:88830/289‑1 (MQ=255)                                                                                                                                                                                                                                                                   cTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATTATTGAAGTTGTAACGCGGCAGCAGCTGGTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGGTCCAGCTCGTCGTCAAGGATTCTCAGCCCAGCTAGaaa                             >  1:346704/1‑288 (MQ=255)                                                                                                                                                                                                                                                                           aCCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGtttt                                                                                                 <  2:287617/213‑1 (MQ=255)                                                                                                                                                                                                                                                                           aCCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGtttt                                                                                                 >  1:287617/1‑213 (MQ=255)                                                                                                                                                                                                                                                                           aCCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGt                     >  1:117859/1‑289 (MQ=255)                                                                                                                                                                                                                                                                             cACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCa                                                                                                                                                                                          <  2:147185/122‑1 (MQ=255)                                                                                                                                                                                                                                                                             cACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCa                                                                                                                                                                                          >  1:147185/1‑122 (MQ=255)                                                                                                                                                                                                                                                                             cACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTc                    >  2:446081/1‑288 (MQ=255)                                                                                                                                                                                                                                                                             cACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCa                   <  2:117859/289‑1 (MQ=255)                                                                                                                                                                                                                                                                             cACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCa                   <  2:399873/289‑1 (MQ=255)                                                                                                                                                                                                                                                                             cACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCa                   <  2:96466/289‑1 (MQ=255)                                                                                                                                                                                                                                                                               cgcgACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCACg                 >  2:145857/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                              t‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCACGGTCCTGTTTACGGaa  >  1:342146/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                              t‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCACGGTCCTGTTTACGGa   >  2:480629/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                 AATGGTCGTGATATCACGTAACAACCCACTACGATCATTAGCTACCACGCGGACCACCAGCGAATATCCGGCGGAGTAGCTCTCACCCCATACCGCGTCAACAATGCGTTCTGGCGCATGGGAGCGCAGTTCCGCCAGTTGTTCGCAATCGGCGCGGTGTACTGAAATACCGCGCCCCTGGGTAATG‑AAGCCGACAATCTCATCTCCAGGAATCGGCTGGCAGCAGCGCGCGATGTGGTGCATCAGGTTGCCAACACCTTCGACTACCACGCGACCGTTATCTTT‑ACTGCGGTTTTGCGGCGTGTAGCTTTTTTGCTGAAGTTGCTTCAGCGCGGCGGCGTCCTGCTCTTCGGCACTCGGCTTATTAAATTGCGATTGCAGGAAGTTCACCATCTGATTGAGACGGATATCCCCGCCACCAATCGCCGCCAGCAACTCGTCGACATCATTGAAGTTGTAACGCGGCAGCAGATGTTTTTCTGCTTCTTTCAGGCTGATCCCCAGATGTTCCAGCTCGTCGTCAAGGATTTGCCGCCCAGCCAGAATGTTTTTGTCACGGTCCTGTTTACGGAA  >  NC_000913/2911592‑2912164

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑