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%	197.5 / NA	61	R226Q (CGG→CAG)	yaiP	putative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/30);  total (31/30)

ACCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAACGGTG  >  NC_000913/383123‑383391                                                                                                                                         |                                                                                                                                     aCCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGAGTCGCCTGCTGa                                                                                                                                    >  1:872237/1‑139 (MQ=255)           cTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGc                                                                                                                          >  2:1038436/1‑139 (MQ=255)                       cTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGGTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTCTGccacca                                                                                                              <  2:126501/139‑1 (MQ=255)                         tttATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAAt                                                                                                            >  2:1282035/1‑139 (MQ=255)                            aTGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                         <  1:1038436/139‑1 (MQ=255)                                aGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCTGCGATTCGCCTGCTGAATATGGTAGCCGTTTGCCACCAATGTCCAgg                                                                                                     <  1:857720/139‑1 (MQ=255)                                 ggCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGt                                                                                                    <  2:480507/139‑1 (MQ=255)                                  gcgcATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTg                                                                                                   >  2:828849/1‑139 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGaa                                                                                                                                   <  2:591122/97‑1 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGaa                                                                                                                                   >  1:591122/1‑97 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccac                                                                                                                <  1:1022857/116‑1 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccac                                                                                                                >  2:1022857/1‑116 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAgg                                                                                                     >  1:1092500/1‑127 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAgg                                                                                                     <  2:1092500/127‑1 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGt                                                                                         >  2:332617/1‑139 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGt                                                                                         >  1:1074277/1‑139 (MQ=255)                                            aCCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGt                                                                                         >  1:623173/1‑139 (MQ=255)                                               gcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCTCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                      <  2:1074277/139‑1 (MQ=255)                                                cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                                     <  2:96039/139‑1 (MQ=255)                                                  tATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCATGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGa                                                                                   >  2:993847/1‑139 (MQ=255)                                                          aaTCCAACGCCGCCAGCGACGCCACTCCTCATGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                           <  1:993847/139‑1 (MQ=255)                                                          aaTCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTa                                                                           <  1:332617/139‑1 (MQ=255)                                                             ccAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGa                                                                        <  2:872237/139‑1 (MQ=255)                                                               aaCGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATc                                                                      >  2:1237248/1‑139 (MQ=255)                                                                              gCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACtt                                                       >  2:181323/1‑139 (MQ=255)                                                                               ccACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                                      <  1:1237248/139‑1 (MQ=255)                                                                                  ctcctcACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGc                                                   >  1:542766/1‑139 (MQ=255)                                                                                    cctcACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                                 <  2:623173/139‑1 (MQ=255)                                                                                             ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGt                                        >  1:1085763/1‑139 (MQ=255)                                                                                               cTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                                      <  2:1238939/123‑1 (MQ=255)                                                                                               cTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                                      >  1:1238939/1‑123 (MQ=255)                                                                                                   tGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGaa                                  <  2:342447/139‑1 (MQ=255)                                                                                                        tCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATcc                             >  2:515316/1‑139 (MQ=255)                                                                                                            gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACg                                                                         >  1:738485/1‑91 (MQ=255)                                                                                                            gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACg                                                                         <  2:738485/91‑1 (MQ=255)                                                                                                            gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCa                         >  2:687511/1‑139 (MQ=255)                                                                                                               ggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                                 <  2:1254838/112‑1 (MQ=255)                                                                                                               ggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                                 >  1:1238701/1‑112 (MQ=255)                                                                                                               ggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                                 <  2:1238701/112‑1 (MQ=255)                                                                                                               ggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                                 >  1:1254838/1‑112 (MQ=255)                                                                                                                   taCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCg                  <  2:1085763/139‑1 (MQ=255)                                                                                                                    aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCg                  >  1:694593/1‑138 (MQ=255)                                                                                                                    aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCg                  <  2:694593/138‑1 (MQ=255)                                                                                                                    aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGc                 <  2:1059986/139‑1 (MQ=255)                                                                                                                    aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGc                 >  1:1224652/1‑139 (MQ=255)                                                                                                                    aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGc                 >  1:1059986/1‑139 (MQ=255)                                                                                                                      gATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTg               <  1:1282035/139‑1 (MQ=255)                                                                                                                       aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                                    >  1:1115211/1‑101 (MQ=255)                                                                                                                       aTGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                                    <  2:1115211/101‑1 (MQ=255)                                                                                                                         gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCg                  <  1:537673/133‑1 (MQ=255)                                                                                                                         gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCg                  >  2:537673/1‑133 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                      >  1:755003/1‑64 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                                      <  2:755003/64‑1 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAAc                                                         >  1:177364/1‑93 (MQ=255)                                                                                                                          cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAAc                                                         <  2:177364/93‑1 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTga                                                                                                  <  1:1010068/51‑1 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTga                                                                                                  >  2:1010068/1‑51 (MQ=255)                                                                                                                           aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGataat          >  1:231785/1‑139 (MQ=255)                                                                                                                             cGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                                    <  2:1011972/95‑1 (MQ=255)                                                                                                                             cGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                                    >  1:1011972/1‑95 (MQ=255)                                                                                                                                   gCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAACGGTg  <  1:687511/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                     ACCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAATAAACGGTG  >  NC_000913/383123‑383391

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