Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 49,644 T→C 40.2% K155K (AAA→AAG ksgA ← S‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE49,6440TC40.2% 66.4 / 112.0 154K155K (AAA→AAGksgAS‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (69/23);  new base C (18/44);  total (87/67)
Fisher's exact test for biased strand distribution p-value = 2.38e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

Only 100 of 154 total aligned reads displayed.

GGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGC  >  minE/49594‑49709
                                                  |                                                                 
gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATTCGCTTTGCTGTTCGGt                                                       <  1:413213/63‑1 (MQ=255)
gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGt                                                       <  1:2124983/63‑1 (MQ=255)
gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGt                                                       <  1:2198403/63‑1 (MQ=255)
gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGt                                                       <  1:231356/63‑1 (MQ=255)
gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGt                                                       <  1:3211533/63‑1 (MQ=255)
gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGt                                                       <  1:2890955/63‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:3090692/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:230960/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:2418390/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:2487608/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:2543669/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:2730536/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:3004064/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:1039260/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:3262452/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:3338927/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:3409274/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:544495/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:572534/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:705668/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:757952/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:976342/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:144890/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:1094665/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:1201445/69‑1 (MQ=255)
    tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:1553924/69‑1 (MQ=255)
     cacgTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:108978/68‑1 (MQ=255)
     cacgTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa                                             <  1:359533/68‑1 (MQ=255)
           gCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGAtt                                      >  1:2991190/1‑69 (MQ=255)
           gCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGAtt                                      >  1:1345286/1‑69 (MQ=255)
                aGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacca                                 >  1:3042113/1‑69 (MQ=255)
                aGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacca                                 >  1:1739034/1‑69 (MQ=255)
                aGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacca                                 >  1:653210/1‑69 (MQ=255)
                aGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacc                                  >  1:1008835/1‑68 (MQ=255)
                 gTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcaccac                                <  1:2294546/69‑1 (MQ=255)
                 gTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcaccac                                <  1:3302852/69‑1 (MQ=255)
                 gTATTGCGCCATGACGCTTAATAGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcaccac                                >  1:2194128/1‑69 (MQ=255)
                      gcgcCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGa                                        <  1:1927108/56‑1 (MQ=255)
                      gcgcCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGa                                        <  1:274593/56‑1 (MQ=255)
                          cATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCttttt                        >  1:1557931/1‑68 (MQ=255)
                          cATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCttttt                        >  1:967326/1‑68 (MQ=255)
                          cATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCttttt                        >  1:553018/1‑68 (MQ=255)
                          cATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCttttt                        >  1:440101/1‑68 (MQ=255)
                           aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTg                       >  1:2143825/1‑68 (MQ=255)
                           aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                      >  1:2433937/1‑69 (MQ=255)
                           aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                      >  1:39457/1‑69 (MQ=255)
                           aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                      >  1:2535788/1‑69 (MQ=255)
                           aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                      >  1:1682682/1‑69 (MQ=255)
                           aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                      >  1:1883802/1‑69 (MQ=255)
                           aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                      >  1:482479/1‑69 (MQ=255)
                           aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                      >  1:1140362/1‑69 (MQ=255)
                           aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                      >  1:3326096/1‑69 (MQ=255)
                           aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc                      >  1:1536620/1‑69 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGc                                                  >  1:1403197/1‑38 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACctc                             >  1:2200310/1‑59 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa                    >  1:3153286/1‑68 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa                    >  1:3363671/1‑68 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa                    >  1:3428436/1‑68 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa                    >  1:3053471/1‑68 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa                    >  1:2968023/1‑68 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa                    >  1:1835420/1‑68 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa                    >  1:2014897/1‑68 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa                    >  1:2068774/1‑68 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa                    >  1:2180163/1‑68 (MQ=255)
                              aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa                    >  1:150572/1‑68 (MQ=255)
                                 cTTAATCGACCATACGCTTTGCTTTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATa                <  1:2967161/69‑1 (MQ=255)
                                 cTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATa                <  1:3242411/69‑1 (MQ=255)
                                       cGACCATACGCTTTGCTGTTCGGTCCTGCACCCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg           <  1:1306263/68‑1 (MQ=255)
                                       cGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg           <  1:2013358/68‑1 (MQ=255)
                                       cGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg           <  1:901566/68‑1 (MQ=255)
                                       cGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg           <  1:1486618/68‑1 (MQ=255)
                                       cGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg           <  1:609941/68‑1 (MQ=255)
                                         aCCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATa                <  1:1218013/61‑1 (MQ=255)
                                          ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc                             >  1:3144255/1‑47 (MQ=39)
                                          ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc                             >  1:2506161/1‑47 (MQ=39)
                                          ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc                             >  1:2640117/1‑47 (MQ=39)
                                          ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc                             >  1:3134878/1‑47 (MQ=39)
                                          ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc                             >  1:1090053/1‑47 (MQ=39)
                                          ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc                             >  1:2806389/1‑47 (MQ=39)
                                          ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc                             >  1:2842489/1‑47 (MQ=39)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                           >  1:3413462/1‑44 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                           >  1:745180/1‑44 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                           >  1:2438474/1‑44 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc  >  1:2416337/1‑69 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc  >  1:130487/1‑69 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc  <  1:506250/69‑1 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc  >  1:2578115/1‑69 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc  >  1:2763902/1‑69 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc  >  1:1230853/1‑69 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc  >  1:1234797/1‑69 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc  >  1:3262207/1‑69 (MQ=255)
                                               cGCTTTGCTGTTCGGTCCTGAAACCAGACGATTCACCACCTCtt                           >  1:370581/1‑44 (MQ=255)
                                               cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc  <  1:313708/69‑1 (MQ=255)
                                               cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc  <  1:2156663/69‑1 (MQ=255)
                                               cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc  <  1:591138/69‑1 (MQ=255)
                                               cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc  <  1:3307524/69‑1 (MQ=255)
                                               cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc  <  1:886154/69‑1 (MQ=255)
                                               cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc  <  1:948451/69‑1 (MQ=255)
                                               cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc  <  1:1391572/69‑1 (MQ=255)
                                               agCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt                           >  1:2925812/2‑44 (MQ=255)
                                                  |                                                                 
GGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGC  >  minE/49594‑49709

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