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

CGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGT  >  NC_000913/4024208‑4024467                                                                                                                                 |                                                                                                                                    cGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGa                                                                                                                           >  1:201929/1‑139 (MQ=255)      cAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCg                                                                                                                      <  2:1536701/139‑1 (MQ=255)        aCGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                    <  1:1208835/139‑1 (MQ=255)        aCGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                    <  2:201929/139‑1 (MQ=255)         cGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTa                                                                                                                   >  2:198229/1‑139 (MQ=255)                  ccATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                          <  2:238132/139‑1 (MQ=255)                   cATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                         <  2:1349315/139‑1 (MQ=255)                           ccATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                    <  1:317872/120‑1 (MQ=255)                           ccATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                    >  2:317872/1‑120 (MQ=255)                              tatTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAAtc                                                                                              >  1:227476/1‑139 (MQ=255)                                 tttGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATctct                                                                                           <  2:490117/139‑1 (MQ=255)                                   tGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTa                                                                                                 >  2:121680/1‑131 (MQ=255)                                   tGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTa                                                                                                 <  1:121680/131‑1 (MQ=255)                                        cGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTAttct                                                                                                               >  2:354132/1‑112 (MQ=255)                                        cGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTAttct                                                                                                               <  1:354132/112‑1 (MQ=255)                                        cGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATtaa                                                                                    <  2:227476/139‑1 (MQ=255)                                                      gCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAAt                                                                      >  2:487767/1‑139 (MQ=255)                                                       cGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcac                                                                                                           >  1:731422/1‑101 (MQ=255)                                                       cGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcac                                                                                                           <  2:731422/101‑1 (MQ=255)                                                                gaTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                          >  2:1074808/1‑93 (MQ=255)                                                                gaTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagc                                                            >  1:320765/1‑139 (MQ=255)                                                                gaTGCCGCATGGCTGGCAGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTtaca                                                                                                          <  1:1074808/93‑1 (MQ=255)                                                                   gCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTa                                                                            <  2:884950/120‑1 (MQ=255)                                                                   gCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTa                                                                            >  1:884950/1‑120 (MQ=255)                                                                      gCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGt                                                                                                                                  <  1:898733/63‑1 (MQ=255)                                                                      gCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGt                                                                                                                                  >  2:898733/1‑63 (MQ=255)                                                                       cATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTc                                                     <  2:320765/139‑1 (MQ=255)                                                                        gtggctggcTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                         >  1:1595686/2‑86 (MQ=255)                                                                        aTGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                         <  2:1595686/86‑1 (MQ=255)                                                                          ggctggctGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                          >  1:1667983/1‑83 (MQ=255)                                                                          ggctggctGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                          <  2:1667983/83‑1 (MQ=255)                                                                            ctggctgCCACCACGGATGCTAATGTCAAAGCACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCttt                                                                                         >  1:1113082/1‑98 (MQ=255)                                                                            ctggctgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCTCACTGTGCTTAATCTCttt                                                                                         <  2:1113082/98‑1 (MQ=255)                                                                            ctggctgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCAt                                                <  1:198229/139‑1 (MQ=255)                                                                                 tgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGt                                           <  1:487767/139‑1 (MQ=255)                                                                                         cGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACt                                                                                                        >  2:717015/1‑70 (MQ=255)                                                                                         cGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACt                                                                                                        <  1:717015/70‑1 (MQ=255)                                                                                          ggATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGaaa                                  <  1:1555206/139‑1 (MQ=255)                                                                                             tGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                    >  1:942001/1‑54 (MQ=255)                                                                                             tGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                    <  2:942001/54‑1 (MQ=255)                                                                                             tGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTa                                                                                                 <  2:230999/73‑1 (MQ=255)                                                                                             tGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTa                                                                                                 >  1:230999/1‑73 (MQ=255)                                                                                               cTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                    >  1:148746/1‑52 (MQ=255)                                                                                               cTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                    <  2:148746/52‑1 (MQ=255)                                                                                                   tGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGa                         <  2:1665445/139‑1 (MQ=255)                                                                                                    gTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                          >  1:1501222/1‑57 (MQ=255)                                                                                                    gTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTcaca                                                                                                          <  2:1501222/57‑1 (MQ=255)                                                                                                      cAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCaa                                                                       <  1:270053/90‑1 (MQ=255)                                                                                                      cAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCaa                                                                       >  2:270053/1‑90 (MQ=255)                                                                                                       aaaaCACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGc                     >  2:1295415/1‑139 (MQ=255)                                                                                                        aaaCACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagc                                                            <  2:861564/99‑1 (MQ=255)                                                                                                        aaaCACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagc                                                            >  1:861564/1‑99 (MQ=255)                                                                                                            acTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTa                                                                                                                   <  1:474263/40‑1 (MQ=255)                                                                                                            acTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTa                                                                                                                   >  2:474263/1‑40 (MQ=255)                                                                                                            acTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCg                <  1:1295415/139‑1 (MQ=255)                                                                                                               gTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATtaat                                                                                   >  1:301121/1‑69 (MQ=255)                                                                                                               gTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATtaat                                                                                   <  2:301121/69‑1 (MQ=255)                                                                                                                      aTTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCtt      >  1:1557861/1‑139 (MQ=255)                                                                                                                        tGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                         >  2:1542807/1‑38 (MQ=255)                                                                                                                        tGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                         <  1:1542807/38‑1 (MQ=255)                                                                                                                        tGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgag                                                             <  1:1333260/82‑1 (MQ=255)                                                                                                                        tGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgag                                                             >  2:1333260/1‑82 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTtta                                                                                        >  1:92905/1‑53 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTtta                                                                                        <  2:92905/53‑1 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATtaa                                                                                    >  1:1374115/1‑57 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATtaa                                                                                    <  2:1374115/57‑1 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAGTAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGt  >  1:1105460/1‑139 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgag                                                             <  1:998144/80‑1 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgag                                                             >  2:998144/1‑80 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGcc                                                       >  2:1599324/1‑86 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGcc                                                       >  2:1593312/1‑86 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGcc                                                       <  1:1593312/86‑1 (MQ=255)                                                                                                                          cGTTTTAAAGTTTGCGGAACTCCGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGcc                                                       <  1:1599324/86‑1 (MQ=255)                                                                                                                              ttAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGcc        <  2:1582736/129‑1 (MQ=255)                                                                                                                                 |                                                                                                                                    CGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGGT  >  NC_000913/4024208‑4024467

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