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%	119.0 / NA	38	*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 (20/18);  total (20/18)

TTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGT  >  NC_000913/1664990‑1665249                                                                                                                                   |                                                                                                                                  ttCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTGCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCg                                                                                                                           >  1:777280/1‑139 (MQ=255)   cTTGTTATTATTATCGAAAAAGCAGATCACCCTTTCTCCCGTGCACTGGGTGAACTGGCTCGCTAAATGCTGGCACAAAGGCAATCACAACTGTTAATTCCTGTCGTTATTAAAACGCTGTTACGATAAGAGCACCGTa                                                                                                                         <  1:307000/139‑1 (MQ=255)    ttGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTaa                                                                                                                        <  2:250279/139‑1 (MQ=255)            tttATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCCATCACAACTGTTAATTCCTGTCGCGGTTCAACCGCTGCTTCGATAAGATCACCGTCATATTGCCg                                                                                                                >  1:764325/1‑139 (MQ=255)            tttATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCg                                                                                                                >  1:537817/1‑139 (MQ=255)                       aGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGa                                                                                                     <  1:321237/139‑1 (MQ=255)                                         ccGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaa                                                                                   <  1:136261/139‑1 (MQ=255)                                           gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                 >  1:1010792/1‑139 (MQ=255)                                           gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                 <  2:869092/139‑1 (MQ=255)                                           gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                 >  2:246444/1‑139 (MQ=255)                                           gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                 >  1:869092/1‑139 (MQ=255)                                           gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                 <  1:638255/139‑1 (MQ=255)                                           gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                 >  1:343418/1‑139 (MQ=255)                                           gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                                 <  1:32986/139‑1 (MQ=255)                                                cTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAgg                                                                            <  1:348596/139‑1 (MQ=255)                                                 tGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGt                                                                           >  2:909712/1‑139 (MQ=255)                                                   ggTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCgg                                                                                                               >  2:1132676/1‑101 (MQ=255)                                                   ggTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCgg                                                                                                               <  1:1132676/101‑1 (MQ=255)                                                             cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGAt                                                                                                                                >  1:124879/1‑74 (MQ=255)                                                             cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGAt                                                                                                                                <  2:124879/74‑1 (MQ=255)                                                             cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACg                                                               >  1:1007315/1‑139 (MQ=255)                                                             cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACg                                                               >  2:458413/1‑139 (MQ=255)                                                               cGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGtt                                                             >  2:20916/1‑139 (MQ=255)                                                                    aaCGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttt                                                        <  1:246444/139‑1 (MQ=255)                                                                    aaCGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttt                                                        <  1:20916/139‑1 (MQ=255)                                                                                           aaCTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCtt                                 >  1:21212/1‑139 (MQ=255)                                                                                                tttgttCCTGTCCGGGTTAAACCGCTGTTTCGGTAAGATCACCGCAATATGTTCGGGCGGTGACGATTACGCGCAATCCGGCAATAATAGGCTACAGCGTAACGTTATTTTATCTTTTAAAGCAGGCACTGCTTATGCg                            >  1:12157/5‑139 (MQ=255)                                                                                                       cTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                     >  2:5138/1‑139 (MQ=255)                                                                                                         gTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCgg                                                                                                               >  2:868417/1‑47 (MQ=255)                                                                                                         gTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCgg                                                                                                               <  1:868417/47‑1 (MQ=255)                                                                                                         gTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                           >  1:556879/1‑131 (MQ=255)                                                                                                         gTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                           <  2:556879/131‑1 (MQ=255)                                                                                                          tCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCgg                                                                                      <  2:3127/71‑1 (MQ=255)                                                                                                          tCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCgg                                                                                      >  1:3127/1‑71 (MQ=255)                                                                                                               gTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCAGGCACTGCTTTTGCGGCTGGCCTCTTTTGc             >  2:823795/1‑139 (MQ=255)                                                                                                                   aaCCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa         <  2:21212/139‑1 (MQ=255)                                                                                                                   aaCCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa         <  2:777280/139‑1 (MQ=255)                                                                                                                          gTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGt  <  2:1007315/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                  TTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGT  >  NC_000913/1664990‑1665249

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