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%	80.2 / NA	22	intergenic (‑155/‑15)	ribA/pgpB	GTP cyclohydrolase 2/phosphatidylglycerophosphatase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/11);  total (11/11)

CAGTTCTTCAAATCCCACCATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTCGCATTAAGGCTGCCT  >  NC_000913/1339083‑1339555                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 cAGTTCTTTAAATCCCACCAGCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCGTCTGCCACACGTTTAAGCTGCATGAAATTCTTCAGATAATTCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTTATCGCGACGCGTCGATCGTTCACGGAAGGGCACAAATCTGTCAATTTTTCGTGGAGCTGGCGGTTTCAGGTAGGATTGTGGGACTTGTTAAAAAGTAGAGGCCATGc                                                                                                                                                                                                                                <  2:292659/251‑1 (MQ=255)        tcAAATCCCACCATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTc                                                                                                                                                                                                                            <  2:102537/248‑1 (MQ=255)                 accaTCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCAATGCGTTCGATTGCCAGACg                                                                                                                                                                                                                >  1:239756/1‑251 (MQ=255)                       aGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGc                                                                                                                                                                                                          >  2:244563/1‑251 (MQ=255)                          aaaTCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCCCACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGt                                                                                                                                                                                                       >  2:169758/1‑251 (MQ=255)                                            ggCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCt                                                                                                                                                                                     <  1:25326/251‑1 (MQ=255)                                                            ccACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGcc                                                                                                                                                                     <  2:49146/251‑1 (MQ=255)                                                                acGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGAATTGGCTTATTTTGCCATAACGAGCCGAGCTTCACCTAATCACGCCGCGTCGATCATTCACGGACCGGCACAAATTTGTCATTTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCGGTGAGAGCTACCCTATTGCTTGTCATGCCAGTAGCCGTa                                                                                                                                                                  >  1:205086/1‑250 (MQ=255)                                                                                  aTTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGtt                                                                                                                                               <  1:33931/251‑1 (MQ=255)                                                                                       ccAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCa                                                                                                                                           <  1:244563/250‑1 (MQ=255)                                                                                                   ggTTCTGTAGTGGCTTAGTTGGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGACGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCATCCTGGAGAACaa                                                                                                                              <  2:267259/251‑1 (MQ=255)                                                                                                                  tATTTTGCCATAACTAGAAGGGGCTCACCTAATCACGACGCGTCGAGCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACGGGCGTTTTCAGTTATGATTGTGGGCCTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGTGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTaaa                                                                                                                <  1:305961/250‑1 (MQ=255)                                                                                                                                    aGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACt                                                                                             <  1:250642/251‑1 (MQ=255)                                                                                                                                      gggTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTACAAGCGGCTTTTTGGGTTACTg                                                                                            >  2:178429/1‑250 (MQ=255)                                                                                                                                         ttCACCTAATAAAGACGCGCAGATCGTTAACGCAATGGCCCAACTCTGTCAATTTTTCCTGGACCTGACGTTTTCCGTTCTGATTGTGGGACTTCTTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACAGCAGTGGGAGCTGCAATCTTGCTTGACATGCAAGGAGCAGTATGGATTTCTGGATGGCGTTGGCAACCTGGAGAACACATTTGGCTACTAAAAACGGCTTTTTGGGTTCCTGCAAc                                                                                        >  2:38930/1‑251 (MQ=255)                                                                                                                                             ccTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTTCCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTc                                                                                    <  1:225475/251‑1 (MQ=255)                                                                                                                                                                  gTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCAtt                                                               >  1:12168/1‑251 (MQ=255)                                                                                                                                                                                                        aCTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCt                         >  2:23389/1‑251 (MQ=255)                                                                                                                                                                                                                tttCAGTTATGAGTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCTATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTTTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCTCCCAGCCCTGGGGCGTCATTACACATTTGATGGTAGTCGGCTGGGTTCTCCGGTGTCTTCGtttg                  <  2:155789/250‑2 (MQ=255)                                                                                                                                                                                                                    aGTTATGATTGTGTGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGAACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACAGTAACCAAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTCGctt             >  2:265986/1‑249 (MQ=255)                                                                                                                                                                                                                     gTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTCGcctt            >  2:36272/1‑248 (MQ=255)                                                                                                                                                                                                                               tgGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTCGCATTAAGGCTGCCt  >  2:85518/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 CAGTTCTTCAAATCCCACCATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTCGCATTAAGGCTGCCT  >  NC_000913/1339083‑1339555

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