Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 51,445 A→G 33.3% A312A (GCT→GCC)  surA ← peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE51,4450AG33.3% 77.4 / 56.3 102A312A (GCT→GCC) surApeptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base A (41/27);  new base G (26/8);  total (67/35)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 104 total aligned reads displayed.

TCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGAA  >  minE/51389‑51512
                                                        |                                                                   
tCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGt                                                         >  1:793695/1‑69 (MQ=255)
       agaACTCTTTCGCTGCTGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGt                                                   <  1:2660545/68‑1 (MQ=255)
       agaACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGt                                                   <  1:2558354/68‑1 (MQ=255)
       agaACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGt                                                   <  1:2725241/68‑1 (MQ=255)
       agaACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGt                                                   <  1:172473/68‑1 (MQ=255)
       agaACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGt                                                   <  1:1175350/68‑1 (MQ=255)
       agaACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGt                                                   <  1:572602/68‑1 (MQ=255)
       agaACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGt                                                   <  1:841822/68‑1 (MQ=255)
       agaACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTAACACTCTTGATATCGGCAGCAATCTGTTCCAGt                                                   <  1:560554/68‑1 (MQ=255)
         aaCTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGttt                                                 >  1:881037/1‑68 (MQ=255)
         aaCTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGttt                                                 >  1:3195293/1‑68 (MQ=255)
           ctctTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTcac                                              >  1:2189258/1‑69 (MQ=255)
           ctctTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTcac                                              >  1:3210793/1‑69 (MQ=255)
           ctctTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTcac                                              >  1:1419066/1‑69 (MQ=255)
           ctctTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTcac                                              >  1:1644716/1‑69 (MQ=255)
           ctctTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTcac                                              >  1:2061721/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTcacgc                                            >  1:3243138/1‑66 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTcacgc                                            >  1:2079854/1‑66 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTcacgc                                            >  1:1675245/1‑66 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTcacgc                                            >  1:3017247/1‑66 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTcacg                                             >  1:1039619/1‑66 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:2740905/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:3028365/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:3036567/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:335465/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:445519/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:2538185/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:2351464/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:2286961/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:710511/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:2086/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:948013/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:2065410/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:2496204/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:181808/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:1766230/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:1256070/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:1450362/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:1527623/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGTTCCAGTTTCAcgcg                                           >  1:1928396/1‑69 (MQ=255)
              tttCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCGGCAGCAATCTGCTCCAGTTTCAcgcg                                           >  1:1028152/1‑69 (MQ=255)
                     gCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGt                                    <  1:2668055/69‑1 (MQ=255)
                     gCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGt                                    <  1:469559/69‑1 (MQ=255)
                             aaaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTg                                     >  1:134705/1‑60 (MQ=255)
                             aaaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagt                            >  1:2001192/1‑69 (MQ=255)
                             aaaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagt                            >  1:1240267/1‑69 (MQ=255)
                              aaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagtc                           >  1:372980/1‑69 (MQ=255)
                              aaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagtc                           >  1:2718441/1‑69 (MQ=255)
                              aaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagtc                           >  1:3033947/1‑69 (MQ=255)
                              aaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagtc                           >  1:1168494/1‑69 (MQ=255)
                              aaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagtc                           >  1:826415/1‑69 (MQ=255)
                              aaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagtc                           >  1:832456/1‑69 (MQ=255)
                              aaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagtc                           >  1:981230/1‑69 (MQ=255)
                              aaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagtc                           >  1:1960550/1‑69 (MQ=255)
                              aaaGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCgtcagt                            >  1:2541246/1‑68 (MQ=255)
                                 gTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATg                        <  1:2202125/69‑1 (MQ=255)
                                 gTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATg                        <  1:2055420/69‑1 (MQ=255)
                                 gTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATg                        <  1:1154370/69‑1 (MQ=255)
                                 gTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATg                        <  1:115449/69‑1 (MQ=255)
                                 gTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATg                        <  1:823118/69‑1 (MQ=255)
                                 gTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATg                        <  1:793057/69‑1 (MQ=255)
                                 gTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATg                        <  1:2556682/69‑1 (MQ=255)
                                 gTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATg                        <  1:1720828/69‑1 (MQ=255)
                                 gTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATg                        <  1:2464781/69‑1 (MQ=255)
                                  tCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATg                        <  1:279062/68‑1 (MQ=255)
                                           cACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGAc               <  1:2126005/68‑1 (MQ=255)
                                           cACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGAc               <  1:1423438/68‑1 (MQ=255)
                                           cACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGAc               <  1:2961366/68‑1 (MQ=255)
                                           cACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGAc               <  1:671415/68‑1 (MQ=255)
                                           cACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGAc               <  1:1310678/68‑1 (MQ=255)
                                           cACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGAc               <  1:501878/68‑1 (MQ=255)
                                             ctctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            <  1:236155/69‑1 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:1575202/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:2994210/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:763832/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:1162619/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:2313675/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:522293/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:867613/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:927508/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:1043655/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:1327822/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:2281619/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGt            >  1:1370146/1‑67 (MQ=255)
                                               ctTGATATCAGCAGCAATCTGTTCCAGTTTAACACGGTCCTGTTCGTCAGTCATGATCGGCGACgg             >  1:2502884/1‑66 (MQ=255)
                                                ttGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTc         >  1:2259623/1‑69 (MQ=255)
                                                ttGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTc         >  1:2136262/1‑69 (MQ=255)
                                                ttGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTc         >  1:2064721/1‑69 (MQ=255)
                                                ttGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTc         >  1:2680194/1‑69 (MQ=255)
                                                ttGATATCAGAAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTc         >  1:2129616/1‑69 (MQ=255)
                                                ttGATATCAACAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTc         >  1:945619/1‑69 (MQ=255)
                                                     atCAGCAGCAATCTGTTCCAGTTTCACAGGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTcagcag    <  1:3301910/69‑1 (MQ=255)
                                                     atCAGCAGCAATCTGTTCCAGTTTCACACGGGGCTGTTCGTCAGTCATGATCGGCGACGGTTTcagcag    <  1:1573129/69‑1 (MQ=255)
                                                     atCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTcagcag    <  1:1155402/69‑1 (MQ=255)
                                                     atCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTcagcag    <  1:2045302/69‑1 (MQ=255)
                                                     atCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTcagcag    <  1:738747/69‑1 (MQ=255)
                                                      tCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTcagcag    >  1:489895/1‑68 (MQ=255)
                                                      tCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTcagcag    <  1:1978674/68‑1 (MQ=255)
                                                        agcagcAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGaa  <  1:785648/68‑1 (MQ=255)
                                                        agcagcAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGaa  <  1:3338375/68‑1 (MQ=255)
                                                        |                                                                   
TCCTGAGAGAACTCTTTCGCTGCGGCAGCAAAAGTCGTTTTACCACTCTTGATATCAGCAGCAATCTGTTCCAGTTTCACACGGGCCTGTTCGTCAGTCATGATCGGCGACGGTTTCAGCAGAA  >  minE/51389‑51512

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: ‑