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%	180.8 / NA	52	G11G (GGA→GGG)	rpiA	ribose 5‑phosphate isomerase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/28);  total (24/28)

GAAGCATCTGAACTGGAAACGGCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCATCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCATCCATAAAAT  >  NC_000913/3059159‑3059424                                                                                                                                       |                                                                                                                                    gAAGCATCTGAACTGGAAACGGCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTa                                                                                                                                 >  2:881127/1‑139 (MQ=255) gAAGCATATGGACTGGAAACGGCCCCTTCAATCTTGCCATTCATTGTACCGAGCGAGTCCATAAAGTGTGCGGCGGTGGAACCTGTACCAACAACAACAATGGCGCCGGGCTGCACACAGTGAAGTTCCGCACAGCCTa                                                                                                                                 >  2:387267/1‑139 (MQ=255)                     ggCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCttttt                                                                                                                        <  1:260407/128‑1 (MQ=255)                     ggCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCttttt                                                                                                                        >  2:260407/1‑128 (MQ=255)                     ggCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATcc                                                                                                             >  1:399973/1‑139 (MQ=255)                                     ccTTTCATTGTACCGAGCGCATCAATAAAGTGAGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGttt                                                                                             >  2:552408/1‑139 (MQ=255)                                                    aGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGaa                                                                              >  2:379497/1‑139 (MQ=255)                                                     gcgcGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGaaa                                                                             >  1:309377/1‑139 (MQ=255)                                                     gcgcGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGaaa                                                                             >  1:865008/1‑139 (MQ=255)                                                             aTAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCGTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                     <  1:552408/139‑1 (MQ=255)                                                             aTAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACAGAAGTGCCGCCCACCCTACTGCTTTTTTCAATCCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                     <  2:653965/139‑1 (MQ=255)                                                              tAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACg                                                                    <  2:762546/139‑1 (MQ=255)                                                                aaGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCa                                                                  <  1:899336/139‑1 (MQ=255)                                                                   tgtgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCt                                                                                        <  1:343686/114‑1 (MQ=255)                                                                   tgtgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCt                                                                                        >  2:343686/1‑114 (MQ=255)                                                                    gtgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcaca                                                                      <  2:790353/131‑1 (MQ=255)                                                                    gtgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcaca                                                                      >  1:790353/1‑131 (MQ=255)                                                                     tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATtata                                                             >  2:669810/1‑139 (MQ=255)                                                                            gTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGc                                                                                          <  2:509036/103‑1 (MQ=255)                                                                            gTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGc                                                                                          >  1:509036/1‑103 (MQ=255)                                                                                aaCCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAg                                                  <  1:379497/139‑1 (MQ=255)                                                                                        ccTACACCAACAATGGTGACGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAg                                          <  1:625206/139‑1 (MQ=255)                                                                                             acCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATcc                                                                                                             >  1:21853/1‑67 (MQ=255)                                                                                             acCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATcc                                                                                                             <  2:21853/67‑1 (MQ=255)                                                                                             acCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCtg                                                                                       >  2:394465/1‑89 (MQ=255)                                                                                             acCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCtg                                                                                       <  1:394465/89‑1 (MQ=255)                                                                                                aacaaTGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCa                                                                  <  1:597123/107‑1 (MQ=255)                                                                                                aacaaTGCTGCCGGGCTGAACATACTGAAGTTCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGCTTCGCCTGTGGTATGAAATTTCACACGCa                                                                  >  2:597123/1‑107 (MQ=255)                                                                                                     tGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                     <  1:320064/99‑1 (MQ=255)                                                                                                     tGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                     >  2:320064/1‑99 (MQ=255)                                                                                                       gTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGc                           >  1:552540/1‑139 (MQ=255)                                                                                                        tGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATTGTTTCGCCTGTGGTATGAAATTTCacac                                                                     >  1:844478/1‑96 (MQ=255)                                                                                                        tGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                     <  2:844478/96‑1 (MQ=255)                                                                                                         gCCGGGCTGAACATACTGCACTGCCGCCCACCCTACTGATTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCcg                         <  2:865008/139‑1 (MQ=255)                                                                                                         gCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCt                                                                                                                            >  2:93208/1‑40 (MQ=255)                                                                                                         gCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCt                                                                                                                            <  1:93208/40‑1 (MQ=255)                                                                                                         gCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGAt                                                                                                  <  2:827023/66‑1 (MQ=255)                                                                                                         gCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGAt                                                                                                  >  1:827023/1‑66 (MQ=255)                                                                                                          ccGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCcgc                        <  1:669810/139‑1 (MQ=255)                                                                                                          ccGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCcg                         <  1:260661/138‑1 (MQ=255)                                                                                                          ccGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCcg                         >  2:260661/1‑138 (MQ=255)                                                                                                            gggCTGAACATACTGATGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGa                                         <  1:278319/120‑1 (MQ=255)                                                                                                            gggCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTAAGa                                         >  2:278319/1‑120 (MQ=255)                                                                                                             ggCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTg                     <  2:94061/139‑1 (MQ=255)                                                                                                                tGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGtt                  <  2:309377/139‑1 (MQ=255)                                                                                                                     aTACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTc             <  1:707875/139‑1 (MQ=255)                                                                                                                       aCTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCAt           <  1:881127/139‑1 (MQ=255)                                                                                                                       aCTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCAt           <  1:930694/139‑1 (MQ=255)                                                                                                                       aCTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCAt           <  2:399973/139‑1 (MQ=255)                                                                                                                       aCTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCAt           >  2:930694/1‑139 (MQ=255)                                                                                                                                ccgccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCAt                                                                 <  1:175682/76‑1 (MQ=255)                                                                                                                                ccgccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCAt                                                                 >  2:175682/1‑76 (MQ=255)                                                                                                                                ccgccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCATCCATAGAat  >  1:600192/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                    GAAGCATCTGAACTGGAAACGGCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCATCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCATCCATAAAAT  >  NC_000913/3059159‑3059424

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