Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,952,884 G→A 12.9% intergenic (‑193/+75) leuQ ← / ← rsmC tRNA‑Leu/16S RNA m2G1207 methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,952,8840GA12.9% 148.4 / 10.8 54intergenic (‑193/+75)leuQ/rsmCtRNA‑Leu/16S RNA m2G1207 methylase
Reads supporting (aligned to +/- strand):  ref base G (29/18);  new base A (3/4);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

GCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGC  >  minE/2952815‑2952951
                                                                     |                                                                   
gCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACg                                                                     >  1:2330371/1‑70 (MQ=255)
      aGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCt                                                              <  1:2009551/71‑1 (MQ=255)
           tGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAActt                                                         >  1:2109310/1‑69 (MQ=255)
             gTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGc                                                       <  1:2008351/71‑1 (MQ=255)
               atgcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAg                                                     >  1:1465149/1‑71 (MQ=255)
               atgcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAg                                                     >  1:2079663/1‑71 (MQ=255)
                tgcattcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAgg                                                                 <  1:2309817/58‑1 (MQ=255)
                  catgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGtt                                                  >  1:464906/1‑71 (MQ=255)
                        ggggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTa                                             <  1:2750459/70‑1 (MQ=255)
                            cagatgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAg                                            <  1:950931/67‑1 (MQ=255)
                                tgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGAt                                     >  1:1144213/1‑70 (MQ=255)
                                 gccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGAt                                     >  1:1388348/1‑69 (MQ=255)
                                 gccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGAt                                     >  1:1406971/1‑69 (MQ=255)
                                   cagatgcGACGCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGTGCAGGttttt                                               >  1:2483765/1‑56 (MQ=18)
                                   cagatgcGACGCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGTGCAGGttttt                                               >  1:2479270/1‑56 (MQ=18)
                                    agatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGc                                <  1:3136668/71‑1 (MQ=255)
                                      atgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGt                              >  1:3801982/1‑71 (MQ=255)
                                          gACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCa                           >  1:1922309/1‑70 (MQ=255)
                                          gACGCTGACGCGTCTTATCTGGCCTACAAAGGGCTa                                                             <  1:3618509/36‑1 (MQ=14)
                                          gACGCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGt                                                         >  1:3460426/1‑40 (MQ=17)
                                          gACGCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGt                                                         <  1:2633690/40‑1 (MQ=17)
                                            cGCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGt                                                         <  1:380235/38‑1 (MQ=17)
                                             gCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCAcgcc                       <  1:1581746/71‑1 (MQ=255)
                                             gCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGTGCAGGttt                                                 <  1:2356655/45‑1 (MQ=25)
                                              cTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGt                                              <  1:1876783/47‑1 (MQ=255)
                                                ggCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAg                                            >  1:3599125/1‑47 (MQ=255)
                                                 gcgcGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCa                    <  1:428120/70‑1 (MQ=255)
                                                  cgcgTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTa                                             <  1:1344986/44‑1 (MQ=255)
                                                  cgcgTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTa                                             <  1:1296708/44‑1 (MQ=255)
                                                   gcgTCTTATCTGGCCTACGAAGGGCTAACGTGCAg                                                     >  1:57537/1‑35 (MQ=255)
                                                   gcgTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATa                                    <  1:1204308/52‑1 (MQ=255)
                                                   gcgTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAg                                  >  1:2583045/1‑54 (MQ=255)
                                                     gTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa               >  1:1176388/1‑71 (MQ=255)
                                                     gTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa               >  1:1478570/1‑71 (MQ=255)
                                                     gTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa               >  1:2556783/1‑71 (MQ=255)
                                                      tCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGa                                      <  1:3059670/47‑1 (MQ=255)
                                                      tCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGa                                      <  1:377552/47‑1 (MQ=255)
                                                      tCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa               <  1:588979/70‑1 (MQ=255)
                                                        ttATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCg                                        >  1:2765398/1‑43 (MQ=255)
                                                        ttATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGacac            >  1:3879891/1‑71 (MQ=255)
                                                        ttATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGacac            >  1:3875302/1‑71 (MQ=255)
                                                          aTCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGtt                             >  1:1878419/1‑52 (MQ=255)
                                                            cTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGt                                          >  1:1270320/1‑37 (MQ=255)
                                                            cTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCAcgct                       <  1:1866687/56‑2 (MQ=255)
                                                            cTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTa        >  1:3728311/1‑71 (MQ=255)
                                                            cTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTa        >  1:3514752/1‑71 (MQ=255)
                                                             tGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTAt       >  1:2218443/1‑71 (MQ=255)
                                                               gCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGt                              <  1:944342/46‑1 (MQ=255)
                                                                ccTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTa        <  1:1741151/67‑1 (MQ=255)
                                                                ccTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCg    >  1:3437085/1‑71 (MQ=255)
                                                                 cTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCgg   >  1:2896109/1‑71 (MQ=255)
                                                                 cTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCgg   >  1:3428572/1‑71 (MQ=255)
                                                                 cTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCg    >  1:2653060/1‑70 (MQ=255)
                                                                   aCGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGc  >  1:1959222/1‑70 (MQ=255)
                                                                    cGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATaa                                   >  1:980247/1‑36 (MQ=255)
                                                                     |                                                                   
GCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGC  >  minE/2952815‑2952951

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