Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,110,806	A→G	G234G (GGT→GGC)	tpiA ←	triosephosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,110,806	0	A	G	100.0%	139.2 / NA	39	G234G (GGT→GGC)	tpiA	triosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/20);  total (19/20)

TGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGATCACTTGTTCAGCGATGTTAGCGTCAACTTTAGCG  >  NC_000913/4110670‑4110926                                                                                                                                         |                                                                                                                         tgAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcgcc                                                                                                                        >  1:1206206/1‑139 (MQ=255)   agaTTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGcgccac                                                                                                                      >  2:1324344/1‑139 (MQ=255)      ttttGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCaccaa                                                                                                                   >  2:1363399/1‑139 (MQ=255)                                    taagtGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCaa                                                                                     >  1:1133853/1‑139 (MQ=255)                                        tGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAg                                                                                 >  1:1155624/1‑139 (MQ=255)                                                       ccGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGa                                                                  <  2:1155624/139‑1 (MQ=255)                                                        cGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGAc                                                                 <  2:1206206/139‑1 (MQ=255)                                                                 tCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACa                                                        >  2:1685290/1‑139 (MQ=255)                                                                  cAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACag                                                       <  1:1363399/139‑1 (MQ=255)                                                                       ttAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAgccg                                                  <  1:1685290/139‑1 (MQ=255)                                                                        tAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAgcgcg                                                                                                           <  1:815185/81‑1 (MQ=255)                                                                        tAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAgcgcg                                                                                                           >  2:815185/1‑81 (MQ=255)                                                                         aaGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAgccgcc                                                <  1:1688765/139‑1 (MQ=255)                                                                         aaGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAgccgcc                                                >  2:327819/1‑139 (MQ=255)                                                                           gCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAgcgcg                                                                                                           <  1:1323026/78‑1 (MQ=255)                                                                           gCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAgcgcg                                                                                                           >  2:1323026/1‑78 (MQ=255)                                                                               gTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTg                                          >  1:647085/1‑139 (MQ=255)                                                                                         ttCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAgcg                                                                                                             <  1:340785/62‑1 (MQ=255)                                                                                         ttCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAgcg                                                                                                             >  2:340785/1‑62 (MQ=255)                                                                                          tCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGc                                                                           <  2:1341208/95‑1 (MQ=255)                                                                                          tCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGc                                                                           >  1:1341208/1‑95 (MQ=255)                                                                                          tCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGATCACtt                               >  2:980300/1‑139 (MQ=255)                                                                                                gCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGATCACTTGTTCAg                         <  1:980300/139‑1 (MQ=255)                                                                                                gCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGATCACTTGTTCAg                         <  2:1133853/139‑1 (MQ=255)                                                                                                   ttAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATg                                      <  1:1608895/123‑1 (MQ=255)                                                                                                   ttAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATg                                      >  2:1608895/1‑123 (MQ=255)                                                                                                           tACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGAt                                    >  2:721654/1‑117 (MQ=255)                                                                                                           tACTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGAt                                    <  1:721654/117‑1 (MQ=255)                                                                                                             cTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAgcgcg                                                                                                           <  1:1408336/44‑1 (MQ=255)                                                                                                             cTGCGAAGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAgcgcg                                                                                                           >  2:1408336/1‑44 (MQ=255)                                                                                                                  aaGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAgcg                                                                                                             >  2:1618584/1‑37 (MQ=255)                                                                                                                  aaGGCGTCAGCTTTCAGAGAAGCGCCACCAACCAgcg                                                                                                             <  1:1618584/37‑1 (MQ=255)                                                                                                                    ggCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGt                                                                                                       <  1:321545/41‑1 (MQ=255)                                                                                                                    ggCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGt                                                                                                       >  2:321545/1‑41 (MQ=255)                                                                                                                     gCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATg                                      <  2:1232116/105‑1 (MQ=255)                                                                                                                     gCGTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATg                                      >  1:1232116/1‑105 (MQ=255)                                                                                                                       gTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGATCAc                                 >  1:671424/1‑108 (MQ=255)                                                                                                                       gTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGATCAc                                 <  2:671424/108‑1 (MQ=255)                                                                                                                       gTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGATCACTTGTTCAGCGATGTTAGCGTCAACTTTAGCg  >  1:160750/1‑139 (MQ=255)                                                                                                                       gTCAGCTTTCAGAGAAGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGATCACTTGTTCAGCGATGTTAGCGTCAACTTTAGCg  <  2:647085/139‑1 (MQ=255)                                                                                                                                      aGCGCCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATg                                      <  1:1304544/88‑1 (MQ=255)                                                                                                                                         |                                                                                                                         TGAGATTTTGCGTTAAGGCGAAGAGTTAAGGAAAGTAAGTGCCGGATATGAAATCCGGCACCTGTCAGACTTAAGCCTGTTTAGCCGCTTCTGCAGCTTTAACGATTACTGCGAAGGCGTCAGCTTTCAGAGAAGCACCACCAACCAGCGCGCCGTCGATATCCGGCTGAGCAAACAGTTCTGCAGCGTTAGACGCGTTTACAGAGCCGCCGTACTGAATGATCACTTGTTCAGCGATGTTAGCGTCAACTTTAGCG  >  NC_000913/4110670‑4110926

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