Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 291,204 T→C 50.0% I274T (ATC→ACC)  ppiD → peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE291,2040TC50.0% ‑5.7 / 70.5 48I274T (ATC→ACC) ppiDpeptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base T (10/14);  new base C (12/12);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTG  >  minE/291137‑291272
                                                                   |                                                                    
atCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCatcatc                                                                     >  1:162115/1‑69 (MQ=255)
 tCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAg                                                                  <  1:118701/71‑1 (MQ=255)
 tCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCACCCa                                                                   <  1:2971418/70‑1 (MQ=255)
           ctacGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCACCCAGa                                                                 >  1:3122051/1‑62 (MQ=255)
           ctacGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCACCCAGa                                                                 >  1:3140693/1‑62 (MQ=255)
           ctacGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACt                                                         >  1:3104991/1‑70 (MQ=255)
            tacGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGAcc                                                               >  1:1850973/1‑63 (MQ=255)
            tacGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCaa                                                             >  1:1496661/1‑65 (MQ=255)
               gACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTg                                                        >  1:1097170/1‑67 (MQ=255)
               gACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTg                                                        <  1:1417619/67‑1 (MQ=255)
                  cAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCAccc                                                                    <  1:3651419/52‑1 (MQ=255)
                   aGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAg                                                                  <  1:1315739/53‑1 (MQ=255)
                   aGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAg                                                                  <  1:724861/53‑1 (MQ=255)
                   aGCATCAGGATAAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATg                                                  >  1:1784397/1‑69 (MQ=255)
                    gCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAg                                               >  1:2526810/1‑71 (MQ=255)
                     cATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGc                                              >  1:905566/1‑71 (MQ=255)
                           gATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGc                                        <  1:2493395/71‑1 (MQ=255)
                           gATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGc                                        <  1:2657966/71‑1 (MQ=255)
                           gATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGc                                        <  1:580998/71‑1 (MQ=255)
                            aTCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCg                                       <  1:214953/71‑1 (MQ=255)
                              cAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCgg                                      >  1:1580820/1‑70 (MQ=255)
                                      ccAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGt                                     <  1:2688169/63‑1 (MQ=255)
                                      ccAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGt                                     <  1:2113661/63‑1 (MQ=255)
                                      ccAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGt                                     <  1:2311941/63‑1 (MQ=255)
                                      ccAGCCGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCgg                                      >  1:3248652/1‑62 (MQ=255)
                                      ccAGCCGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCgg                                      >  1:3355430/1‑62 (MQ=255)
                                      ccAGCCGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCgg                                      >  1:2241791/1‑62 (MQ=255)
                                      ccAGCCGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCgg                                      >  1:2247751/1‑62 (MQ=255)
                                       cAGCCGCAGCGTACCCACTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGt                                     <  1:264108/62‑1 (MQ=255)
                                        aGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCg                                       >  1:1510301/1‑59 (MQ=255)
                                        aGCCGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGa                                                 >  1:919965/1‑49 (MQ=255)
                                        aGCCGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAg                           >  1:1306146/1‑71 (MQ=255)
                                          ccGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCt                         <  1:465310/71‑1 (MQ=255)
                                          ccGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGc                          <  1:3538670/70‑1 (MQ=255)
                                           cGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTg                        <  1:704939/71‑1 (MQ=255)
                                           cGCAGCGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTg                        <  1:1649511/71‑1 (MQ=255)
                                             cAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGaa                      >  1:2163287/1‑71 (MQ=255)
                                             cAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGaa                      >  1:3035796/1‑71 (MQ=255)
                                               gCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAAt                     >  1:2370012/1‑70 (MQ=255)
                                                cGTACCCGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATaa                   <  1:1716466/71‑1 (MQ=255)
                                                  tACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAg                 >  1:2005572/1‑71 (MQ=255)
                                                    cccGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAgg                >  1:145010/1‑70 (MQ=255)
                                                     ccGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGt                                     <  1:144999/48‑1 (MQ=255)
                                                     ccGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCg              <  1:2110621/71‑1 (MQ=255)
                                                     ccGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCg              <  1:3437545/71‑1 (MQ=255)
                                                     ccGCTACAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCg              <  1:416443/71‑1 (MQ=255)
                                                           cAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCg              <  1:3478859/65‑1 (MQ=255)
                                                           cAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGAtt        <  1:1673759/71‑1 (MQ=255)
                                                           cAGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGAtt        <  1:681586/71‑1 (MQ=255)
                                                            aGCATCACCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGAttt       <  1:3653721/71‑1 (MQ=255)
                                                              catcatCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTtg     <  1:3392292/71‑1 (MQ=255)
                                                              catcacCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGAtttt      <  1:2119150/70‑1 (MQ=255)
                                                              catcacCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTtg     <  1:1522133/71‑1 (MQ=255)
                                                              catcacCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTtg     <  1:3120611/71‑1 (MQ=255)
                                                                tcaccCAGACCAAAACTGAAGCTGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTtgct   >  1:1167343/1‑71 (MQ=255)
                                                                tcaccCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTtgct   >  1:1705182/1‑71 (MQ=255)
                                                                  atcCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTtgctg  <  1:1871853/70‑1 (MQ=255)
                                                                   |                                                                    
ATCCAGAGCTACTACGACCAGCATCAGGATCAATTCACCCAGCCGCAGCGTACCCGCTACAGCATCATCCAGACCAAAACTGAAGATGAAGCGAAAGCGGTACTTGATGAGCTGAATAAAGGCGGTGATTTTGCTG  >  minE/291137‑291272

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: