Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,883,198 A→G 50.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,1980AG50.4% 4.2 / 177.7 115F429L (TTC→CTC) rumA23S rRNA (uracil‑5)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base G (32/26);  total (60/55)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 130 total aligned reads displayed.

AACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGT  >  minE/1883128‑1883264
                                                                      |                                                                  
aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                    <  1:1096725/71‑1 (MQ=255)
aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                    <  1:2738239/71‑1 (MQ=255)
aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaa                                                                    <  1:2111209/71‑1 (MQ=255)
 aCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaaa                                                                   <  1:1759678/71‑1 (MQ=255)
 aCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaaa                                                                   <  1:2541283/71‑1 (MQ=255)
  ccATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGaaa                                                                   <  1:435983/70‑1 (MQ=255)
   cATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGagaa                                                                  <  1:287681/70‑1 (MQ=255)
   cATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAg                                                                 >  1:960709/1‑71 (MQ=255)
     tCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTa                                                               <  1:372776/71‑1 (MQ=255)
      cGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAc                                                              >  1:4102333/1‑71 (MQ=255)
      cGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAc                                                              >  1:3011788/1‑71 (MQ=255)
       gTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTAc                                                              >  1:2385400/1‑70 (MQ=255)
       gTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTAc                                                              >  1:1850570/1‑70 (MQ=255)
             tcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGAt                                                       >  1:881587/1‑71 (MQ=255)
             tcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGAt                                                       >  1:2071470/1‑71 (MQ=255)
               cTTTAGGGACTAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTc                                                     >  1:2760422/1‑71 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCa                                                            >  1:4206801/1‑64 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCa                                                            >  1:1304308/1‑64 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCa                                                            >  1:1228623/1‑64 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGAt                                                       <  1:1747893/69‑1 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGAt                                                       <  1:3379766/69‑1 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGAt                                                       <  1:1561432/69‑1 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGAt                                                       <  1:1466135/69‑1 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGAt                                                       <  1:3822309/69‑1 (MQ=255)
               cTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGAt                                                       <  1:4035607/69‑1 (MQ=255)
                tttAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGAt                                                       <  1:904737/68‑1 (MQ=255)
                  tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGGGCGAGAGAAGTAc                                                              >  1:215283/1‑59 (MQ=255)
                  tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTAc                                                              >  1:2365470/1‑59 (MQ=255)
                  tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTc                                                     >  1:3020873/1‑68 (MQ=255)
                         cAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGa                                                 >  1:3586535/1‑65 (MQ=255)
                         cAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGa                                                 >  1:1953290/1‑65 (MQ=255)
                         cAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGa                                                 >  1:3745872/1‑65 (MQ=255)
                          aGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAc                                                              >  1:1898781/1‑51 (MQ=255)
                          aGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTAc                                                              >  1:2234742/1‑51 (MQ=255)
                          aGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                          <  1:2857786/71‑1 (MQ=255)
                          aGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                          <  1:2926657/71‑1 (MQ=255)
                          aGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                          <  1:1890419/71‑1 (MQ=255)
                          aGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                          <  1:1943121/71‑1 (MQ=255)
                           gACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGagaa                                                                  >  1:3227766/1‑46 (MQ=255)
                           gACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGagaa                                                                  >  1:2759870/1‑46 (MQ=255)
                            aCCTGCCGAAATTGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                        <  1:1245877/71‑1 (MQ=255)
                            aCCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                        <  1:1467590/71‑1 (MQ=255)
                            aCCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                        <  1:2063280/71‑1 (MQ=255)
                            aCCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                        <  1:1729484/71‑1 (MQ=255)
                            aCCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtg                                        <  1:3147256/71‑1 (MQ=255)
                             ccTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGAt                                                >  1:1444942/1‑62 (MQ=255)
                                gCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                    >  1:1046224/1‑71 (MQ=255)
                                 ccGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                    <  1:2676433/70‑1 (MQ=255)
                                 ccGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                    <  1:696391/70‑1 (MQ=255)
                                 ccGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGtgtgtg                                    <  1:1935455/70‑1 (MQ=255)
                                   gAAATCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGa                                 <  1:1507478/71‑1 (MQ=255)
                                   gAAATCGGCAAATAGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTggg                                  >  1:351598/1‑70 (MQ=255)
                                    aaaTCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGt                                                                <  1:2511762/39‑1 (MQ=255)
                                    aaaTCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                <  1:3622815/71‑1 (MQ=255)
                                    aaaTCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                <  1:968648/71‑1 (MQ=255)
                                    aaaTCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                <  1:2946517/71‑1 (MQ=255)
                                     aaTCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAAc                               <  1:3119784/71‑1 (MQ=255)
                                     aaTCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAAc                               <  1:2441850/71‑1 (MQ=255)
                                      aTCGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTAc                                                              >  1:1434089/1‑39 (MQ=255)
                                        cGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAg                                                  <  1:3205699/49‑1 (MQ=255)
                                        cGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAg                                                  <  1:27633/49‑1 (MQ=255)
                                        cGGCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACat                             >  1:2874985/1‑70 (MQ=255)
                                         ggCAAATCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAg                                                  <  1:4006520/48‑1 (MQ=255)
                                          gCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                >  1:3239766/1‑65 (MQ=255)
                                           cAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATcc                         <  1:227766/71‑1 (MQ=255)
                                           cAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATcc                         <  1:348674/71‑1 (MQ=255)
                                            aaaTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATc                          >  1:3407659/1‑69 (MQ=255)
                                             aaTCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAg                       >  1:479698/1‑71 (MQ=255)
                                             aaTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCa                                                   <  1:910428/43‑1 (MQ=255)
                                             aaTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCa                                                   <  1:2048908/43‑1 (MQ=255)
                                              aTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                <  1:4004507/61‑1 (MQ=255)
                                              aTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                <  1:3345591/61‑1 (MQ=255)
                                              aTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAg                       >  1:1660430/1‑70 (MQ=255)
                                              aTCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAg                       >  1:1138343/1‑70 (MQ=255)
                                               tCGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCa                     >  1:2536194/1‑71 (MQ=255)
                                                cGCAACTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCAt                    >  1:4100126/1‑71 (MQ=255)
                                                   aaCTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGc                 <  1:3524598/71‑1 (MQ=255)
                                                    aCTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                <  1:795321/55‑1 (MQ=255)
                                                     cTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                <  1:2615068/54‑1 (MQ=255)
                                                     cTATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGc                 >  1:76360/1‑69 (MQ=255)
                                                     cTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCa                        >  1:1294399/1‑62 (MQ=255)
                                                     cTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCa                        >  1:3017133/1‑62 (MQ=255)
                                                      tATTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                <  1:3928323/53‑1 (MQ=255)
                                                       aTTTAACGCGCGAGAGAAGTACCATCGATTCCAGATGTcc                                            >  1:939601/1‑40 (MQ=255)
                                                       aTTTAACGCGCGAGAGAAGTAACATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt             >  1:3577081/1‑71 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtgtg                                      >  1:3281225/1‑46 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtgtg                                      >  1:2152146/1‑46 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCgtgtg                                      >  1:3640327/1‑46 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt             >  1:1844611/1‑71 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt             >  1:2414187/1‑71 (MQ=255)
                                                       aTTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt             >  1:3721528/1‑71 (MQ=255)
                                                           aaCGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgc          <  1:1259808/70‑1 (MQ=255)
                                                           aaCGCGCGAGAGAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgc          <  1:809938/70‑1 (MQ=255)
                                                           aaCGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg         >  1:1283976/1‑71 (MQ=255)
                                                              gcgcgAGAAAAGTACCATCGATTCCAGCTGTCCCgtgtgt                                     <  1:3813447/40‑1 (MQ=255)
                                                              gcgcgAGAAAAGTACCATCGATTCCAGATGTCCCgtgtgt                                     <  1:2723203/40‑1 (MQ=255)
                                                              gcgcgAGAAAAGTACCATCGATTCCAGATGTCCCgtgtgt                                     <  1:1123863/40‑1 (MQ=255)
                                                                gcgAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGCAACATATCCAGCATCGCCAGTCGCGCAATg    >  1:4141444/1‑71 (MQ=255)
                                                                     aaaaGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGt  <  1:2663417/68‑1 (MQ=255)
                                                                     aaaaGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGt  <  1:3820064/68‑1 (MQ=255)
                                                                      |                                                                  
AACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGT  >  minE/1883128‑1883264

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: