Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,059,293	T→C	G11G (GGA→GGG)	rpiA ←	ribose 5‑phosphate isomerase, constitutive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,059,293	0	T	C	100.0%	131.2 / NA	36	G11G (GGA→GGG)	rpiA	ribose 5‑phosphate isomerase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/16);  total (20/16)

GAAACGGCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCATCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCATCCATA  >  NC_000913/3059174‑3059420                                                                                                                        |                                                                                                                                gAAACGGCCCCTTCAATCCGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAAtt                                                                                                              >  1:444557/1‑139 (MQ=255)      ggCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATcc                                                                                                         >  1:35349/1‑139 (MQ=255)      ggCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATcc                                                                                                         >  2:596745/1‑139 (MQ=255)                    ggCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGt                                                                                           >  1:394655/1‑139 (MQ=255)                            aTTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCtg                                                                                   <  1:487218/139‑1 (MQ=255)                                   cGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATg                                                                            >  1:709358/1‑139 (MQ=255)                                    gAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGa                                                                           >  1:678499/1‑139 (MQ=255)                                     aGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAAt                                                                                                               >  1:682568/1‑102 (MQ=255)                                     aGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAAt                                                                                                               <  2:682568/102‑1 (MQ=255)                                     aGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGaa                                                                          >  2:706281/1‑139 (MQ=255)                                        gggTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAAtt                                                                       >  1:15108/3‑139 (MQ=255)                                            cAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcac                                                                   <  1:387892/139‑1 (MQ=255)                                             aaTAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGCTTCGCCTGTGGTATGAAATTTcaca                                                                  <  1:596745/139‑1 (MQ=255)                                              aTAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                 <  2:709358/139‑1 (MQ=255)                                              aTAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                 <  2:576133/139‑1 (MQ=255)                                              aTAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                 <  2:466018/139‑1 (MQ=255)                                              aTAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                 <  2:98306/139‑1 (MQ=255)                                                   gtgtGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCAtt                                                            <  1:118534/139‑1 (MQ=255)                                                      tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCtg                                                                                   <  2:237845/113‑1 (MQ=255)                                                      tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCtg                                                                                   >  1:237845/1‑113 (MQ=255)                                                        cggcggTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTa                                                                                                                             >  1:105901/1‑69 (MQ=255)                                                        cggcggTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTa                                                                                                                             <  2:105901/69‑1 (MQ=255)                                                         ggcggTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATAc                                                      >  1:490253/1‑139 (MQ=255)                                                            ggTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACaaa                                                   <  1:706281/139‑1 (MQ=255)                                                                                    aaTGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCt                                                                                                        <  1:728123/62‑1 (MQ=255)                                                                                    aaTGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCt                                                                                                        >  2:728123/1‑62 (MQ=255)                                                                                    aaTGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTg                                                                                                       <  2:482530/63‑1 (MQ=255)                                                                                    aaTGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTg                                                                                                       >  1:482530/1‑63 (MQ=255)                                                                                        gTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGc                       >  1:629782/1‑139 (MQ=255)                                                                                        gTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGc                       >  1:399211/1‑139 (MQ=255)                                                                                                   aaCATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTaa            >  2:274557/1‑139 (MQ=255)                                                                                                             aGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCATCCATa  >  2:434185/1‑139 (MQ=255)                                                                                                               tGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACaaa                                                   >  1:417057/1‑88 (MQ=255)                                                                                                               tGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACaaa                                                   <  2:417057/88‑1 (MQ=255)                                                                                                                   gccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                 <  2:680992/70‑1 (MQ=255)                                                                                                                   gccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                 >  1:680992/1‑70 (MQ=255)                                                                                                                        |                                                                                                                                GAAACGGCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCATCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCATCCATA  >  NC_000913/3059174‑3059420

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