Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,339,315	C→T	intergenic (‑155/‑15)	ribA ← / → pgpB	GTP cyclohydrolase 2/phosphatidylglycerophosphatase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,339,315	0	C	T	100.0%	43.6 / NA	13	intergenic (‑155/‑15)	ribA/pgpB	GTP cyclohydrolase 2/phosphatidylglycerophosphatase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/6);  total (7/6)

CAAATCCCACCATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCT  >  NC_000913/1339091‑1339523                                                                                                                                                                                                                                 |                                                                                                                                                                                                                 cAAATCCCACCATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGATTTATTAAAAAGGAGAGGCCATGCGTTCGAt                                                                                                                                                                                         <  2:178487/250‑1 (MQ=255)                   aaTCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTg                                                                                                                                                                      <  1:3130/251‑1 (MQ=255)                                                            ttAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAACCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGa                                                                                                                              <  1:191495/250‑1 (MQ=255)                                                                              tcCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGc                                                                                                            >  1:69193/1‑250 (MQ=255)                                                                               ccAGATAATGCTGGTTCTGTATTGGCTGATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACGGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCa                                                                                                           <  2:61288/250‑1 (MQ=255)                                                                                                               tGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAACGTTGGCTACTaaagacggc                                                                          >  2:9679/1‑245 (MQ=255)                                                                                                                        gagaAGGGGTTCACCTAATCACGACGAGTCGATCGTTCACGGAATGGCACAACTCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGACCCTGCGTTCGATTGACAGACGAACAGAAGTGGGAGCTGCACTATAGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGTTGGCGTTGGCAACCTGAAGAACAAAGTTGGCTACTAAAAGCGGATTTTTGGGt                                                                 >  1:248041/1‑251 (MQ=255)                                                                                                                                  tCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGaacct                                                       >  1:7984/1‑248 (MQ=255)                                                                                                                                         aTCACGACGCTTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGAGTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTGCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCAccc                                                <  1:152813/251‑1 (MQ=255)                                                                                                                                                                      gCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGAtt                   >  1:191183/1‑251 (MQ=255)                                                                                                                                                                         caAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGAtttt                 >  1:168424/1‑250 (MQ=255)                                                                                                                                                                              cTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCgg           >  1:7640/1‑251 (MQ=255)                                                                                                                                                                                        ttCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGTCTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGTGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATGTGATTTTATTCGGCTGGTTtct  <  2:166167/250‑1 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                 CAAATCCCACCATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCT  >  NC_000913/1339091‑1339523

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