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%	180.5 / NA	53	A267T (GCA→ACA)	purU	formyltetrahydrofolate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/29);  total (24/29)

GCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGT  >  NC_000913/1287690‑1287960                                                                                                                                         |                                                                                                                                       gccgcTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      <  1:1027792/139‑1 (MQ=255) gccgcTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      <  1:1277220/139‑1 (MQ=255)    gcTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTa                                                                                                                                   >  2:960021/1‑139 (MQ=255)      tGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACt                                                                                                                                 <  2:73510/139‑1 (MQ=255)                  cGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCa                                                                                                                     >  1:137910/1‑139 (MQ=255)                       agaAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      >  2:1403179/1‑117 (MQ=255)                                    ggTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtc                                                                                                   >  1:1453523/1‑139 (MQ=255)                                               tttCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      >  1:114881/1‑93 (MQ=255)                                               tttCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGt                                                                                                                                    >  2:199043/1‑95 (MQ=255)                                               tttCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGt                                                                                                                                    <  1:199043/95‑1 (MQ=255)                                                 tCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGAc                                                                                      <  1:1313972/139‑1 (MQ=255)                                                        tAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGc                                                                                                                      <  1:1032987/100‑1 (MQ=255)                                                        tAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGc                                                                                                                      >  2:1032987/1‑100 (MQ=255)                                                        tAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgcg                                                                               >  2:1040038/1‑139 (MQ=255)                                                        tAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCG‑AAACAAAGACGCGCTGTGTCa                                                                                                                                      <  1:61993/83‑1 (MQ=255)                                                            ggATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatc                                                                           <  2:173934/139‑1 (MQ=255)                                                              aTTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatcat                                                                         >  1:979094/1‑139 (MQ=255)                                                               ttAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCata                                                                        >  2:392735/1‑139 (MQ=255)                                                                    ggTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTc                                                                   <  2:1453523/139‑1 (MQ=255)                                                                      tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  2:979094/139‑1 (MQ=255)                                                                      tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  1:1040038/139‑1 (MQ=255)                                                                      tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  1:535289/139‑1 (MQ=255)                                                                      tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  2:344071/139‑1 (MQ=255)                                                                      tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  2:120149/139‑1 (MQ=255)                                                                        aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      >  2:737210/1‑68 (MQ=255)                                                                        aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGta                                                                                                                            <  2:250333/78‑1 (MQ=255)                                                                        aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGta                                                                                                                            >  1:250333/1‑78 (MQ=255)                                                                        aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCt                                                               <  1:495431/139‑1 (MQ=255)                                                                              tATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAg                                                         >  2:423285/1‑139 (MQ=255)                                                                                    tGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGa                                                   <  2:557321/139‑1 (MQ=255)                                                                                        ggCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      >  1:345667/1‑52 (MQ=255)                                                                                        ggCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      <  2:345667/52‑1 (MQ=255)                                                                                         gCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCa                                                                                                                     <  1:520907/68‑1 (MQ=255)                                                                                         gCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCa                                                                                                                     >  2:520907/1‑68 (MQ=255)                                                                                                     taatCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTc                                  <  1:392735/139‑1 (MQ=255)                                                                                                           ttCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCat                            >  1:1324491/1‑139 (MQ=255)                                                                                                           ttCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCat                            >  1:1416017/1‑139 (MQ=255)                                                                                                            tCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  1:355591/101‑1 (MQ=255)                                                                                                            tCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 >  2:355591/1‑101 (MQ=255)                                                                                                            tCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatg                           <  1:310349/139‑1 (MQ=255)                                                                                                            tCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatg                           <  1:960021/139‑1 (MQ=255)                                                                                                             cGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatga                          >  2:785025/1‑139 (MQ=255)                                                                                                                             aGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatc                                                                           >  2:1317653/1‑74 (MQ=255)                                                                                                                             aGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatc                                                                           <  1:1317653/74‑1 (MQ=255)                                                                                                                              gACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcg                                                                                 >  2:235815/1‑67 (MQ=255)                                                                                                                              gACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcg                                                                                 <  1:235815/67‑1 (MQ=255)                                                                                                                              gACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  2:532616/83‑1 (MQ=255)                                                                                                                              gACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 >  1:532616/1‑83 (MQ=255)                                                                                                                                  cgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACAt                                            >  2:584473/1‑100 (MQ=255)                                                                                                                                  cgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACAt                                            <  1:584473/100‑1 (MQ=255)                                                                                                                                   gcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcg                                                                                 >  2:850055/1‑62 (MQ=255)                                                                                                                                   gcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcg                                                                                 <  1:850055/62‑1 (MQ=255)                                                                                                                                     tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  <  2:137910/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       GCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGT  >  NC_000913/1287690‑1287960

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