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

CGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGT  >  NC_000913/1287692‑1287960                                                                                                                                       |                                                                                                                                       cgcTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGt                                                                                                                                    >  2:859354/1‑139 (MQ=255)             ttACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGt                                                                                                                        >  1:1117846/1‑139 (MQ=255)             ttACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGt                                                                                                                        >  2:887704/1‑139 (MQ=255)             ttACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGAAGCGCTGTGTCAGTACTTTGTATAGt                                                                                                                        >  1:3879/1‑139 (MQ=255)                     agaAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACt                                                                                                                >  1:27154/1‑139 (MQ=255)                                  ggTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtc                                                                                                   >  1:1408727/1‑139 (MQ=255)                                      gCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCAc                                                                                                                    >  2:584477/1‑118 (MQ=255)                                      gCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCAc                                                                                                                    <  1:584477/118‑1 (MQ=255)                                       cGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGAc                                                                                              >  2:1425311/1‑139 (MQ=255)                                               tCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGAc                                                                                      <  2:678765/139‑1 (MQ=255)                                                cATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGAcc                                                                                     <  2:27154/139‑1 (MQ=255)                                                 aTCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCt                                                                                    >  2:205572/1‑139 (MQ=255)                                                  tCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACG‑TTTTCTCGACGTCACGACCTgc                                                                                  <  1:156680/139‑1 (MQ=255)                                                            aTTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGcatcat                                                                         <  1:859354/139‑1 (MQ=255)                                                                    tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  1:754172/139‑1 (MQ=255)                                                                    tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  2:1117846/139‑1 (MQ=255)                                                                    tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  2:701549/139‑1 (MQ=255)                                                                    tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  2:462983/139‑1 (MQ=255)                                                                    tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  1:1421485/139‑1 (MQ=255)                                                                    tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                                 <  1:205572/139‑1 (MQ=255)                                                                      aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtct                                                                                                  <  2:662043/104‑1 (MQ=255)                                                                      aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtct                                                                                                  >  1:662043/1‑104 (MQ=255)                                                                      aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCt                                                               >  1:716071/1‑139 (MQ=255)                                                                      aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCt                                                               >  2:661309/1‑139 (MQ=255)                                                                      aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCt                                                               >  2:200113/1‑139 (MQ=255)                                                                      aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCt                                                               >  1:14446/1‑139 (MQ=255)                                                                        cAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgt                                                             <  2:1205614/139‑1 (MQ=255)                                                                         aGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgtg                                                            <  1:1413568/139‑1 (MQ=255)                                                                         aGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGAAGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgtg                                                            <  2:3879/139‑1 (MQ=255)                                                                           ttATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      <  1:616816/63‑1 (MQ=255)                                                                           ttATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                                      >  2:616816/1‑63 (MQ=255)                                                                                   gCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGAt                                                  <  1:661309/139‑1 (MQ=255)                                                                                   gCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGAt                                                  >  2:166092/1‑139 (MQ=255)                                                                                    cTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATc                                                 <  2:14446/139‑1 (MQ=255)                                                                                       gCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGAc                                              <  1:166092/139‑1 (MQ=255)                                                                                                      tCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGGTGTGTAGGTATGATCGACATGAATa                                       >  1:1264121/1‑131 (MQ=255)                                                                                                      tCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGGTGTGTAGGTATGATCGACATGAATa                                       <  2:1264121/131‑1 (MQ=255)                                                                                                        gTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCa                             <  1:1010188/139‑1 (MQ=255)                                                                                                           cGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatga                          >  1:315717/1‑139 (MQ=255)                                                                                                           cGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatga                          >  1:1298088/1‑139 (MQ=255)                                                                                                             aTTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGTGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatgatg                        >  1:222963/1‑139 (MQ=255)                                                                                                                    tAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgc                                                                                <  1:965243/76‑1 (MQ=255)                                                                                                                    tAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgc                                                                                >  2:965243/1‑76 (MQ=255)                                                                                                                    tAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgtg                                                            >  2:73475/1‑96 (MQ=255)                                                                                                                    tAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgtg                                                            <  1:73475/96‑1 (MQ=255)                                                                                                                    tAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGc                 >  1:537116/1‑139 (MQ=255)                                                                                                                    tAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGc                 >  2:646436/1‑139 (MQ=255)                                                                                                                           aGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTc          <  2:1044725/139‑1 (MQ=255)                                                                                                                           aGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTc          >  1:681263/1‑139 (MQ=255)                                                                                                                            gACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTcc         >  1:973270/1‑139 (MQ=255)                                                                                                                            gACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTcc         >  2:1027134/1‑139 (MQ=255)                                                                                                                                cgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatg                           >  2:939162/1‑117 (MQ=255)                                                                                                                                cgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatg                           <  1:939162/117‑1 (MQ=255)                                                                                                                                cgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGAt     >  2:395423/1‑139 (MQ=255)                                                                                                                                   tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  <  2:973270/139‑1 (MQ=255)                                                                                                                                   tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  <  2:315717/139‑1 (MQ=255)                                                                                                                                   tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  <  2:1298088/139‑1 (MQ=255)                                                                                                                                   tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  <  1:887704/139‑1 (MQ=255)                                                                                                                                   tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  <  1:395423/139‑1 (MQ=255)                                                                                                                                   tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  <  1:200113/139‑1 (MQ=255)                                                                                                                                   tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  <  1:1425311/139‑1 (MQ=255)                                                                                                                                   tgtgTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGt  <  1:1027134/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       CGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGT  >  NC_000913/1287692‑1287960

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