Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,883,198 A→G 26.4% F429L (TTC→CTC)  rumA ← 23S rRNA (uracil‑5)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,883,1980AG26.4% 187.4 / 79.7 110F429L (TTC→CTC) rumA23S rRNA (uracil‑5)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (36/45);  new base G (24/5);  total (60/50)
Fisher's exact test for biased strand distribution p-value = 4.20e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 138 total aligned reads displayed.

AACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATC  >  minE/1883128‑1883269
                                                                      |                                                                       
aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                         <  1:172995/71‑1 (MQ=255)
aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                         <  1:1913461/71‑1 (MQ=255)
aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                         <  1:2844726/71‑1 (MQ=255)
 aCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaaa                                                                        <  1:1763435/71‑1 (MQ=255)
 aCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaaa                                                                        <  1:3327997/71‑1 (MQ=255)
 aCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                         >  1:2698469/1‑70 (MQ=255)
    aTCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAg                                                                      >  1:3552685/1‑70 (MQ=255)
     tCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGt                                                                     <  1:2965868/70‑1 (MQ=255)
       gTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTa                                                                    >  1:2977156/1‑69 (MQ=255)
       gTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAc                                                                   >  1:3080079/1‑70 (MQ=255)
       gTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAc                                                                   <  1:2831944/70‑1 (MQ=255)
       gTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAc                                                                   >  1:1148648/1‑70 (MQ=255)
       gTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAc                                                                   >  1:3569135/1‑70 (MQ=255)
          ctctCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATc                                                               <  1:3601750/71‑1 (MQ=255)
          ctctCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATc                                                               <  1:3616077/71‑1 (MQ=255)
           tctcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                         <  1:2092425/60‑1 (MQ=255)
            ctcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaaaa                                                                       >  1:1004336/1‑61 (MQ=255)
            ctcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaaaa                                                                       >  1:2420987/1‑61 (MQ=255)
            ctcCTATAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaaaa                                                                       >  1:765367/1‑61 (MQ=255)
             tcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATc                                                               <  1:3147319/68‑1 (MQ=255)
             tcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGAt                                                            >  1:54833/1‑71 (MQ=255)
              ccTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGt                                                                     >  1:483494/1‑61 (MQ=255)
              ccTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGAtt                                                           >  1:3583269/1‑71 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAg                                                                      >  1:3713502/1‑59 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAg                                                                      >  1:2721518/1‑59 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAg                                                                      >  1:2443891/1‑59 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAg                                                                      >  1:422402/1‑59 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTc                                                          >  1:2759620/1‑71 (MQ=255)
                tttAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                         <  1:2022116/55‑1 (MQ=255)
                  tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCa                                                        <  1:2573488/70‑1 (MQ=255)
                         cAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTccc                                                <  1:1166397/71‑1 (MQ=255)
                         cAGACCTGACGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTccc                                                <  1:1408172/71‑1 (MQ=255)
                              cTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTccc                                                >  1:418916/1‑66 (MQ=255)
                              cTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTccc                                                >  1:3840882/1‑66 (MQ=255)
                                gCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                         >  1:1061738/1‑71 (MQ=255)
                                gCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAcct                                                                 >  1:826961/1‑46 (MQ=255)
                                gCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCAt                                                                >  1:2317905/1‑48 (MQ=255)
                                gCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCAt                                                                >  1:1720380/1‑48 (MQ=255)
                                gCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCAt                                                                >  1:2313045/1‑48 (MQ=255)
                                gCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCAt                                                                >  1:42771/1‑48 (MQ=255)
                                 ccGAAATCTGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAcc                                                                  <  1:4107075/45‑1 (MQ=255)
                                 ccGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAcc                                                                  <  1:2248518/45‑1 (MQ=255)
                                 ccGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTgg                                        >  1:75959/1‑71 (MQ=255)
                                 ccGAAATCGGCAAATCGCAACTATTTAAAGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTgg                                        >  1:3936513/1‑71 (MQ=255)
                                    aaaTCGGCAAATCTCAACTATTTAACGCGCGAGaaa                                                                        <  1:171436/36‑1 (MQ=255)
                                    aaaTCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                         <  1:1268460/35‑1 (MQ=255)
                                    aaaTCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                         <  1:2749984/35‑1 (MQ=255)
                                      aTCGGCAAATCGCAACTATTTAACGCGCGAGaaaa                                                                       >  1:1476626/1‑35 (MQ=255)
                                      aTCGGCAAATCGCAACTATTTAACGCGCGAGaaaa                                                                       >  1:1475390/1‑35 (MQ=255)
                                         ggCAGATCGCAACTATTTAACGCGCGAGAGAAGTAc                                                                   >  1:555120/1‑36 (MQ=37)
                                         ggCAAATCGCAACTATTTAACGCGCGAGAGAAGTAc                                                                   >  1:679391/1‑36 (MQ=255)
                                           cAAATCGCAACTATTTAACGCGGGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                     >  1:3143299/1‑64 (MQ=255)
                                           cAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGa                                                             <  1:636552/40‑1 (MQ=255)
                                           cAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGa                                                             <  1:1584545/40‑1 (MQ=255)
                                           cAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGa                                                             <  1:2165741/40‑1 (MQ=255)
                                           cAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGa                                                             <  1:3328761/40‑1 (MQ=255)
                                           cAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                     >  1:1517817/1‑64 (MQ=255)
                                           cAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                     >  1:304719/1‑64 (MQ=255)
                                           cAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                     >  1:2020223/1‑64 (MQ=255)
                                             aaTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                             <  1:2072864/54‑1 (MQ=255)
                                             aaTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                             <  1:3446343/54‑1 (MQ=255)
                                             aaTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                             <  1:1336229/54‑1 (MQ=255)
                                             aaTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAg                            <  1:3801319/71‑1 (MQ=255)
                                              aTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTcc                                                 >  1:1945203/1‑49 (MQ=255)
                                               tCGCATCTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCa                          <  1:2955511/71‑1 (MQ=255)
                                               tCGCAACTATTTAACGCGCGATAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCa                          <  1:2761401/71‑1 (MQ=255)
                                               tCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACat                                  <  1:3081834/63‑1 (MQ=255)
                                               tCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCa                          <  1:1251609/71‑1 (MQ=255)
                                               tCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCa                          <  1:1789939/71‑1 (MQ=255)
                                               tCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCa                          <  1:1502728/71‑1 (MQ=255)
                                                    aCTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGa                                      <  1:2082277/54‑1 (MQ=255)
                                                     cTATTTAACGCGCGAGAAAAGTACCATCGATTCCAg                                                       <  1:3236106/36‑1 (MQ=255)
                                                     cTATTTAACGCGCGAGAAAAGTACCATCGATTCCAg                                                       <  1:1922720/36‑1 (MQ=255)
                                                     cTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGc                      >  1:953991/1‑69 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCTGAt                                                     >  1:257435/1‑36 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGAt                                                     >  1:3548010/1‑36 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGAt                                                     >  1:342067/1‑36 (MQ=255)
                                                        tttAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt                  <  1:3201247/70‑1 (MQ=255)
                                                         ttAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCa                             >  1:1300179/1‑58 (MQ=255)
                                                         ttAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCa                             >  1:1318970/1‑58 (MQ=255)
                                                          tAACGCGCGAGAGAAGTACCATCGATTCCAGATGt                                                   >  1:349423/1‑35 (MQ=255)
                                                          tAACGCGCGAGAGAAGTACCATCGATTCCAGATGt                                                   >  1:718859/1‑35 (MQ=255)
                                                          tAACGCGCGAGAGAAGTACCATCGATTCCAGATGt                                                   >  1:3974881/1‑35 (MQ=255)
                                                          tAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgc               >  1:3685684/1‑71 (MQ=255)
                                                            aCGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCgcgc             <  1:1425152/71‑1 (MQ=255)
                                                                 cgAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATgg        <  1:1187363/71‑1 (MQ=255)
                                                                  gagagaAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGt       <  1:416098/71‑1 (MQ=255)
                                                                  gagagaAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGt       <  1:3318030/71‑1 (MQ=255)
                                                                  gagagaAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGt       <  1:2162145/71‑1 (MQ=255)
                                                                  gagagaAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGt       <  1:655732/71‑1 (MQ=255)
                                                                  gagaAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATc                               >  1:2373090/1‑47 (MQ=255)
                                                                  gagaAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATc                               >  1:600146/1‑47 (MQ=255)
                                                                  gagaAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATc                               >  1:3801298/1‑47 (MQ=255)
                                                                  gagaAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGt       <  1:3400839/71‑1 (MQ=255)
                                                                  gagaAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGt       <  1:922314/71‑1 (MQ=255)
                                                                   agaAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGta      <  1:3161308/71‑1 (MQ=255)
                                                                   agaAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGta      <  1:1388778/71‑1 (MQ=255)
                                                                    gaAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGtat     <  1:2688164/71‑1 (MQ=255)
                                                                     aaaaGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGtata    >  1:1033843/1‑71 (MQ=255)
                                                                     aaaaGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGT‑‑CGCAATGGTATATc  >  1:3509586/1‑71 (MQ=255)
                                                                      |                                                                       
AACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATC  >  minE/1883128‑1883269

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: