Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,665,120	G→A	*205* (TAG→TAA)	dmsD →	twin‑argninine leader‑binding protein for DmsA and TorA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,665,120	0	G	A	100.0%	129.1 / NA	40	*205* (TAG→TAA)	dmsD	twin‑argninine leader‑binding protein for DmsA and TorA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/18);  total (22/18)

ACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCG  >  NC_000913/1664985‑1665255                                                                                                                                        |                                                                                                                                        acGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGAt                                                                                                                                      <  1:1381789/139‑1 (MQ=255)           gATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAata                                                                                                                            <  1:1084137/139‑1 (MQ=255)              gTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTg                                                                                                                         >  1:379845/1‑139 (MQ=255)                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                       >  1:931559/1‑139 (MQ=255)                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                       <  2:379845/139‑1 (MQ=255)                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                       <  1:915399/139‑1 (MQ=255)                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaa                                                                                         >  1:1460859/1‑137 (MQ=255)                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaa                                                                                         <  2:1460859/137‑1 (MQ=255)                                                         gTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTAc                                                                              <  2:931559/139‑1 (MQ=255)                                                                  cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCgg                                                                                            <  1:148799/116‑1 (MQ=255)                                                                  cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCgg                                                                                            >  2:148799/1‑116 (MQ=255)                                                                  cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaataa                                                                                      <  2:923996/122‑1 (MQ=255)                                                                  cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaa                                                                                         <  2:475911/119‑1 (MQ=255)                                                                  cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaa                                                                                         >  1:475911/1‑119 (MQ=255)                                                                  cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAACATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaataa                                                                                      >  1:923996/1‑122 (MQ=255)                                                                           cGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATcc                                                                                              >  1:251540/1‑105 (MQ=255)                                                                           cGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATcc                                                                                              <  2:251540/105‑1 (MQ=255)                                                                                 cgcAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTa                                                      >  1:1140684/1‑139 (MQ=255)                                                                                           aTCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTACCTCTTAAAgcacgcac                                            >  2:1431291/1‑139 (MQ=255)                                                                                                aaCTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaa                                                                                         >  2:464386/1‑89 (MQ=255)                                                                                                aaCTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaa                                                                                         <  1:464386/89‑1 (MQ=255)                                                                                                  cTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTAt                                                            <  1:1044513/116‑1 (MQ=255)                                                                                                  cTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTAt                                                            >  2:1044513/1‑116 (MQ=255)                                                                                                       aaTTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                >  2:850348/1‑139 (MQ=255)                                                                                                          tCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                       <  1:228277/81‑1 (MQ=255)                                                                                                          tCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                       >  2:228277/1‑81 (MQ=255)                                                                                                            cTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTc                                                                        >  1:631401/1‑94 (MQ=255)                                                                                                            cTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTc                                                                        <  2:631401/94‑1 (MQ=255)                                                                                                             tGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCgg                                                                                                                     <  2:1203186/48‑1 (MQ=255)                                                                                                             tGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCgg                                                                                                                     >  1:1203186/1‑48 (MQ=255)                                                                                                               tCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTaaa                                                    >  2:1497672/1‑111 (MQ=255)                                                                                                               tCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTaaa                                                    <  1:1497672/111‑1 (MQ=255)                                                                                                               tCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCttt                                      >  1:19921/1‑125 (MQ=255)                                                                                                               tCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCttt                                      <  2:19921/125‑1 (MQ=255)                                                                                                               tCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                        >  2:358571/1‑139 (MQ=255)                                                                                                               tCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                        >  1:1343896/1‑139 (MQ=255)                                                                                                                               gTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAg                                                   <  2:654558/96‑1 (MQ=255)                                                                                                                               gTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAg                                                   >  1:654558/1‑96 (MQ=255)                                                                                                                               gTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGt        >  1:464679/1‑139 (MQ=255)                                                                                                                                     aTAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCg  >  2:20321/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        ACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCG  >  NC_000913/1664985‑1665255

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