Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,091,787 A→G G218G (GGA→GGG)  rsmE → 16S rRNA m(3)U1498 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,091,7870AG92.1% 96.9 / ‑5.3 38G218G (GGA→GGG) rsmE16S rRNA m(3)U1498 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (17/18);  minor base C (0/3);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

ATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCATCGCAAACATCAACATCAAGAA  >  NC_000913/3091642‑3091930
                                                                                                                                                 |                                                                                                                                               
aTCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTcg                                                                                                                                            >  2:55065/1‑151 (MQ=255)
    tCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTcgc                                                                                                                                           <  2:542678/148‑1 (MQ=255)
                      cAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACt                                                                                                                      <  1:98604/151‑1 (MQ=255)
                       aGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTg                                                                                                                     <  1:126061/151‑1 (MQ=255)
                          aTCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTgcgc                                                                                                                  <  1:308349/151‑1 (MQ=255)
                            cAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCt                                                                                                                 <  2:661112/150‑1 (MQ=255)
                            cAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCt                                                                                                                 <  1:216067/150‑1 (MQ=255)
                            cAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTc                                                                                                                >  2:516701/1‑151 (MQ=255)
                            cAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTc                                                                                                                >  1:581015/1‑151 (MQ=255)
                             aaTACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTc                                                                                                                >  2:196579/1‑150 (MQ=255)
                              aTACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCAc                                                                                                              >  1:373337/1‑151 (MQ=255)
                                    tGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCAt                                                                                                        >  2:176974/1‑151 (MQ=255)
                                            ccGGTTGAACGCGCCCGCCTGCTGTTTGCCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGCCTGCCCGCTATAAATTTACTGATACCCTGTTGGGCCCTCGCTTTTTGCTTACAGAGACACCTGCACTCACCGCCATTACcg                                                                                                   <  1:370123/148‑1 (MQ=255)
                                                   accgcgTCCGCCTGCTGATTGGCCGGGAGGGCGGTTTATCGGCAGATAAAATTGCCATGACTGCCCGCTATCAATTTCCTGATATCCTGTTGGGCCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACa                                                                                            <  1:524882/146‑1 (MQ=255)
                                                      gcgcCCGCCGGCTGATGGGCCGGGAAGGCGTTTTATCGGCAAATGAATTGCCCCTGCCGGCCCAATATCATTTTGCGGATATCCTGTTGGGGCCGCGGGTTTTGGGTACAGAGACACTTGCGCCCACCGCCATTACCGCGCTACAAGTAc                                                                                       <  1:94589/150‑1 (MQ=255)
                                                      gcgTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATAAAATTGCCATGACTGCCCGCTATCATTTTACTGATATCCTGTTGGGGCCTCGCTTTTTGCGTACAGAGAAAACTGCGCTCACCGCCATTACCGCGCTACAAGTAc                                                                                       <  2:209172/150‑1 (MQ=255)
                                                        gTCCCCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGAACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcga                                                                                     <  1:148359/150‑1 (MQ=255)
                                                          ccGCCTGCTGATTGGCCCGGAAGGCGTTTTATCGGCAGATGAAATTCCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatt                                                                                   <  1:176973/150‑1 (MQ=255)
                                                                            ggAGGGCGGTTTAGCGAAAGAAAAATTTGACATACCCGCCCGCTTCCATTTTACGATTGCCCTTTTGGGGTCTGGCTTTTGGCGTGCAGAGAAAACGGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTa                                                                    <  1:595761/147‑1 (MQ=255)
                                                                                ggCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGgaga                                                             >  2:173665/1‑150 (MQ=255)
                                                                                        aTCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATg                                                    >  2:198320/1‑151 (MQ=255)
                                                                                        aTCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATg                                                    >  2:231314/1‑151 (MQ=255)
                                                                                              aGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCa                                                >  1:33265/1‑149 (MQ=255)
                                                                                                  gAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCg                                          >  2:382067/1‑151 (MQ=255)
                                                                                                   aaaTTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTGGGGGCCTCGCTTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCaa                                               <  2:592686/145‑1 (MQ=255)
                                                                                                         ccATGACTGCCCGCTATCAATTTACTGATATCCTGTGGGGCCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCg                                    <  1:324670/150‑1 (MQ=255)
                                                                                                            tGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGCGAAGAATAATGATCAAGCTCGGCATCGTGa                                 >  1:472509/1‑150 (MQ=255)
                                                                                                               cTGCCCGCTATCAATTTACTGATATCCTGTGGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATgg                              <  2:149982/150‑1 (MQ=255)
                                                                                                                  cccGCTATCAATTTACTGATATCCTGTGGGGGCCTCGCGTTTTGCGTAAAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGAc                            <  1:194492/149‑1 (MQ=255)
                                                                                                                       tATCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCAt                       <  2:10072/149‑1 (MQ=255)
                                                                                                                        ataaaTTTACTGATACCCTGTTGGGGCCTCGGGTTTGGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCATc                      <  1:609509/146‑1 (MQ=255)
                                                                                                                        aTCAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCATCGc                    >  2:394114/1‑151 (MQ=255)
                                                                                                                          cAATTTACTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCATCGCa                   >  1:410948/1‑150 (MQ=255)
                                                                                                                           attttACTTATACCCTGTTGGGGCCTGGCGTTTTGGGTACAGAGACAACTGCGCCCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCATc                      <  1:261695/144‑1 (MQ=255)
                                                                                                                                 cTGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCATCGCAaacatca            >  2:454172/1‑150 (MQ=255)
                                                                                                                                  tGATATCCTGTTGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCATCGCAaacatcaac          >  1:265906/1‑151 (MQ=255)
                                                                                                                                    atatCCTGTTGGGGCCTCGCGTTTTGCGTACAGAAACAACTGCCCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCATCGCAaacatca            <  1:454167/147‑1 (MQ=255)
                                                                                                                                             tGGGGCCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCATCGCAAACATCAACATCaagaa  >  2:416106/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTGGCGATTTGGGCTAACGGAGAAGAATAATGATCAAGCTCGGCATCGTGATGGACCCCATCGCAAACATCAACATCAAGAA  >  NC_000913/3091642‑3091930

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