Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,263,269 T→C 44.9% Q80Q (CAA→CAG rplA ← 50S ribosomal subunit protein L1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,263,2690TC44.9% 14.1 / 73.9 89Q80Q (CAA→CAGrplA50S ribosomal subunit protein L1
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base C (19/21);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

CAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGTTGCACCACG  >  minE/2263205‑2263331
                                                                |                                                              
cAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTgg                                                             >  1:3060086/1‑68 (MQ=255)
cAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTgg                                                             >  1:3228753/1‑68 (MQ=255)
cAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTgg                                                             >  1:501751/1‑68 (MQ=255)
cAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTgg                                                             >  1:1722264/1‑68 (MQ=255)
cAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTgg                                                             >  1:1642942/1‑68 (MQ=255)
cAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTgg                                                             >  1:1377808/1‑68 (MQ=255)
cAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTgg                                                             >  1:1227692/1‑68 (MQ=255)
cAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTg                                                              >  1:1050826/1‑67 (MQ=255)
cAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTg                                                              >  1:353336/1‑67 (MQ=255)
cAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGAGTTTGCACCTTgg                                                             >  1:1090727/1‑68 (MQ=255)
   ccAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTa                                                          <  1:634333/68‑1 (MQ=255)
    cAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTaaa                                                        <  1:93831/69‑1 (MQ=255)
    cAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTaaa                                                        <  1:556559/69‑1 (MQ=255)
    cAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTaaa                                                        <  1:599145/69‑1 (MQ=255)
    cAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTaaa                                                        <  1:2322458/69‑1 (MQ=255)
    cAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTaaa                                                        <  1:1792361/69‑1 (MQ=255)
    cAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTaaa                                                        <  1:1871894/69‑1 (MQ=255)
    cAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTaaa                                                        <  1:1381561/69‑1 (MQ=255)
    cAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTaaa                                                        <  1:3244365/69‑1 (MQ=255)
    cAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTaaa                                                        <  1:822340/69‑1 (MQ=255)
     aGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGTGTa                                                          <  1:381432/66‑1 (MQ=255)
        tCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACg                                                    >  1:3395368/1‑69 (MQ=255)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:1954318/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:527564/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:1287770/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:1974849/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:1174630/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:2326567/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:2381610/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:2518060/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:351930/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:3351404/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:1947676/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:1341713/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:822361/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:1406386/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:819446/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:782693/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:1510994/1‑66 (MQ=38)
            ccATACCTACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:1197461/1‑66 (MQ=38)
            ccATACATACCAGTTCAGCACCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCCTGGGTAAATACgg                                                   >  1:4728/1‑66 (MQ=37)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:2338463/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:235917/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:3020236/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:2417953/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:2970115/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:1029714/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:1639735/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:1699078/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:2589711/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:321082/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:1790761/1‑69 (MQ=255)
                 cctaccAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCga                                           >  1:1898423/1‑69 (MQ=255)
                       aGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCgaacgg                                       <  1:193650/67‑1 (MQ=255)
                       aGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCgaacgg                                       <  1:2730652/67‑1 (MQ=255)
                       aGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCgaacgg                                       <  1:1206960/67‑1 (MQ=255)
                       aGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCgaacgg                                       <  1:1667336/67‑1 (MQ=255)
                             gcgcCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGcc                                <  1:2482362/68‑1 (MQ=255)
                             gcgcCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGcc                                <  1:1321310/68‑1 (MQ=255)
                             gcgcCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGcc                                <  1:2208872/68‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:3114096/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:1170789/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:829836/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:649385/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:16967/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:1740841/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:207667/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:513702/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:2189332/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:432679/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:39614/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:2669841/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:2869215/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:2877345/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:3345973/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:3066005/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:3187385/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:3264117/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:2906126/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGtt          <  1:2924595/69‑1 (MQ=255)
                                                  tCAGCGTTTGCACCCTGGGTAAATACGGCTACGCGAACGGAACGGCCAGCACCGTGCGGCAGTACAGtt          <  1:2472345/69‑1 (MQ=255)
                                                     gCGTTTGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGTTg         <  1:1496767/67‑1 (MQ=255)
                                                        tttGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGTTGcacca    <  1:2899949/69‑1 (MQ=255)
                                                        tttGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGTTGcacca    <  1:972567/69‑1 (MQ=255)
                                                        tttGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGTTGcacca    <  1:1211087/69‑1 (MQ=255)
                                                        tttGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGTTGcacca    <  1:556595/69‑1 (MQ=255)
                                                        tttGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGTTGcacca    <  1:2579434/69‑1 (MQ=255)
                                                         ttGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGTTGcacca    <  1:1472681/68‑1 (MQ=255)
                                                          tGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACAGTGCGGCAGTACAGTTGCACCacg  >  1:1768427/1‑69 (MQ=255)
                                                                |                                                              
CAGCCAGATCTTCCATACCTACCAGTTCTGCGCCTGCAGCTTTAGCAGCTTCAGCGTTTGCACCTTGGGTAAATACGGCTACGCGAACGGAACGGCCAGTACCGTGCGGCAGTACAGTTGCACCACG  >  minE/2263205‑2263331

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