Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,923,070 G→A A204V (GCA→GTA)  rsmG ← 16S rRNA m(7)G527 methyltransferase, SAM‑dependent, glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,923,0700GA100.0% 81.1 / NA 26A204V (GCA→GTA) rsmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent, glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/13);  total (13/13)

TATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCA  >  NC_000913/3922939‑3923190
                                                                                                                                   |                                                                                                                        
tATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAAt                                                                                                                   <  2:386846/139‑1 (MQ=255)
          aTGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGa                                                                                                         >  1:533725/1‑139 (MQ=255)
            gCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATg                                                                                                       >  2:744287/1‑139 (MQ=255)
                            ggCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAgg                                                                                       >  1:674295/1‑139 (MQ=255)
                                              ttATATCAGTAGAAAGCCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCg                                                                                                <  2:237993/112‑1 (MQ=255)
                                              ttATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCACTTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTATTTTAATCACCACCAGATGACGTTCg                                                                                                >  1:237993/1‑112 (MQ=255)
                                                      gTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGTCGTTCGCCATCAAGGGCTGGAACCTGAAGTTTAACCACTGa                                                             >  1:82228/1‑139 (MQ=255)
                                                       tAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGAt                                                            >  2:793996/1‑139 (MQ=255)
                                                       tAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTCTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGa                                                                                                         >  1:249013/1‑94 (MQ=255)
                                                       tAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTCTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGa                                                                                                         <  2:249013/94‑1 (MQ=255)
                                                        aGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCt                                                                             <  1:1030915/121‑1 (MQ=255)
                                                        aGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCt                                                                             >  2:1030915/1‑121 (MQ=255)
                                                            aaCCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGAc                                                       >  2:851653/1‑139 (MQ=255)
                                                            aaCCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGAc                                                       <  2:370644/139‑1 (MQ=255)
                                                                  gttgttgttAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATcc                                                                                          >  1:444019/1‑98 (MQ=255)
                                                                  gttgttgttAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATcc                                                                                          <  2:444019/98‑1 (MQ=255)
                                                                           aaCAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCggg                                        <  1:793996/139‑1 (MQ=255)
                                                                                  tAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATAttcttc                                           >  2:484828/1‑129 (MQ=255)
                                                                                  tAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATAttcttc                                           <  1:115916/129‑1 (MQ=255)
                                                                                  tAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATAttcttc                                           >  2:115916/1‑129 (MQ=255)
                                                                                  tAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATAttcttc                                           <  1:484828/129‑1 (MQ=255)
                                                                                  tAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACaaa                                 <  2:533725/139‑1 (MQ=255)
                                                                                         tCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCggg                                        <  1:887346/125‑1 (MQ=255)
                                                                                         tCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCggg                                        >  2:887346/1‑125 (MQ=255)
                                                                                                     gATTTTTTTGATAAAAATTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCa              <  1:851653/139‑1 (MQ=255)
                                                                                                                 aaaaaTTAAATTTTATTTACTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCa  <  2:82228/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                        
TATCTTGTTGATGCAAAAGAGTGAACGTGGCGTTAAATGTAACCAGTTATATCAGTAGAAAACCTGGTTGTTGTTAACAGTCTAACCGGTCAATTTTTTATGATTTTTTTGATAAAAATTAAATTTTATTTGCTTTAATCACCACCAGATGACGTTCGCCATCCAGGGCTGGAACCTGAAGTTTAACCACTGATTCGACCTGATATTCTTCGGGCAACAAAGCGATTTCATCTTCCGGCATTTGCCCTTTCA  >  NC_000913/3922939‑3923190

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: