Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,410,053 (T)9→8 intergenic (+47/‑80) rlmB → / → yjfI 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,410,0450T.94.5% 216.8 / 12.5 57intergenic (+39/‑88)rlmB/yjfI23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent/DUF2170 family protein
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base . (29/23);  total (34/23)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAT  >  NC_000913/4409947‑4410139
                                                                                                  |                                                                                              
cTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtt                                                                                               >  2:776830/1‑100 (MQ=255)
 tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttt                                                                                              >  2:161460/1‑100 (MQ=255)
 tGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttt                                                                                              >  2:545040/1‑100 (MQ=255)
   aaCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGctttt                                                                                            <  2:1476797/100‑1 (MQ=38)
    aCGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGtttttt                                                                                           <  2:1533500/100‑1 (MQ=255)
     cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                          >  2:113089/1‑100 (MQ=255)
     cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttt                                                                                          <  1:315673/100‑1 (MQ=255)
     cGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttn                                                                                          <  2:1398282/100‑2 (MQ=255)
       tttCGGTTGCGACCGTAATTTGCTTATTTGATGCTGTGCGCCAGCTCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                        <  2:1083463/100‑2 (MQ=255)
       tttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                        <  1:41280/100‑2 (MQ=255)
        ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgt                                                                                       <  1:380413/100‑3 (MQ=255)
        ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttg                                                                                        <  1:374287/99‑2 (MQ=255)
        ttCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAACCATCCAAATCTGGATGGcttttcataa                                                                                       <  2:298822/100‑6 (MQ=37)
         tCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                      <  2:1553353/100‑4 (MQ=255)
         tCGGTTGCGACCGGAATTTGCTTAGTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGttttttttgtc                                                                                      <  1:1118265/100‑4 (MQ=255)
           ggTTGCGACCGTAAGTTGCTTATTTAAAGCGGTTCGCCATCGCAGCTAATTTCTCATAATTATGAAGAGCCATCCAGATTTGGATGGTTTTTTTTGTCta                                                                                   <  1:1088589/100‑1 (MQ=255)
                cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                              >  2:1622706/1‑100 (MQ=255)
                cGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCt                                                                              >  2:343683/1‑100 (MQ=255)
                 gACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTg                                                                             >  2:840717/1‑100 (MQ=255)
                   ccGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGt                                                                           >  2:115064/1‑100 (MQ=255)
                    cGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTa                                                                          <  1:790430/100‑1 (MQ=255)
                      gAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAg                                                                        <  1:963939/100‑1 (MQ=255)
                        aTTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGAt                                                                      >  2:328797/1‑100 (MQ=255)
                           tGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAAt                                                                   <  1:47405/100‑1 (MQ=255)
                            gCTTATTTGAAGCGGTTCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAAtt                                                                  <  2:1762386/100‑1 (MQ=255)
                             cTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTa                                                                 <  2:901339/100‑1 (MQ=255)
                             attATTTGAAGCGGTGCGCAAACGCAGCTAATTTCTCAGAATTAAGAAAAGACATCCAGATTTGGATGATTTTTTTTGTCTATAGATGGTAAGATAATTa                                                                 >  1:1252044/2‑100 (MQ=255)
                               tATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACg                                                               >  2:419653/1‑100 (MQ=255)
                                aTTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGt                                                              >  2:81484/1‑100 (MQ=255)
                                 tttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATATCTGGTAAGATAATTACGTa                                                             >  1:787703/1‑100 (MQ=255)
                                  ttGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTAt                                                            >  2:191059/1‑100 (MQ=255)
                                      aGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                        >  2:8056/1‑100 (MQ=255)
                                      aGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCa                                                        >  1:1510195/1‑100 (MQ=255)
                                          gTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAata                                                    >  2:1453049/1‑100 (MQ=255)
                                            gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                                  >  2:20702/1‑100 (MQ=255)
                                            gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                                  <  1:1622695/100‑1 (MQ=255)
                                            gcgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATAtt                                                  >  2:278546/1‑100 (MQ=255)
                                             cgcCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTc                                                 >  2:1369504/1‑100 (MQ=255)
                                             cgcCAGCGCAGCTAATTTCTCAGAAATATGAAAAGCCATCCAGATTTAGATGGTTTTTTTTTTCTATCGCTCATAAGATCATTACGTTCTCCCAATATTc                                                 >  1:1620764/1‑100 (MQ=255)
                                             cgcCAGCGCAGCTAATTTCTAAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAatat                                                   >  1:364905/1‑98 (MQ=255)
                                                 agtgcAGCTAATTTCTCAGAAAAAAGAAAAGCAAACCAGTTTTGTTTGATTTTATTTGTATATAACTAGTAAGATACTTACGTATTTAAAat                                                     >  1:711236/4‑92 (MQ=255)
                                                      aGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTTGTAAGATAATTACGTATTGCAAATATTCAtt                                              <  1:1383451/94‑1 (MQ=255)
                                                           aTTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAtt                                   >  2:937257/1‑100 (MQ=255)
                                                            tttCTCAGAATTATGAAAAGCCATCCTGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATAttt                                  <  2:1373830/100‑1 (MQ=255)
                                                                tcAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAAc                              <  1:841854/100‑1 (MQ=255)
                                                                 cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                             >  2:12355/1‑100 (MQ=255)
                                                                 cAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACt                             <  1:158410/100‑1 (MQ=255)
                                                                  aGAATTATGAAAAGCCATCCAGAAATGGATGGATTTTTTTGTTTATATCTGGTAAGATAATTAAGTATTTCCAATATTCATTTTATTAATATTTAAACtt                            >  1:709226/1‑100 (MQ=255)
                                                                       tATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAATTTAATCt                       <  2:956107/100‑1 (MQ=255)
                                                                         aaaaaaACCATCCAAATCTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCta                      >  2:843228/3‑100 (MQ=39)
                                                                          gAAAAGCCATCCAGATTTGTATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                    <  1:336043/100‑1 (MQ=255)
                                                                          gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                    >  1:1292994/1‑100 (MQ=255)
                                                                          gAAAAGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                    >  2:252680/1‑100 (MQ=255)
                                                                          gAAAAGCCATCCAGATTTGGATGGTTTTTTTTCTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCtata                    >  2:1724348/1‑100 (MQ=255)
                                                                          gAAAAGCCATCCAGATTTGGAAGGATTTTTTTGTATATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATTta                      >  1:855442/1‑98 (MQ=255)
                                                                           aaaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTatat                   >  2:1523593/1‑100 (MQ=255)
                                                                            aaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                  <  1:343681/100‑1 (MQ=255)
                                                                            aaaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATg                  <  2:263807/100‑1 (MQ=255)
                                                                             aaGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATgg                 <  1:161460/100‑1 (MQ=255)
                                                                              aGCCATCCAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGa                <  1:1506184/100‑1 (MQ=255)
                                                                                 cATCCAGATTTGGATGGATTTATTAGTCTATAGCTGGTAAGATAATAACGTATCGCAAATATTCATTTTATTTATATTTAATCTTAATCTATATGGaata             >  1:760583/1‑100 (MQ=255)
                                                                                     cAGATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAataata          <  1:1106666/99‑1 (MQ=255)
                                                                                      agttttGGATTGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata        <  2:1735977/97‑1 (MQ=255)
                                                                                      aaatctGGATGGTTTTTTTTGTTTATATCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatata        >  1:65258/6‑100 (MQ=255)
                                                                                       gATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat       <  1:780928/100‑1 (MQ=255)
                                                                                       gATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat       >  2:1263347/1‑100 (MQ=255)
                                                                                       gATTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAatatat       <  1:12355/100‑1 (MQ=255)
                                                                                        aTTTGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg      <  2:692224/100‑1 (MQ=255)
                                                                                         tttGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa     >  2:664247/1‑100 (MQ=255)
                                                                                         tttGGATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGa     <  1:1688537/100‑1 (MQ=255)
                                                                                          ttttaatatttttttttGTCTATAGCTAGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTTATCTATATGGAATAATATATgc     <  1:762528/92‑2 (MQ=255)
                                                                                          ttGGATGGTTTTTTTTTTTTTTATCTGATAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATg      >  1:768108/1‑98 (MQ=255)
                                                                                            ggATGGTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGAAATAATATATGAAAt  >  1:1740201/1‑100 (MQ=255)
                                                                                                  |                                                                                              
CTGAACGTTTCGGTTGCGACCGGAATTTGCTTATTTGAAGCGGTGCGCCAGCGCAGCTAATTTCTCAGAATTATGAAAAGCCATCCAGATTTGGATGGTTTTTTTTTGTCTATAGCTGGTAAGATAATTACGTATTGCAAATATTCATTTTATTAATATTTAAACTTAATCTATATGGAATAATATATGACAT  >  NC_000913/4409947‑4410139

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