Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,146,011	C→T	*195* (TAG→TAA)	yceF ←	m(7)GTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,146,011	0	C	T	100.0%	120.6 / NA	38	*195* (TAG→TAA)	yceF	m(7)GTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/17);  total (21/17)

CTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCAC  >  NC_000913/1145882‑1146138                                                                                                                                  |                                                                                                                                cTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATc                                                                                                                        >  2:582248/1‑139 (MQ=255)                   tACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTttcccttcc                                                                                                      <  1:671566/139‑1 (MQ=255)                            ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAAc                                                                                             >  1:107848/1‑139 (MQ=255)                            ataAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAAc                                                                                             >  2:830656/1‑139 (MQ=255)                              aaaCCTGATAATGGGCGGTTGGCCTGATAAGGCATTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACAt                                                                                           >  1:2682/1‑139 (MQ=255)                                   tGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGac                                                                                      >  1:846980/1‑139 (MQ=255)                                    gATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                               <  2:743345/129‑1 (MQ=255)                                    gATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                               >  1:743345/1‑129 (MQ=255)                                    gATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGaca                                                                                     >  2:237978/1‑139 (MQ=255)                                     aTAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCGGCGATTGCCGGATGCGACGTAACCCCATACGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCAACGTAACCTCTGacac                                                                                    >  2:852859/1‑139 (MQ=255)                                     aTAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGacac                                                                                    <  1:573746/139‑1 (MQ=255)                                     aTAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGacac                                                                                    >  2:261220/1‑139 (MQ=255)                                         tGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAg                                                                                                                      <  1:379648/101‑1 (MQ=255)                                         tGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAg                                                                                                                      >  2:379648/1‑101 (MQ=255)                                           ggCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGcgcg                                                                              >  2:473371/1‑139 (MQ=255)                                                             gCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAaaccaacc                                                            <  2:107848/139‑1 (MQ=255)                                                              cGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAaccaacca                                                           <  1:830656/139‑1 (MQ=255)                                                              cGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGACTAGCGGTAAaccaacca                                                           <  1:582248/139‑1 (MQ=255)                                                               gTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAg                                                          <  2:102815/139‑1 (MQ=255)                                                                tttACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGc                                                         <  2:846980/139‑1 (MQ=255)                                                                       cgcATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAg                                                  >  1:896649/1‑139 (MQ=255)                                                                               gccgccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGc                                                         >  1:726917/1‑124 (MQ=255)                                                                               gccgccAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGc                                                         <  2:726917/124‑1 (MQ=255)                                                                                           gATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGc                              <  2:896649/139‑1 (MQ=255)                                                                                                         aCGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAaccaacca                                                           <  1:598110/96‑1 (MQ=255)                                                                                                         aCGTAACCACATCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAaccaacca                                                           >  2:598110/1‑96 (MQ=255)                                                                                                                  caTCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATcc       >  1:457216/1‑139 (MQ=255)                                                                                                                   aTCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACg                                                                                                 >  1:576438/1‑48 (MQ=255)                                                                                                                   aTCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACg                                                                                                 <  2:576438/48‑1 (MQ=255)                                                                                                                    tCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                               <  1:725153/49‑1 (MQ=255)                                                                                                                    tCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTa                                                                                               >  2:725153/1‑49 (MQ=255)                                                                                                                    tCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAgc                                                                                 >  2:338224/1‑63 (MQ=255)                                                                                                                    tCCGGCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAgc                                                                                 <  1:338224/63‑1 (MQ=255)                                                                                                                       ggCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGa                                                                                       <  1:753871/54‑1 (MQ=255)                                                                                                                       ggCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGa                                                                                       >  2:753871/1‑54 (MQ=255)                                                                                                                       ggCGCAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCAc  >  1:371862/1‑139 (MQ=255)                                                                                                                        gcgcAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAaccaacca                                                           <  1:458970/81‑1 (MQ=255)                                                                                                                        gcgcAAAAAATTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAaccaacca                                                           >  2:458970/1‑81 (MQ=255)                                                                                                                                  |                                                                                                                                CTGCGTAACGCGCGCTAATACAAGCATATAAACCTGATAATGGGCGGTTGGCCTGATAAGGCGTTTACGCCGCATCAGGCCGCCAGCACCGATTGCCGGATGCGACGTAACCACATCCGGCGCAAAAAACTATCCCATCAGCGGGTTTTTCCCTTCCCGACGTAACATCTGACACAGCGCGATTAGCGGTAAACCAACCAGCGTGTTAGGATCACGCCCCTCTAAACGCTCAAACAGCGTAATGCCAAATCCTTCAC  >  NC_000913/1145882‑1146138

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