Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,928,279 C→T 48.5% D226D (GAC→GAT yjgP → conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,928,2790CT48.5% 7.8 / 126.8 140D226D (GAC→GATyjgPconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (43/29);  new base T (51/17);  total (94/46)
Fisher's exact test for biased strand distribution p-value = 7.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 140 total aligned reads displayed.

CGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCG  >  minE/2928213‑2928340
                                                                  |                                                             
cgtcACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACt                                                              >  1:2293910/1‑68 (MQ=255)
cgtcACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACt                                                              >  1:1282367/1‑68 (MQ=255)
cgtcACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACt                                                              >  1:2448523/1‑68 (MQ=255)
    aCTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAg                                                         <  1:1747245/69‑1 (MQ=255)
    aCTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAg                                                         <  1:3384950/69‑1 (MQ=255)
       ctcAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGAt                                                      <  1:3145898/69‑1 (MQ=255)
       ctcAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGAt                                                      <  1:2844681/69‑1 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCa                                                          >  1:2376978/1‑60 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTAt                                                   >  1:1020749/1‑67 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATc                                                  >  1:1335076/1‑68 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATc                                                  >  1:770228/1‑68 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:2628036/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:2782598/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:296372/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:2158023/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:1897049/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:687140/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:1817771/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:3026762/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:1028983/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:75219/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:124667/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCa                                                 >  1:1305856/1‑69 (MQ=255)
            ccAGGGAACGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTAGAGGATTATCa                                                 >  1:2811875/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCa                                                 >  1:3387847/1‑63 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGAt                                           >  1:2943134/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGAt                                           >  1:1072498/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGAt                                           >  1:3252975/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGAt                                           >  1:1771553/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGAt                                           >  1:752144/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGAt                                           >  1:1430282/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCTAt                                           >  1:2367003/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGa                                            >  1:2433104/1‑68 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:523670/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2717678/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2648081/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2729541/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:173895/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:150207/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:3294262/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2750275/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2457912/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:1388902/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2409671/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:735904/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:1841482/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:594761/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2201298/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2199766/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:1162131/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:1688011/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2061272/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2045357/1‑69 (MQ=255)
                  aaCGCGCTTCGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:1366422/1‑69 (MQ=255)
                  aaCGCGCTTAGAAGGCACTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGAt                                           >  1:2662412/1‑69 (MQ=255)
                   aCGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATc                                          <  1:1234938/69‑1 (MQ=255)
                   aCGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATc                                          <  1:1082185/69‑1 (MQ=255)
                       gcTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTg                                      <  1:839/69‑1 (MQ=255)
                       gcTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATAAGGCGATCATTg                                      <  1:621075/69‑1 (MQ=255)
                         ttCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGt                                    <  1:627744/69‑1 (MQ=255)
                          tCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATGATTGGt                                    <  1:1202348/68‑1 (MQ=255)
                                gCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAgg                             >  1:2644571/1‑69 (MQ=255)
                                gCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAgg                             >  1:673194/1‑69 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:2604952/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:3229095/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:869248/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:424736/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:1171751/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:2204523/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:2100801/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:2079061/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:1871463/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:2761111/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCgg                          <  1:1345777/69‑1 (MQ=255)
                                   cTGCACTGTTACGTGATTTCCGCATTACGGATTTCCAGGAATATCAGGCGATCATTGGTCACCAGGCgg                          <  1:3330668/69‑1 (MQ=255)
                                     gCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGt                         >  1:703303/1‑68 (MQ=255)
                                     gCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTg                        >  1:2608995/1‑69 (MQ=255)
                                     gCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTg                        >  1:1016396/1‑69 (MQ=255)
                                     gCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTg                        >  1:3317875/1‑69 (MQ=255)
                                     gCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTg                        >  1:2576425/1‑69 (MQ=255)
                                     gCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTg                        >  1:2487879/1‑69 (MQ=255)
                                             aCGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGa                >  1:3096613/1‑69 (MQ=255)
                                                    ttCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAAc         <  1:69135/69‑1 (MQ=255)
                                                    ttCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAAc         <  1:3058320/69‑1 (MQ=255)
                                                     tCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAAc         <  1:1594638/68‑1 (MQ=255)
                                                         cATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATa     <  1:2382283/68‑1 (MQ=255)
                                                         cATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATa     <  1:2980005/68‑1 (MQ=255)
                                                         cATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATa     <  1:901703/68‑1 (MQ=255)
                                                         cATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATa     <  1:936167/68‑1 (MQ=255)
                                                         cATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATa     <  1:973465/68‑1 (MQ=255)
                                                          aTTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATa     <  1:3158479/67‑1 (MQ=255)
                                                          aTTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATAcc   <  1:2278406/69‑1 (MQ=255)
                                                          aTTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATAcc   <  1:3324430/69‑1 (MQ=255)
                                                          aTTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATAcc   <  1:1822260/69‑1 (MQ=255)
                                                             aCGGACTTCCAGGATTATCAGGCGATCATTGGTCACCa                               >  1:1721201/1‑38 (MQ=255)
                                                             aCGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCg  >  1:665289/1‑67 (MQ=255)
                                                             aCGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCg  >  1:2408184/1‑67 (MQ=255)
                                                             aCGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCg  >  1:748457/1‑67 (MQ=255)
                                                             aCGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCg  >  1:1538358/1‑67 (MQ=255)
                                                             aCGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCg  >  1:2989125/1‑67 (MQ=255)
                                                                  |                                                             
CGTCACTCTCAACCAGGGAACGCGCTTCGAAGGCACTGCATTGTTACGTGATTTCCGCATTACGGACTTCCAGGATTATCAGGCGATCATTGGTCACCAGGCGGTGGCGCTCGACCCGAACGATACCG  >  minE/2928213‑2928340

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: