Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	383,259	C→T	R226Q (CGG→CAG)	yaiP ←	putative family 2 glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	383,259	0	C	T	100.0%	203.5 / NA	62	R226Q (CGG→CAG)	yaiP	putative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/32);  total (30/32)

ACCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAA  >  NC_000913/383123‑383386                                                                                                                                         |                                                                                                                                aCCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGa                                                                                                                               >  2:904752/1‑139 (MQ=255)  ccGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGaa                                                                                                                              <  1:904752/139‑1 (MQ=255)                       cTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                         <  2:971968/139‑1 (MQ=255)                                  gcgcATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGc                                                                                                              <  2:254261/123‑1 (MQ=255)                                  gcgcATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGc                                                                                                              >  1:254261/1‑123 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGc                                                                                                              >  2:980888/1‑113 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGc                                                                                                              <  1:980888/113‑1 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGt                                                                                    >  2:89222/1‑139 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGt                                                                                    >  1:965906/1‑139 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGG‑ATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTc                                                                                   >  2:153953/1‑139 (MQ=255)                                              cgcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                         <  2:116991/116‑1 (MQ=255)                                              cgcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                         >  1:116991/1‑116 (MQ=255)                                               gcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                 <  1:548991/139‑1 (MQ=255)                                                cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGt                                                                                               <  1:260125/124‑1 (MQ=255)                                                cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGt                                                                                               >  2:260125/1‑124 (MQ=255)                                                         caATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGt                                                                       <  2:48713/139‑1 (MQ=255)                                                            tCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCg                                                                                                                   >  2:202659/1‑92 (MQ=255)                                                            tCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCg                                                                                                                   <  1:202659/92‑1 (MQ=255)                                                             ccAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAACGTCTTCGACTTTAGTACGa                                                                   >  1:23391/1‑139 (MQ=255)                                                               aaCGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATc                                                                 >  2:287441/1‑139 (MQ=255)                                                                 cgccgcCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                            >  1:33067/1‑126 (MQ=255)                                                                 cgccgcCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                            <  2:33067/126‑1 (MQ=255)                                                                    cgccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGt                                                                                                                        <  2:460091/79‑1 (MQ=255)                                                                    cgccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGt                                                                                                                        >  1:460091/1‑79 (MQ=255)                                                                    cgccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCgaga                                                            >  1:569143/1‑139 (MQ=255)                                                                      ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGtt                                                                                                                 >  2:592776/1‑84 (MQ=255)                                                                      ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGtt                                                                                                                 <  1:592776/84‑1 (MQ=255)                                                                      ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGaaa                                                          >  2:64096/1‑139 (MQ=255)                                                                                 aCTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                               <  2:415631/139‑1 (MQ=255)                                                                                 aCTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                               <  2:1060105/139‑1 (MQ=255)                                                                                    cctcACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                            <  2:70605/139‑1 (MQ=255)                                                                                     ctcACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                            <  2:497792/106‑1 (MQ=255)                                                                                     ctcACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                            >  1:497792/1‑106 (MQ=255)                                                                                         cgcgGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTGCGATCCGAGAAACCGAACTTACGCAATACAt                                       <  2:636960/139‑1 (MQ=255)                                                                                         cgcgGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                       <  2:691818/139‑1 (MQ=255)                                                                                         cgcgGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                       <  2:1125235/139‑1 (MQ=255)                                                                                             ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCa                                                                                        <  1:1212912/86‑1 (MQ=255)                                                                                             ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCa                                                                                        >  2:1212912/1‑86 (MQ=255)                                                                                             ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCga                                                              >  2:764358/1‑112 (MQ=255)                                                                                             ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCga                                                              <  1:764358/112‑1 (MQ=255)                                                                                             ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGt                                   <  1:795490/139‑1 (MQ=255)                                                                                               cTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTAc                                 >  2:394117/1‑139 (MQ=255)                                                                                                    gCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAAc                            >  1:125630/1‑139 (MQ=255)                                                                                                        tCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCaa                                                                                       >  2:1155001/1‑76 (MQ=255)                                                                                                        tCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCaa                                                                                       <  1:1155001/76‑1 (MQ=255)                                                                                                            gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                                >  1:1028128/1‑79 (MQ=255)                                                                                                            gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                                <  2:1028128/79‑1 (MQ=255)                                                                                                            gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCa                    >  2:62898/1‑139 (MQ=255)                                                                                                                   taCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCg             <  1:64096/139‑1 (MQ=255)                                                                                                                    aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGc            <  1:153953/139‑1 (MQ=255)                                                                                                                       aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                         <  1:782072/43‑1 (MQ=255)                                                                                                                       aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                         >  2:782072/1‑43 (MQ=255)                                                                                                                       aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                         >  2:86887/1‑43 (MQ=255)                                                                                                                       aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                         <  1:86887/43‑1 (MQ=255)                                                                                                                         gCAGCGATTCGCCTGCTGAATACGGTAGCCCATTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTCTAGTACGATCCGAGAAACGGAACTTACGCAATTCATCAGTACGGAACATCCCGCAGGCACGGCTGaaa       >  1:122487/1‑137 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                                <  2:1130051/65‑1 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                                >  1:1130051/1‑65 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacc                                                                                                          >  1:811577/1‑38 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacc                                                                                                          <  2:811577/38‑1 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                      <  1:837284/74‑1 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                      >  2:837284/1‑74 (MQ=255)                                                                                                                              gATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATaaa  <  1:287441/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                ACCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAA  >  NC_000913/383123‑383386

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