Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	666,217	G→T	L33M (CTG→ATG)	mrdB ←	SEDS family protein MrdB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	666,217	0	G	T	100.0%	104.6 / NA	31	L33M (CTG→ATG)	mrdB	SEDS family protein MrdB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/18);  total (13/18)

GGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAGGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATCCCAGAATGTTTTTTTATTCGGATTATCCGTCATGATTAATGGTCCTCCGC  >  NC_000913/666087‑666330                                                                                                                                   |                                                                                                                  ggggCCCAGCTTTCATAAACGCGTAGAGGAATCTGCGCCATGACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAtgg                                                                                                                 <  2:210298/133‑1 (MQ=255) ggggCCCAGCCTTCATAAACGCGTGGAGGAATTTGTGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAtgg                                                                                                                 >  2:212390/1‑133 (MQ=255) ggggCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAtgg                                                                                                                 >  1:311472/1‑133 (MQ=255) ggggCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCAGGCGCTCCAGATAACCAtgg                                                                                                                 <  1:91257/133‑1 (MQ=255)    gCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGgcgc                                                                                                              <  1:23933/133‑1 (MQ=255)    gCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGgcgc                                                                                                              <  2:143862/133‑1 (MQ=255)    gCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCAAGCCAATATCCTGGCCGGTGGCGCTCCAGATAACCATGgcgc                                                                                                              >  2:263441/1‑133 (MQ=255)     cccAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCt                                                                                                             >  1:482629/1‑133 (MQ=255)     cccAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCt                                                                                                             >  2:448488/1‑133 (MQ=255)          ccTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTaaa                                                                                                        >  2:439689/1‑133 (MQ=255)                 aaaCGCGTGGAGGAATTTGCGCCATCACCCCCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCaa                                                                                                 <  1:215540/133‑1 (MQ=255)                    cgcgTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGc                                                                                              >  2:337209/1‑133 (MQ=255)                      cgTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCa                                                                                            <  2:307302/133‑1 (MQ=255)                           aggaATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAg                                                                                       <  1:34406/133‑1 (MQ=255)                             gaATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGAt                                                                                     <  2:311472/133‑1 (MQ=255)                                tttGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATcag                                                                                  <  1:439689/133‑1 (MQ=255)                                 ttGCGCCAGCACCACCATCATGCCCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATcagc                                                                                 <  1:263441/133‑1 (MQ=255)                                   gcgcCATCGCCACCATGACGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATcagcag                                                                               <  1:290587/133‑1 (MQ=255)                                         tCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATtgt                                                                         <  1:212390/133‑1 (MQ=255)                                             caccaTGACGACCAGTCCGATCACGATTTGACCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGa                                                                     <  1:285857/133‑1 (MQ=255)                                               ccaTGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATc                                                                   <  2:519129/133‑1 (MQ=255)                                                     tgaCCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATg                                                             <  1:534834/133‑1 (MQ=255)                                                             cccATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTAt                                                     >  1:369656/1‑133 (MQ=255)                                                               cATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATcc                                                   <  2:22806/133‑1 (MQ=255)                                                                 tCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGGAGCTTTGAGGGATCGAAATGGTCTTTATCCCa                                                 >  2:237268/1‑133 (MQ=255)                                                                 tCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATCCCa                                                 <  2:91871/133‑1 (MQ=255)                                                                 tCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATCCCa                                                 >  1:530071/1‑133 (MQ=255)                                                                          ggCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGAGTTGATCCCAGAATGtttt                                        >  1:286306/1‑133 (MQ=255)                                                                              gATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATCCCAGAATGTTTTTTTa                                    <  2:369656/133‑1 (MQ=255)                                                                                ttttACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATCCCAGAATGTTTTTTTAtt                                  <  1:400003/133‑1 (MQ=255)                                                                                          catcatGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGTGATGGACTTTATACCAAAATGTTTTATTATTCGGATTATcc                        >  2:240136/1‑133 (MQ=255)                                                                                          catcatGCCAATATCCTGACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATCCCAGAATGTTTTTTTATTCGGATTATcc                        >  1:409163/1‑133 (MQ=255)                                                                                              atGCCAATAACCTGAACGCTGGCGTTCCAGATAACCATGGCCCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATCCCAGAATGTTTTTTTATTCGGATTATCCGTCa                    <  2:286306/133‑1 (MQ=255)                                                                                                           gACCGCTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATCCCAGAATGTTTTTTTATTCGGATTATCCGTCATGATTAATGGtcc       >  2:440363/1‑133 (MQ=255)                                                                                                                cTGGCGCTCCAGATAACCATGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATCCCAGAATGTTTTTTTATTCGGATTATCCGTCATGATTAATGGTCCTCCGc  >  2:124469/1‑133 (MQ=255)                                                                                                                                   |                                                                                                                  GGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGATGACCAGACCCATCGCGATTTGGCCGATTTTACGCTCCATCATGCCAATATCCTGACCGCTGGCGCTCCAGATAACCAGGGCGCTGTAAACCAGCAATGCCAGTAAGATCAGCAGCATTGTGGGATCGAGATGGACTTTATCCCAGAATGTTTTTTTATTCGGATTATCCGTCATGATTAATGGTCCTCCGC  >  NC_000913/666087‑666330

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