Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,024,336	G→A	*261* (TAG→TAA)  L162L (CTC→CTT)	tatD → rfaH ←	quality control of Tat‑exported FeS proteins, Mg‑dependent cytoplasmic DNase transcription antitermination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,024,336	0	G	A	100.0%	241.5 / NA	74	*261* (TAG→TAA)  L162L (CTC→CTT)	tatD rfaH	quality control of Tat‑exported FeS proteins, Mg‑dependent cytoplasmic DNase transcription antitermination protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/40);  total (34/40)

GGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGG  >  NC_000913/4024209‑4024466                                                                                                                                |                                                                                                                                   ggCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGaa                                                                                                                         <  2:517040/139‑1 (MQ=255)       aCGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                   <  1:554296/139‑1 (MQ=255)       aCGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                   <  2:858737/139‑1 (MQ=255)        cGAGCCAGCCCATCTGCCCCATCTTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTATGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTa                                                                                                                  >  1:97413/1‑139 (MQ=255)        cGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTa                                                                                                                  >  1:1314685/1‑139 (MQ=255)        cGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTa                                                                                                                  <  1:1572641/139‑1 (MQ=255)               gccCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTca                                                                                                           <  2:1314685/139‑1 (MQ=255)                 ccATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         <  2:1150404/139‑1 (MQ=255)                 ccATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         <  1:761167/139‑1 (MQ=255)                 ccATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         <  2:509793/139‑1 (MQ=255)                  cATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                        <  1:307104/139‑1 (MQ=255)                     cTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                   >  1:377022/1‑125 (MQ=255)                     cTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                   <  2:377022/125‑1 (MQ=255)                     cTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACtgt                                                                                                     >  2:877347/1‑139 (MQ=255)                          ccATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACt                                                                                                                       >  2:391617/1‑116 (MQ=255)                          ccATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACt                                                                                                                       <  1:391617/116‑1 (MQ=255)                                tttGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATctct                                                                                          <  2:97413/139‑1 (MQ=255)                                  tGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCttt                                                                                        >  2:32724/1‑139 (MQ=255)                                       cGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGTTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATtaa                                                                                   >  1:430522/1‑139 (MQ=255)                                       cGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACGCTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCg                                                                                                                     <  1:234478/105‑1 (MQ=255)                                       cGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACGCTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCg                                                                                                                     >  2:234478/1‑105 (MQ=255)                                       cGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         <  1:1442198/117‑1 (MQ=255)                                       cGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         >  2:1442198/1‑117 (MQ=255)                                       cGTATCGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATtaa                                                                                   <  2:757599/139‑1 (MQ=255)                                              cgcACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCg                                                                                                                            >  1:989541/1‑91 (MQ=255)                                              cgcACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGAtt                                                                            >  1:644617/1‑139 (MQ=255)                                              cgcACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCCAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGct                                                                                                                            <  2:989541/91‑2 (MQ=255)                                                 aCTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAg                                                                         <  2:644617/139‑1 (MQ=255)                                                   tGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACt                                                                                                                       >  2:340810/1‑91 (MQ=255)                                                   tGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACt                                                                                                                       <  1:340810/91‑1 (MQ=255)                                                     gCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTAt                                                                                                                 <  1:1265827/95‑1 (MQ=255)                                                     gCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTAt                                                                                                                 >  2:1265827/1‑95 (MQ=255)                                                               gaTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         <  2:1600651/93‑1 (MQ=255)                                                               gaTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         >  1:1600651/1‑93 (MQ=255)                                                               gaTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTa                                                                           <  1:1149738/123‑1 (MQ=255)                                                               gaTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTa                                                                           >  2:1149738/1‑123 (MQ=255)                                                                 tGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagcga                                                         >  1:660026/1‑139 (MQ=255)                                                                       aTGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACt                                                                                                                       <  2:976061/71‑1 (MQ=255)                                                                       aTGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACt                                                                                                                       >  1:976061/1‑71 (MQ=255)                                                                                   caccacGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTAttcttc                                                                                                            >  2:1087352/1‑70 (MQ=255)                                                                                   caccacGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTAttcttc                                                                                                            <  1:1087352/70‑1 (MQ=255)                                                                                    accacGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCg                                                                                                                     >  2:1399576/1‑60 (MQ=255)                                                                                    accacGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCg                                                                                                                     <  1:1399576/60‑1 (MQ=255)                                                                                        cGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACtg                                                                                                      >  2:1502581/1‑71 (MQ=255)                                                                                        cGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACtg                                                                                                      <  1:1502581/71‑1 (MQ=255)                                                                                                  tGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGa                        <  2:430522/139‑1 (MQ=255)                                                                                                  tGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGa                        <  1:877347/139‑1 (MQ=255)                                                                                                       aaaCACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCa                                                                       >  1:1513914/1‑87 (MQ=255)                                                                                                       aaaCACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCa                                                                       <  2:1513914/87‑1 (MQ=255)                                                                                                       aaaCACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCgg                                      <  2:41087/120‑1 (MQ=255)                                                                                                       aaaCACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCgg                                      >  1:41087/1‑120 (MQ=255)                                                                                                              gTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTa                                                                                                <  2:1218678/55‑1 (MQ=255)                                                                                                              gTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTa                                                                                                >  1:1218678/1‑55 (MQ=255)                                                                                                              gTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAGTAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAAc            >  2:119219/1‑139 (MQ=255)                                                                                                              gTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAAc            >  2:1596318/1‑139 (MQ=255)                                                                                                               tttGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACt                                                                                                       <  2:1057755/47‑1 (MQ=255)                                                                                                               tttGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACt                                                                                                       >  1:1057755/1‑47 (MQ=255)                                                                                                                 tGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAAcgcg         <  2:660026/139‑1 (MQ=255)                                                                                                                    gATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         >  1:1200832/1‑40 (MQ=255)                                                                                                                    gATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         <  2:1200832/40‑1 (MQ=255)                                                                                                                       tGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCg   >  1:1189700/1‑139 (MQ=255)                                                                                                                        gCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAGTAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCgg  <  1:119219/139‑1 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         <  2:575907/35‑1 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                         >  1:575907/1‑35 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAg                                                                         >  1:885616/1‑67 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAg                                                                         <  2:885616/67‑1 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCa                                                                       <  1:28039/69‑1 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCa                                                                       >  2:28039/1‑69 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgag                                                            >  1:1317192/1‑80 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgag                                                            <  2:1317192/80‑1 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGcc       >  1:939397/1‑133 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGcc       <  2:939397/133‑1 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCgg  <  2:417627/138‑1 (MQ=255)                                                                                                                         cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCgg  >  1:417627/1‑138 (MQ=255)                                                                                                                                |                                                                                                                                   GGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGG  >  NC_000913/4024209‑4024466

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