Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	383,477	G→A	G153G (GGC→GGT)	yaiP ←	putative family 2 glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	383,477	0	G	A	100.0%	145.7 / NA	41	G153G (GGC→GGT)	yaiP	putative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/25);  total (16/25)

TAATAAACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACA  >  NC_000913/383380‑383573                                                                                                  |                                                                                                 taataaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCAcc                                                                                                <  1:917046/100‑1 (MQ=255)     aaaCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTc                                                                                            >  2:1241060/1‑100 (MQ=255)      aaCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCa                                                                                           >  2:1253531/1‑100 (MQ=255)       aCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCaa                                                                                          <  2:782969/100‑1 (MQ=255)                    ccggACCTGCTGTAGCGTGCGCTTCATAACAAACATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCAGCACCTTTCAACGCAGTAGAggg                                                                              >  2:138770/3‑99 (MQ=255)                    cAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAggg                                                                              <  2:605819/99‑1 (MQ=255)                      ggAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAAt                                                                           <  1:146937/100‑1 (MQ=255)                       gAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATg                                                                          >  2:32074/1‑100 (MQ=255)                           tgctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATgccgc                                                                      <  1:1075754/100‑1 (MQ=255)                                  gCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACgg                                                                <  2:343433/99‑1 (MQ=255)                                      gcgcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCg                                                           <  1:1308035/100‑1 (MQ=255)                                      gcgcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCg                                                           >  1:692425/1‑100 (MQ=255)                                        gcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGc                                                         <  2:1512874/100‑1 (MQ=255)                                           tCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCatc                                                      >  2:959821/1‑98 (MQ=255)                                            cATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcg                                                     >  2:877356/1‑100 (MQ=255)                                               aaCAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGTTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCgcgc                                                  <  2:858652/100‑1 (MQ=255)                                               aaCAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCgcgc                                                  <  1:1223395/100‑1 (MQ=255)                                                      atcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAAttt                                           <  1:1335459/100‑1 (MQ=255)                                                      atcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAAttt                                           <  2:279281/100‑1 (MQ=255)                                                       tcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTc                                          <  1:613050/100‑1 (MQ=255)                                                       tcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTc                                          >  2:1189002/1‑100 (MQ=255)                                                         ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                        >  2:907406/1‑100 (MQ=255)                                                         ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                        <  1:953141/100‑1 (MQ=255)                                                         ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                        <  2:203442/100‑1 (MQ=255)                                                           cgaCTTTACGGTCACGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGcc                                      >  1:1063493/3‑100 (MQ=255)                                                             aCTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGcccg                                    >  1:325234/1‑98 (MQ=255)                                                             aCTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAg                                    >  1:881183/1‑100 (MQ=255)                                                                 tACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCata                                <  1:760644/100‑1 (MQ=255)                                                                    ggTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAg                             <  1:474256/100‑1 (MQ=255)                                                                     gTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGc                            <  1:673927/100‑1 (MQ=255)                                                                     gTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGc                            >  1:605674/1‑100 (MQ=255)                                                                                gtgCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCt                 <  2:1433021/100‑1 (MQ=255)                                                                                 tgCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCtt                >  2:952876/1‑100 (MQ=255)                                                                                  gCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTg               <  1:133235/100‑1 (MQ=255)                                                                                  gCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTg               <  1:295399/100‑1 (MQ=255)                                                                                    ggTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGAt             <  1:1186043/100‑1 (MQ=255)                                                                                     gTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATc            <  1:576798/100‑1 (MQ=255)                                                                                      tAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCg           <  2:499911/100‑1 (MQ=255)                                                                                       aaCAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCgg          >  1:966874/1‑100 (MQ=255)                                                                                            aCCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGa     <  2:692423/100‑1 (MQ=255)                                                                                             cccGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGaa    >  2:211123/1‑100 (MQ=255)                                                                                               cgcaCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGTGCGGAACa  >  1:1423906/1‑100 (MQ=255)                                                                                                  |                                                                                                 TAATAAACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACA  >  NC_000913/383380‑383573

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