Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,405,622	C→T	D312D (GAC→GAT)	purA →	adenylosuccinate synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,405,622	0	C	T	100.0%	51.3 / NA	19	D312D (GAC→GAT)	purA	adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/8);  total (11/8)

CCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAA  >  NC_000913/4405417‑4405846                                                                                                                                                                                                              |                                                                                                                                                                                                                                 cccctgctgGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGATACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGc                                                                                                                                                                                     <  1:307447/248‑1 (MQ=255)          gtggcgtggcGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGTGGCGTCGTCGTCTTACCGGCTGGCTGGATACCTTTGCCGTTCGGCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTg                                                                                                                                                                            <  1:229959/251‑1 (MQ=255)                   cgACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGATACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTg                                                                                                                                                                   >  1:97272/1‑251 (MQ=255)                     aCCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTGCTGTGTATCCTCAAAGCTCACTCCACTCGTGTAGGTGCAGGTCCTTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCACGGTAACGAATTCGGCGCAACTGCTGGCCGTCGTCGTCGTACCGGCTGGCTGGATACCGTGGCCGTTCGTCGTGCg                                                                                                                                                                                                             <  1:330949/207‑1 (MQ=255)                      ccGGTTCCGGCCTGGGCCCGCGTTATNTTGATTNNGTNCTGGGTATCCTCAAAGCTTNCTNNNCNNNGNTCGGTNCANNTCCGTTCCCNACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTNGGCGCAACTNNGGGGCGGCGTCGTCGTACCGGNGGGNNGNNTNCCgnggn                                                                                                                                                                                                                          >  2:77591/1‑189 (MQ=255)                      ccGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGGTCCTCTGCAAGCAGGGTAACGAAGTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGATACCGTTGc                                                                                                                                                                                                                          <  1:77591/193‑1 (MQ=255)                                     gCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGATACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTg                                                                                                                                                 >  1:237786/1‑251 (MQ=255)                                      cccGCGTTATGTTGATTACGTTCTGGGTATCTTCAAAGCTTCCTCCACTCGTGTCGGTGCAGGTACGTTCCCGACCGAACTGTTTGATAAAACTGGCGAGTTCCTCTGCAAGAATGATAACGAATTCGGCGCAACTACGGGGCGTAGGCGTCGTACCGACTGGCTGAATACCGTTGCCGTTCGTCGTGCGGTACCGCTGAACTCCCTTTCTGGCTTCTGCCGGACGCAACTGGACGGTCTGTATGGcatg                                                                                                                                                 >  1:271326/1‑247 (MQ=255)                                                                                cTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAAGTTTTTGATGAACATGGAGAGATCCTTTGCAAGCCGGGTAACGAATTCGGGGCAACTAAGGGGCGTCGTCGTAGTACCGGCTGGATGGGTACCGTTGGCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGGTTCTGCCTGACTAAACTGGAAGTTCTGGGTGGCCTGAAAGAGGTTACACTCTGCGTGGCTTACCGTATGCCGGATGGTc                                                                                                      >  1:279306/1‑251 (MQ=255)                                                                                                   ggTCCGTTCCCGACCGAACTGTTTGATGAAACTGNCTANGTCCNCTNCAAGCAGGGTAACGANTTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGANNCNNNTNCCGTtcgt                                                                                                                                                                                                                   <  2:233660/123‑1 (MQ=255)                                                                                                   ggTCCGTTCCCGACCGAACTGTTGGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACGACGGGGCGTCGTCGTCGTACCGGCTGGCTGGATACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAagan                                                                                                                                           <  1:27854/195‑2 (MQ=255)                                                                                                   ggTCCGTTCCCGACAGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTAGTCGTGGTACCGGCTGGCTGGATACCGTTGCCGTtcgt                                                                                                                                                                                                                   >  1:233660/1‑123 (MQ=255)                                                                                                          tCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGATACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCTTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTggn                                                                             <  1:62103/250‑2 (MQ=255)                                                                                                                  gaACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGATACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTg                                                                                              >  1:6324/1‑225 (MQ=255)                                                                                                                             tgaAACTGGCGAGTTACTCTGCAAGCAGGGCAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGATAACGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTg                                                                                                                            >  1:271401/1‑184 (MQ=255)                                                                                                                                                             aCGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGATACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGCGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGACAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTgg                         >  1:134135/1‑251 (MQ=255)                                                                                                                                                                                ggggCGTCGTCGTCGTACCGGCTGGCTGGATACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCGGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATg                                                                                                         <  1:48765/152‑1 (MQ=255)                                                                                                                                                                                     gtcgtcgtcgtACCGGCTGGCTGGATACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGaa  >  1:87449/1‑250 (MQ=255)                                                                                                                                                                                      tcgtcnncgnNNCGNNTGGCTGGATACCGTTGCCGTTCGNCGTGCGGTACAGCTGGACTCCCTGTCTGGNTNCNGCCTNNCTANACTGGACGTTNTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGc                                                                                                                          >  2:239620/1‑129 (MQ=255)                                                                                                                                                                                      tcgtcgtcgtACCGGCTGGCTGGATACCGTTGCCGTTCGGCGTGCTGTACAGCTGAACTCCCTTTCTTGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACGCTGCGTGGc                                                                                                                          <  1:239620/129‑1 (MQ=255)                                                                                                                                                                                                              |                                                                                                                                                                                                                                 CCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGATTTACGAAACCATGCCGGGCTGGTCTGAATCCACCTTCGGCGTGAA  >  NC_000913/4405417‑4405846

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