Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,287,826	C→T	A267T (GCA→ACA)	purU ←	formyltetrahydrofolate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,826	0	C	T	100.0%	114.1 / NA	34	A267T (GCA→ACA)	purU	formyltetrahydrofolate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/16);  total (18/16)

TTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGT  >  NC_000913/1287704‑1287960                                                                                                                           |                                                                                                                                       ttACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGt                                                                                                                        >  2:526914/1‑139 (MQ=255)  tACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTg                                                                                                                       >  2:789073/1‑139 (MQ=255)         agaAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACt                                                                                                                                 >  2:413682/1‑122 (MQ=255)         agaAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACt                                                                                                                                 <  1:413682/122‑1 (MQ=255)         agaAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACt                                                                                                                >  2:164061/1‑139 (MQ=255)         agaAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACt                                                                                                                >  2:47019/1‑139 (MQ=255)                       gTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      >  2:180630/1‑103 (MQ=255)                       gTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      <  1:180630/103‑1 (MQ=255)                                                aTTAAAGGGAACCAGTTATGTTTGCAGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatcat                                                                         <  1:333063/139‑1 (MQ=255)                                                  tAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCatat                                                                       >  1:769654/1‑139 (MQ=255)                                                        tAACCGGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  1:278621/139‑1 (MQ=255)                                                        tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCATGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  1:268039/139‑1 (MQ=255)                                                        tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  1:789073/139‑1 (MQ=255)                                                                      tGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  1:610768/125‑1 (MQ=255)                                                                      tGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 >  2:610768/1‑125 (MQ=255)                                                                      tGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGa                                                   <  1:643815/139‑1 (MQ=255)                                                                      tGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGa                                                   <  1:526914/139‑1 (MQ=255)                                                                        cTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATc                                                 <  2:769654/139‑1 (MQ=255)                                                                         tGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  2:461639/122‑1 (MQ=255)                                                                         tGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 >  1:461639/1‑122 (MQ=255)                                                                                             ttCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCat                            >  1:353787/1‑139 (MQ=255)                                                                                              tCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatg                           <  2:353787/139‑1 (MQ=255)                                                                                               cGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatga                          >  1:144234/1‑139 (MQ=255)                                                                                               cGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatga                          >  1:82235/1‑139 (MQ=255)                                                                                               cGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatga                          >  1:580653/1‑139 (MQ=255)                                                                                                        tAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGc                 <  1:79343/139‑1 (MQ=255)                                                                                                        tAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGc                 >  1:621759/1‑139 (MQ=255)                                                                                                        tAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGAACTGCGCGCATCATATCTTAAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGc                 >  1:783247/1‑139 (MQ=255)                                                                                                                   gcgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGa      <  2:82235/139‑1 (MQ=255)                                                                                                                    cgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTc                                  <  2:130860/110‑1 (MQ=255)                                                                                                                    cgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTc                                  >  1:130860/1‑110 (MQ=255)                                                                                                                    cgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGAt     >  2:160192/1‑139 (MQ=255)                                                                                                                       tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  >  1:768999/1‑139 (MQ=255)                                                                                                                       tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  <  1:42274/139‑1 (MQ=255)                                                                                                                        gtgtCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgtg                                                            >  1:632866/1‑80 (MQ=255)                                                                                                                           |                                                                                                                                       TTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGT  >  NC_000913/1287704‑1287960

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