Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,245,522	G→A	G100D (GGC→GAC)	exuT →	hexuronate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,245,522	0	G	A	100.0%	117.3 / NA	35	G100D (GGC→GAC)	exuT	hexuronate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/18);  total (17/18)

GCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGGCCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTCCATCGCAGTAGGCTACTTTAACG  >  NC_000913/3245389‑3245644                                                                                                                                      |                                                                                                                           gctgctTATACGGTCATGCAACCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTgg                                                                                                                       <  1:769861/139‑1 (MQ=255)   tgctTATACGGTCATGCAACCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCt                                                                                                                     >  1:314479/1‑139 (MQ=255)   tgctTATACGGTCATGCAACCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCt                                                                                                                     <  1:358325/139‑1 (MQ=255)         taCGGTCATGCAACCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCt                                                                                                               <  1:11982/139‑1 (MQ=255)               cATGCAACCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCgtggt                                                                                                         >  2:635785/1‑139 (MQ=255)               cATGCAACCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCgtggt                                                                                                         >  1:709001/1‑139 (MQ=255)                 tGCAACCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGTTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTgg                                                                                                                       >  2:381442/1‑123 (MQ=255)                 tGCAACCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGTTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTgg                                                                                                                       <  1:381442/123‑1 (MQ=255)                    aaCCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGcggt                                                                                                    <  1:51420/139‑1 (MQ=255)                        ggTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGcggtcggt                                                                                                >  1:669179/1‑139 (MQ=255)                                gTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGa                                                                                        <  1:635785/139‑1 (MQ=255)                                    tGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAagca                                                                                    <  2:709001/139‑1 (MQ=255)                                    tGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAagca                                                                                    >  2:340558/1‑139 (MQ=255)                                           gATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCgatga                                                                             >  2:755290/1‑139 (MQ=255)                                           gATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCgatga                                                                             >  2:701638/1‑139 (MQ=255)                                                           aaaTCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGaa                                                             >  1:37539/1‑139 (MQ=255)                                                              tCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCt                                                                                                                     <  1:114505/80‑1 (MQ=255)                                                              tCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCt                                                                                                                     >  2:114505/1‑80 (MQ=255)                                                              tCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCt                                                                                                               <  1:358906/86‑1 (MQ=255)                                                              tCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCt                                                                                                               >  2:358906/1‑86 (MQ=255)                                                              tCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGc                                                          <  1:630080/139‑1 (MQ=255)                                                                ggTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCa                                                        >  2:738604/1‑139 (MQ=255)                                                                            tttGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGt                                            >  2:45454/1‑139 (MQ=255)                                                                             ttGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGtt                                           <  1:340558/139‑1 (MQ=255)                                                                                            ccGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCg                            <  1:768065/139‑1 (MQ=255)                                                                                             cGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGt                           <  1:701638/139‑1 (MQ=255)                                                                                               tgtTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTc                         >  1:456397/1‑139 (MQ=255)                                                                                                gtTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTcc                        <  2:314479/139‑1 (MQ=255)                                                                                                      tggtgCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTCCATCGCa                  <  1:738604/139‑1 (MQ=255)                                                                                                          gCAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTCCATCGCAGTAg              <  1:755290/139‑1 (MQ=255)                                                                                                           cAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTCCATCGCAGTAgg             <  1:37598/139‑1 (MQ=255)                                                                                                           cAACCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTCCATCGCAGTAgg             >  1:109285/1‑139 (MQ=255)                                                                                                            aaCCGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTCCATCGCAGTAGGc            <  2:109285/139‑1 (MQ=255)                                                                                                              ccGCGCTGGCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTCCATCGCAGTAGGCTa          >  2:164743/1‑139 (MQ=255)                                                                                                                      gCAGGTAGCTGGGGTGACCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTCCATCGCAGTAGGCTACTTTAACg  >  1:517360/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                           GCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACTGGATGTGCTGGGTACGAAAATCGGTTATGCAATGTTTGCTGTGCTGTGGGCCGTGTTCTGTGGTGCAACCGCGCTGGCAGGTAGCTGGGGTGGCCTGGCTGTTGCTCGTGGTGCGGTCGGTGCCGCGGAAGCAGCGATGATCCCGGCGGGTCTGAAAGCCAGCTCCGAATGGTTCCCGGCGAAAGAGCGTTCCATCGCAGTAGGCTACTTTAACG  >  NC_000913/3245389‑3245644

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