Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,061,553 G→A 5.7% K20K (AAG→AAA rpsU → 30S ribosomal subunit protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,061,5530GA 5.7% 290.5 / 5.6 87K20K (AAG→AAArpsU30S ribosomal subunit protein S21
Reads supporting (aligned to +/- strand):  ref base G (46/36);  new base A (3/2);  total (49/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCG  >  minE/2061485‑2061623
                                                                    |                                                                      
gagaGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCg                                                                      >  1:3287774/1‑71 (MQ=255)
gagaGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCg                                                                      >  1:3612611/1‑71 (MQ=255)
gagaGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGACGTTCGACGTAGCTCTGCGTCGCTTCAAGCg                                                                      >  1:625022/1‑71 (MQ=255)
 agagGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGt                                                                     >  1:3955954/1‑71 (MQ=255)
 agagGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGt                                                                     >  1:1975976/1‑71 (MQ=255)
  gagGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGtt                                                                    >  1:2975369/1‑71 (MQ=255)
   agGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTc                                                                   >  1:1294462/1‑71 (MQ=255)
       acaTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGc                                                               <  1:1018230/71‑1 (MQ=255)
       acaTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGc                                                               <  1:932962/71‑1 (MQ=255)
         aTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGa                                                             <  1:4697634/71‑1 (MQ=255)
         aTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGa                                                             <  1:2514086/71‑1 (MQ=255)
         aTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGa                                                             <  1:4650931/71‑1 (MQ=255)
         aTGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGa                                                             <  1:1768505/71‑1 (MQ=255)
          tGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGa                                                             >  1:542875/1‑70 (MQ=255)
           gCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaa                                                            >  1:281357/1‑70 (MQ=255)
            ccGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaaa                                                          <  1:2203690/71‑1 (MQ=255)
            ccGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaa                                                           <  1:2549376/70‑1 (MQ=255)
            ccGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaa                                                           >  1:1994760/1‑70 (MQ=255)
            ccGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaa                                                           <  1:622552/70‑1 (MQ=255)
             cGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaa                                                           >  1:720970/1‑69 (MQ=255)
             cGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaa                                                           >  1:2831723/1‑69 (MQ=255)
                   ttAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAgg                                                    <  1:1864572/70‑1 (MQ=255)
                   ttAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgt                                                   >  1:711071/1‑71 (MQ=255)
                    tAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtg                                                  >  1:1850508/1‑71 (MQ=255)
                    tAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtg                                                  >  1:307214/1‑71 (MQ=255)
                     aaaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                 >  1:364969/1‑71 (MQ=255)
                     aaaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                 >  1:4416356/1‑71 (MQ=255)
                     aaaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                 >  1:3040003/1‑71 (MQ=255)
                      aaGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                 <  1:2317958/70‑1 (MQ=255)
                           cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGc                                                       >  1:127102/1‑59 (MQ=255)
                           cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGc                                                       >  1:2404828/1‑59 (MQ=255)
                           cGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTgg                                            <  1:1992416/70‑1 (MQ=255)
                             tGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                 <  1:2518852/63‑1 (MQ=255)
                             tGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                 <  1:1463528/63‑1 (MQ=255)
                             tGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                 <  1:880993/63‑1 (MQ=255)
                              gAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCgg                                         >  1:2811013/1‑70 (MQ=255)
                               aaaaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGgtgt                                                 <  1:2022568/61‑1 (MQ=255)
                               aaaaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCgg                                         >  1:285253/1‑69 (MQ=255)
                               aaaaCGAGCCGTTCGAAGTCGCAATGCGTCGCTTCAAACGCACCATCGAAAAAACTGGTCTGCTGACGGaa                                       >  1:1449614/1‑71 (MQ=255)
                               aaaaCGAGCCGTTCGAAGTCGCAATGCGTCGCTTCAAACGCACCATCGAAAAAACTGGTCTGCTGACGGaa                                       >  1:1019933/1‑71 (MQ=255)
                               aaaaCGAGCCGTTCGAAGTCGCAATGCGTCGCTTCAAACGCACCATCGAAAAAACTGGTCTGCTGACGGa                                        >  1:432435/1‑70 (MQ=255)
                                aaaCGAGCCGTTCGAAGTCGCAATGCGTCGCTTCAAACGCACCATCGAAAAAACTGGTCTGCTGACGGaac                                      <  1:3206717/71‑2 (MQ=255)
                                aaaCGAGCCGTTCGAAGTCGCAATGCGTCGCTTCAAACGCACCATCGAAAAAACTGGTCTGCTGACGGaac                                      <  1:1410787/71‑2 (MQ=255)
                                 aaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGt                                     <  1:987020/71‑1 (MQ=255)
                                 aaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGt                                     <  1:4660032/71‑1 (MQ=255)
                                 aaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGt                                     <  1:4741940/71‑1 (MQ=255)
                                 aaCGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGt                                     >  1:283093/1‑71 (MQ=255)
                                     aGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaaaa                                                         >  1:217485/1‑47 (MQ=255)
                                     aGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaaaaa                                                         >  1:4322648/1‑47 (MQ=255)
                                         gTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                             <  1:1336043/71‑1 (MQ=255)
                                         gTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                             <  1:1479754/71‑1 (MQ=255)
                                          ttCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGaa                                                            <  1:4548268/39‑1 (MQ=255)
                                            cGACGTAGCTCTGCGTCGCTTCAAGCTTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                           >  1:4333830/1‑70 (MQ=255)
                                            cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTc                                               >  1:375935/1‑50 (MQ=255)
                                            cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGaa                                       >  1:544467/1‑58 (MQ=255)
                                            cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGaa                                       >  1:353370/1‑58 (MQ=255)
                                            cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGaa                                       >  1:2890048/1‑58 (MQ=255)
                                            cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGaa                                       >  1:946537/1‑58 (MQ=255)
                                            cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                           >  1:2437457/1‑70 (MQ=255)
                                            cGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                           <  1:3036209/70‑1 (MQ=255)
                                             gACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTg                          <  1:4286524/70‑1 (MQ=255)
                                                 tAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGtt                                                <  1:2920429/44‑1 (MQ=255)
                                                  aGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGt                       >  1:3257521/1‑68 (MQ=255)
                                                  aGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTc                     >  1:1869581/1‑70 (MQ=255)
                                                   gCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCg                                          <  1:3799372/48‑1 (MQ=255)
                                                   gCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTa                   <  1:1819439/71‑1 (MQ=255)
                                                   gCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGAGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCt                    <  1:4684452/70‑1 (MQ=255)
                                                     tctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTAt                  <  1:718194/70‑1 (MQ=255)
                                                      ctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                             <  1:2316137/58‑1 (MQ=255)
                                                      ctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                             <  1:1307541/58‑1 (MQ=255)
                                                      ctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                             <  1:1663780/58‑1 (MQ=255)
                                                      ctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                             <  1:173323/58‑1 (MQ=255)
                                                      ctGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTtcgtcgtc                             <  1:867195/58‑1 (MQ=255)
                                                         cGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa              >  1:2984343/1‑70 (MQ=255)
                                                             gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTcgtcgtcgt                           >  1:2857210/1‑53 (MQ=255)
                                                             gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTTGTGAGTTCTATGAAAAACCg         <  1:4737592/71‑1 (MQ=255)
                                                             gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa              >  1:4551479/1‑66 (MQ=255)
                                                             gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa              >  1:735017/1‑66 (MQ=255)
                                                             gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa              >  1:1500118/1‑66 (MQ=255)
                                                             gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa              >  1:4228349/1‑66 (MQ=255)
                                                             gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGaaa              >  1:3080315/1‑66 (MQ=255)
                                                             gCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGCAAAAcc          >  1:102703/1‑70 (MQ=255)
                                                               ttCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAcc          <  1:1057783/68‑1 (MQ=255)
                                                                 cAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTa     >  1:422836/1‑71 (MQ=255)
                                                                 cAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTa     >  1:1024470/1‑71 (MQ=255)
                                                                   aGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAAc           >  1:4381512/1‑63 (MQ=255)
                                                                   aGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTAcc   >  1:1270394/1‑71 (MQ=255)
                                                                    gCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCg  <  1:926129/71‑1 (MQ=255)
                                                                    |                                                                      
GAGAGGCACATGCCGGTAATTAAAGTACGTGAAAACGAGCCGTTCGACGTAGCTCTGCGTCGCTTCAAGCGTTCCTGCGAAAAAGCAGGTGTTCTGGCGGAAGTTCGTCGTCGTGAGTTCTATGAAAAACCGACTACCG  >  minE/2061485‑2061623

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: