Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,091,709 G→A 13.5% L192L (CTG→CTA)  rsmE → 16S rRNA m(3)U1498 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,091,7090GA13.5% 75.7 / 4.1 36L192L (CTG→CTA) rsmE16S rRNA m(3)U1498 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (17/14);  new base A (3/2);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

ATCCGTCCAGCGATGGATCTGGAAGCCTGGTGTGCAGAGCAGGATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTG  >  NC_000913/3091581‑3091851
                                                                                                                                |                                                                                                                                              
aTCCGTCCAGCGATGGATCTGGAAGCCTGGTGTGCAGAGCAGGATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAgg                                                                                                                                  <  1:341870/143‑1 (MQ=255)
  ccGTGCAGCGATGGATCTGGAAGCCTGGTGTGCAGAGCAGGATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCg                                                                                                                                <  1:577619/143‑1 (MQ=255)
    gTCCAGCGATGGATCTGGAAGCCTGGTGTGCAGAGCAGGATGAAGGACTGAACCTGAATCTTCACCCGCGCGCCAGTAACAGCCTCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCGCCGAATGGGCCCGGAAGGCgag                                                                                                                              >  2:111663/1‑141 (MQ=255)
            aTGGATCTGGAAGCCTGGTGTGCAGAGCAGGATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTAATTGGCCCGGAAGGCGGTTTATCGGc                                                                                                                      >  1:154594/1‑143 (MQ=255)
            aTGGATCTGGAAGCCTGGTGTGCAGAGCAGGATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTAATTGGCCCGGAAGGCGGTTTATCGGc                                                                                                                      >  1:288186/1‑143 (MQ=255)
                  cTGGAAGCCTGGTGTGCAGAGCAGGATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACCCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGAtcc                                                                                                                >  1:410373/1‑141 (MQ=255)
                   tGGAAGCCTGGTGTGCAGAGCAGGATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTAGCGGCAGATGaa                                                                                                               >  1:182796/1‑143 (MQ=255)
                   gtggaGCGTGGTGTGCAGAGAACGATGAAGGAGTGATAGTGAGTCTTCAGCTGCGCGCCAGTAACAGCCTCAATACGTTGCCGTTACCGGTTGAACGCGGCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGaa                                                                                                               <  2:281542/139‑1 (MQ=255)
                        gCCTGGTGTGCAGAGCAGGATGAAGGGCTGAAACTCAGTCTTCACCCGCGCGCCAGTAACAGCATCAATACGTGGCCGTTACCGGTTGAACGCGTCCGCCTGATGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGc                                                                                                          <  2:105400/143‑1 (MQ=255)
                           tGGTGTGCAGAGCAGGATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCAt                                                                                                       >  2:799087/1‑143 (MQ=255)
                                 gCAGAGCAGGATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGCTCCCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTCCCATGACTGc                                                                                                 >  2:155514/1‑143 (MQ=255)
                                      gCAGGAAGAAGGACTGAAACTGAAGCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGTTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCt                                                                                            <  1:107449/143‑1 (MQ=255)
                                         ggATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTCTCGGCAGATGAAATTGCCATGACTGACTGGTATc                                                                                         >  1:158447/1‑143 (MQ=255)
                                          gATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCGGCTATCa                                                                                        >  1:11724/1‑143 (MQ=255)
                                                 gACTTAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACGGCCCGCTATCAATTTACt                                                                                 >  2:372283/1‑143 (MQ=255)
                                                      aaactgaaTCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCCGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGatat                                                                            <  2:484473/143‑1 (MQ=255)
                                                         ctgaaTCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCt                                                                         >  1:773742/1‑143 (MQ=255)
                                                          tgaaTCGGCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTg                                                                        <  2:430444/143‑1 (MQ=255)
                                                          ggaaTCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTg                                                                        <  1:876780/142‑1 (MQ=255)
                                                              tCTTCACCCGCGCGGCAGTAACAGCAGCAATACGTGGCCGTGACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTgg                                                                    <  1:568238/143‑1 (MQ=255)
                                                              tCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTACCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGCAACTGCCATGACTTCCCGCCATCACTTTACTGATATGCTGTTgg                                                                    >  1:579494/1‑143 (MQ=255)
                                                                    cccGCGCGCCAGGAACAGCATCATTACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTc                                                              <  2:800827/143‑1 (MQ=255)
                                                                       gcgcgcCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTcgcg                                                           <  1:382655/143‑1 (MQ=255)
                                                                                 aaCAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGGCCTCGGGTTTGGCGtcc                                                 >  1:216070/1‑141 (MQ=255)
                                                                                  aCAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTCATATCCTGTTGGGACCTCGCGTTTTGCGCAcc                                                >  1:720361/1‑142 (MQ=255)
                                                                                  aCAGCATCAATACGTTGCCGTTACCGGGTGAACGCGTCCGCGTGCTGATTGGCCCGGAAGGCGGTTTATCTGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACa                                                <  2:408877/143‑1 (MQ=255)
                                                                                    aGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGAATACAGGCAGTTGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACaga                                              >  1:800661/1‑143 (MQ=255)
                                                                                              cGTTGCCGTTACCGGTTGGACGCGGCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTgcg                                    <  2:519429/143‑1 (MQ=255)
                                                                                                  gCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGCGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCa                                >  2:512603/1‑143 (MQ=255)
                                                                                                        aCCGGTTGAACGCGTCCGCCTGCTAATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCa                          <  2:288186/143‑1 (MQ=255)
                                                                                                        aCCGGTTGAACGCGTCCGCCTGCTAATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCa                          <  2:154594/143‑1 (MQ=255)
                                                                                                                agcgctTCTGCCTGCTGCCTGGCCCGGTCGGCGGTTTATCGGCAGATGACATTGCCATGATTTCCCGCTCGCAATTTACTCATATCCTGTTGTGTCCTCGCGTTTTGCGTACAGAGACACCTGCGCTCACCGCCATTagctcc                  <  2:833915/141‑6 (MQ=255)
                                                                                                                aaCGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACcgcg                  >  2:87406/1‑143 (MQ=255)
                                                                                                                 aCGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCgcgc                 >  2:343246/1‑143 (MQ=255)
                                                                                                                    cgTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCACTTTACTGATATCCTGTTGGGACCTCGCGGTTTGCGTACAGAGACAACTGCGCGCACCGCCAGCACCGGGCTAc              >  1:732850/1‑143 (MQ=255)
                                                                                                                     gTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCCTTTCcgcgttttt             >  1:322992/1‑138 (MQ=255)
                                                                                                                          tctgctgAGTGGCCCGGATGGCGGTTTATCGGAGGATGAAATTGCCATGACTGCGCGCTGTCAATTTACTGATATGCTGTTGGGACCTCGCGTTTTGCGTCCAGTGACTGCTGCGCTCTCCGCCATTACCGCGCTACAAGTAc        <  2:780208/142‑1 (MQ=255)
                                                                                                                                gATTGGCCCGGAAGGCGGTTTCTCTGCAGATGAAATTGCCCTGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatttg  <  2:78240/143‑1 (MQ=255)
                                                                                                                                gATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTAcgatttg  <  2:158447/143‑1 (MQ=255)
                                                                                                                                |                                                                                                                                              
ATCCGTCCAGCGATGGATCTGGAAGCCTGGTGTGCAGAGCAGGATGAAGGACTGAAACTGAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCCTGTTGGGACCTCGCGTTTTGCGTACAGAGACAACTGCGCTCACCGCCATTACCGCGCTACAAGTACGATTTG  >  NC_000913/3091581‑3091851

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑