Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,931,057	T→C	I59I (ATT→ATC)	purT →	phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,931,057	0	T	C	100.0%	81.3 / NA	24	I59I (ATT→ATC)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/15);  total (9/15)

GGCGTAGAGGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATCGAA  >  NC_000913/1930983‑1931132                                                                           |                                                                            ggCGTAGAGGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAat                                                                          <  2:306795/78‑1 (MQ=255)  gCGTAGAGGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGc                                                                      <  1:306796/81‑1 (MQ=255)      agagGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGa                                                                  >  2:306833/1‑81 (MQ=255)         ggTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATgg                                                               >  2:306843/1‑81 (MQ=255)                gCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGc                                                         <  2:306797/80‑1 (MQ=255)                 ccGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCAt                                                       <  1:306798/81‑1 (MQ=255)                   gTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTa                                                     >  1:306834/1‑81 (MQ=255)                    tCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTAc                                                    <  1:306799/81‑1 (MQ=255)                       aTCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGcc                                                 >  2:306844/1‑81 (MQ=255)                         cGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCgt                                               >  2:306842/1‑81 (MQ=255)                          gCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCgt                                               <  2:306800/80‑1 (MQ=255)                                ccGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTg                                        <  2:306801/81‑1 (MQ=255)                                     gCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTg                                   <  1:306802/81‑1 (MQ=255)                                            ccATGCATGTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAAc                            >  2:306845/1‑81 (MQ=255)                                             catgcatgTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAAc                            >  2:306835/1‑80 (MQ=255)                                                  atgTCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACAt                       >  2:306841/1‑80 (MQ=255)                                                     tCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATtata                   <  1:306804/81‑1 (MQ=255)                                                     tCGCGCATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATtat                    <  2:306803/80‑1 (MQ=255)                                                       gcgcATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATc                 <  1:306805/81‑1 (MQ=255)                                                         gcATCGCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGt               <  1:306806/81‑1 (MQ=255)                                                              gCTCCCATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCgg          <  2:306807/81‑1 (MQ=255)                                                                  ccATGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGgag        <  1:306808/79‑1 (MQ=255)                                                                    aTGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATc     >  1:306848/1‑80 (MQ=255)                                                                    aTGTCATCAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATCg    <  2:306809/81‑1 (MQ=255)                                                                       tcatcaATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATCGaa  <  2:306810/80‑1 (MQ=255)                                                                           |                                                                            GGCGTAGAGGTGATTGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATCGAA  >  NC_000913/1930983‑1931132

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