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%	200.2 / 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)

GCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAACGGTGCGC  >  NC_000913/383132‑383394                                                                                                                                |                                                                                                                                        gCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAg                                                                                                                              <  1:437213/139‑1 (MQ=255)            gTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccac                                                                                                                   >  1:917028/1‑139 (MQ=255)              cTTTTATGCAGGCGCATACAGACCGCGTATCCCACAAGCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                                 <  2:1106959/139‑1 (MQ=255)                       aGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAgg                                                                                                        >  2:644833/1‑139 (MQ=255)                         gcgcATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTg                                                                                                      <  2:879920/139‑1 (MQ=255)                                   aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGaa                                                                                                                                      >  2:832580/1‑97 (MQ=255)                                   aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGaa                                                                                                                                      <  1:832580/97‑1 (MQ=255)                                   aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                                 <  1:720818/118‑1 (MQ=255)                                   aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                                 >  2:720818/1‑118 (MQ=255)                                      gcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                         <  1:450933/139‑1 (MQ=255)                                      gcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                         <  1:439710/139‑1 (MQ=255)                                       cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAACGTCCAGGTGAGATCAAGGTCTTc                                                                                        <  2:918908/139‑1 (MQ=255)                                                   tCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACg                                                                            >  2:760837/1‑139 (MQ=255)                                                      aaCGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATc                                                                         >  1:202799/1‑139 (MQ=255)                                                             ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGaaa                                                                  >  1:56199/1‑139 (MQ=255)                                                             ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGaaa                                                                  >  2:1225670/1‑139 (MQ=255)                                                                      ccACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                                         <  2:917028/139‑1 (MQ=255)                                                                      ccACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                                         <  2:1466364/139‑1 (MQ=255)                                                                        aCTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                                       <  1:1225670/139‑1 (MQ=255)                                                                          tcctcACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                                 >  1:1542603/1‑79 (MQ=255)                                                                          tcctcACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                                 <  2:1542603/79‑1 (MQ=255)                                                                                cgcgGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                               <  2:202799/139‑1 (MQ=255)                                                                                 gcgGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAg                                                                                >  1:316701/1‑105 (MQ=255)                                                                                 gcgGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAg                                                                                <  2:316701/105‑1 (MQ=255)                                                                                 gcgGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATc                                              >  1:1577143/1‑139 (MQ=255)                                                                                    ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                            >  1:836169/1‑74 (MQ=255)                                                                                    ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                            <  2:836169/74‑1 (MQ=255)                                                                                               tCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATcc                                >  2:324289/1‑139 (MQ=255)                                                                                                 cTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCg                              >  2:653625/1‑139 (MQ=255)                                                                                                       gTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCa                        <  1:653625/139‑1 (MQ=255)                                                                                                       gTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCa                        <  2:56199/139‑1 (MQ=255)                                                                                                       gTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCa                        <  1:644833/139‑1 (MQ=255)                                                                                                        tataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                                 <  2:1572073/49‑1 (MQ=255)                                                                                                        tataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                                 >  1:1572073/1‑49 (MQ=255)                                                                                                         ataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCAcc                      >  1:268860/1‑139 (MQ=255)                                                                                                         ataCGATGCAGCAATTCGCCTGCTGAATACGGTAGCCGTTTGTCACCAATGCCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCAcc                      >  1:1251219/1‑139 (MQ=255)                                                                                                           aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGc                    <  2:466522/139‑1 (MQ=255)                                                                                                           aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCG‑TTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCgag                                                                     <  2:1134138/89‑1 (MQ=255)                                                                                                           aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCG‑TTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCgag                                                                     >  1:1134138/1‑89 (MQ=255)                                                                                                              aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATAc                                                 >  2:23196/1‑107 (MQ=255)                                                                                                              aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATAc                                                 <  1:23196/107‑1 (MQ=255)                                                                                                               tGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                                       <  1:321671/100‑1 (MQ=255)                                                                                                               tGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                                       >  2:321671/1‑100 (MQ=255)                                                                                                                gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                                 >  1:541133/1‑41 (MQ=255)                                                                                                                gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                                 <  2:541133/41‑1 (MQ=255)                                                                                                                gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                         <  2:944682/65‑1 (MQ=255)                                                                                                                gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                         >  1:944682/1‑65 (MQ=255)                                                                                                                 cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAAc                                                            >  1:851062/1‑93 (MQ=255)                                                                                                                 cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAAc                                                            <  2:851062/93‑1 (MQ=255)                                                                                                                 cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGataa              >  2:819610/1‑139 (MQ=255)                                                                                                                  aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                              >  1:1077337/1‑74 (MQ=255)                                                                                                                  aGCGATTCGCCTGCTGAATACGGTAGCCGTATGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                              <  2:1077337/74‑1 (MQ=255)                                                                                                                      aTTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATccc                               >  2:1299131/1‑117 (MQ=255)                                                                                                                      aTTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATccc                               <  1:1299131/117‑1 (MQ=255)                                                                                                                       ttCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                            >  2:25375/1‑39 (MQ=255)                                                                                                                       ttCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                            <  1:25375/39‑1 (MQ=255)                                                                                                                         cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                                    >  2:461855/1‑61 (MQ=255)                                                                                                                         cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                                    <  1:461855/61‑1 (MQ=255)                                                                                                                         cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAAc                                                                 >  2:133170/1‑80 (MQ=255)                                                                                                                         cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAAc                                                                 <  1:133170/80‑1 (MQ=255)                                                                                                                         cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTAGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAACGGt      >  1:1065647/1‑139 (MQ=255)                                                                                                                             tGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGATATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACTGAACATCCCGCAGGCACCGCTGATAATAAACGGTgcgc  <  1:324289/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                        GCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAACGGTGCGC  >  NC_000913/383132‑383394

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: ‑