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%	123.7 / NA	40	*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 (16/24);  total (16/24)

ATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGG  >  NC_000913/4024204‑4024466                                                                                                                                     |                                                                                                                                   atcCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTg                                                                                                                              <  2:389461/139‑1 (MQ=255) atcCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTg                                                                                                                              <  1:784406/139‑1 (MQ=255)        cgcAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACt                                                                                                                       <  2:408131/139‑1 (MQ=255)        cgcAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACt                                                                                                                       <  1:913785/139‑1 (MQ=255)          cAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCg                                                                                                                     <  1:422001/139‑1 (MQ=255)            aCGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                   <  1:139949/139‑1 (MQ=255)            aCGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                   <  2:280750/139‑1 (MQ=255)                       cATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                        <  2:319155/139‑1 (MQ=255)                       cATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                        <  2:288331/139‑1 (MQ=255)                         tCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                   >  2:893675/1‑126 (MQ=255)                         tCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                   <  1:893675/126‑1 (MQ=255)                                cATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                        <  2:958250/130‑1 (MQ=255)                                cATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                        >  1:958250/1‑130 (MQ=255)                                            cGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTa                                                                                                >  1:560086/1‑126 (MQ=255)                                            cGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTa                                                                                                <  2:560086/126‑1 (MQ=255)                                                      aCTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAg                                                                         >  2:785678/1‑139 (MQ=255)                                                        tGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCa                                                                       >  2:304719/1‑139 (MQ=255)                                                                  aagaTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGga                                                             <  1:159289/139‑1 (MQ=255)                                                                     aTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                        >  2:508826/1‑93 (MQ=255)                                                                     aTGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                        <  1:508826/93‑1 (MQ=255)                                                                               gctggctgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                        <  2:941256/83‑1 (MQ=255)                                                                               gctggctgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCacac                                                                                                        >  1:941256/1‑83 (MQ=255)                                                                               gctggctgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATctct                                                                                          <  2:316004/97‑1 (MQ=255)                                                                               gctggctgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATctct                                                                                          >  1:316004/1‑97 (MQ=255)                                                                                ctggctgCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCAt                                               >  2:766899/1‑139 (MQ=255)                                                                                         accacGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCgg                                      >  1:435070/1‑139 (MQ=255)                                                                                          ccacGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCttt                                                                                        <  1:121634/88‑1 (MQ=255)                                                                                          ccacGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCttt                                                                                        >  2:121634/1‑88 (MQ=255)                                                                                                 tGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATg                              <  2:435070/139‑1 (MQ=255)                                                                                                   cTAATGTCAAAACACTGTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGc                            >  1:304347/1‑139 (MQ=255)                                                                                                                   gTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                   >  1:178782/1‑36 (MQ=255)                                                                                                                   gTTTGGGATTGCGTTTTAAAGTTTGCGGAACTCGGt                                                                                                                   <  2:178782/36‑1 (MQ=255)                                                                                                                          aTTGCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCtt     <  1:766899/139‑1 (MQ=255)                                                                                                                             gCGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCgg  <  1:556555/139‑1 (MQ=255)                                                                                                                              cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTa                                                                                                >  2:710453/1‑44 (MQ=255)                                                                                                                              cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTa                                                                                                <  1:710453/44‑1 (MQ=255)                                                                                                                              cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATtaa                                                                                   >  2:643924/1‑57 (MQ=255)                                                                                                                              cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATtaa                                                                                   <  1:643924/57‑1 (MQ=255)                                                                                                                              cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagcg                                                          <  2:796284/82‑1 (MQ=255)                                                                                                                              cGTTTTAAAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGgagcg                                                          >  1:796284/1‑82 (MQ=255)                                                                                                                                     |                                                                                                                                   ATCCCGGCGCAACGAGCCAGCCCATCTGCCCCATATTTTGCAACGTATTGCGCACTGGCGTGGAGAAGATGCCGCATGGCTGGCTGCCACCACGGATGCTAATGTCAAAACACTGTTTGGGATTGCGTTTTAGAGTTTGCGGAACTCGGTATTCTTCACACTGTGCTTAATCTCTTTATTAATAAGATTAAGCAATAGCATGGAGCGAGCCTCACCATCGGGTTCGGTGAAAATGGCCTGAAAGCCTTCGAACGCGCCTTCGG  >  NC_000913/4024204‑4024466

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