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

GGTTGCCAGTTCTTCAAATCCCACCATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTC  >  NC_000913/1339077‑1339540                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                  ggTTGCCAGTTCTTCAAATCCCACCATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAgg                                                                                                                                                                                                                       >  1:117535/1‑251 (MQ=255)     gCCAGTTCTTCAAATCCCACCAGCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCAt                                                                                                                                                                                                                   <  2:149375/251‑1 (MQ=255)                                  aTCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTgg                                                                                                                                                                                      >  1:90409/1‑251 (MQ=255)                                        ccATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTg                                                                                                                                                                                <  2:183405/251‑1 (MQ=255)                                                       tttGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTc                                                                                                                                                                  >  1:421455/1‑250 (MQ=255)                                                             ttCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGcc                                                                                                                                                             >  1:338850/1‑249 (MQ=255)                                                                  ccACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGCGCTGCACTATTGCTTGTCATGCCAGTAGcc                                                                                                                                                      >  2:155518/1‑251 (MQ=255)                                                                         tttAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATgg                                                                                                                                                >  2:192072/1‑250 (MQ=255)                                                                           tAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGAtt                                                                                                                                             >  1:71091/1‑251 (MQ=255)                                                                                 gCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTctggc                                                                                                                                       <  2:117535/251‑1 (MQ=255)                                                                                                gATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCt                                                                                                                        <  1:194924/251‑1 (MQ=255)                                                                                                   aaTGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGTCAACCTGga                                                                                                                     <  1:87986/251‑1 (MQ=255)                                                                                                                 aTTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGct                                                                                                       <  2:99920/251‑1 (MQ=255)                                                                                                                        tATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTaaa                                                                                                 <  2:412136/250‑1 (MQ=255)                                                                                                                          ttttGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGc                                                                                              >  2:17719/1‑251 (MQ=255)                                                                                                                           tttGCCATAACGAGAAGGGGTTCTCCTAATCTCGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGc                                                                                              <  1:106088/250‑1 (MQ=255)                                                                                                                                cATAACGAGAAGGGGATCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCttt                                                                                        >  1:322718/1‑251 (MQ=255)                                                                                                                                 aTAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCtttt                                                                                       >  1:305487/1‑251 (MQ=255)                                                                                                                                 aTAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGCACATGCGTTCGATTGACAGACGAACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGTTTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTCCTCAAAGCGGCtttt                                                                                       >  2:33961/1‑251 (MQ=255)                                                                                                                                                                    gATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGt                                                    >  2:175659/1‑251 (MQ=255)                                                                                                                                                                                     cacaAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGTGTTGGCAACCTGGAGAACAAAGTTGGCTTCTAAAAGCGGCTTTTTTGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGAttt                                   <  1:271740/251‑1 (MQ=255)                                                                                                                                                                                     cacaAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGCCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGAttt                                   >  1:156527/1‑251 (MQ=255)                                                                                                                                                                                          atcggTCAATTTTTCCTGGAATTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGGTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCTGCTTTTTGGGTTACTTAAACTGTCACCCAGCCCTGGGGCGTCATTGCACATTTGATTTTATTc                              <  2:185304/247‑1 (MQ=255)                                                                                                                                                                                              gTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCt                          <  1:155518/251‑1 (MQ=255)                                                                                                                                                                                                    ttttCCTGGAACTGGCGTTTTCAGTTAGGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGttt                     <  2:90409/250‑1 (MQ=255)                                                                                                                                                                                                                 ggCGTTTTCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGGTGGCAACCGGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGc        <  1:77798/250‑1 (MQ=255)                                                                                                                                                                                                                       ttCAGTTATGATTGTGGGACTTATTAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTc  <  1:334370/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                  GGTTGCCAGTTCTTCAAATCCCACCATCAGGAAATCGCCCCATGGGGTTGGCAGTTTGGCTTCTGCCACACGTTTAAGCTGCATGAAATTCTCCAGATAATGCTGGTTCTGTATTGGCTTATTTTGCCATAACGAGAAGGGGTTCACCTAATCACGACGCGTCGATCGTTCACGGAATGGCACAAATCTGTCAATTTTTCCTGGAACTGGCGTTTTCAGTTATGATTGTGGGACTTATCAAAAAGGAGAGGCCATGCGTTCGATTGCCAGACGTACCGCAGTGGGAGCTGCACTATTGCTTGTCATGCCAGTAGCCGTATGGATTTCTGGCTGGCGTTGGCAACCTGGAGAACAAAGTTGGCTACTAAAAGCGGCTTTTTGGGTTACTGAAACTGTCACCCAGCCCTGGGGCGTCATTACACATTTGATTTTATTCGGCTGGTTTCTCTGGTGTCTGCGTTTTC  >  NC_000913/1339077‑1339540

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