Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,325,765 A→T 27.1% I107I (ATT→ATA cpxP ← periplasmic protein combats stress

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,325,7650AT27.1% 113.4 / 46.2 81I107I (ATT→ATAcpxPperiplasmic protein combats stress
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base T (9/13);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GCTGCTCCGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGT  >  minE/2325701‑2325830
                                                                |                                                                 
gctgctCCGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAAttt                                                               >  1:749216/1‑69 (MQ=255)
gctgctCCGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAAttt                                                               >  1:2220879/1‑69 (MQ=255)
gctgctCCGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAAttt                                                               >  1:2448932/1‑69 (MQ=255)
gctgctCCGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAAttt                                                               >  1:2927368/1‑69 (MQ=255)
gctgctCCGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAAttt                                                               >  1:30385/1‑69 (MQ=255)
   gctCCGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCt                                                            <  1:2616825/69‑1 (MQ=255)
      ccGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:1305876/68‑1 (MQ=255)
      ccGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:3117549/68‑1 (MQ=255)
      ccGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:2779153/68‑1 (MQ=255)
      ccGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:2128146/68‑1 (MQ=255)
             tAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCAt                                                  <  1:1258858/69‑1 (MQ=255)
                 aGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCAttttt                                              <  1:764056/69‑1 (MQ=255)
                 aGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCAttttt                                              <  1:2314530/69‑1 (MQ=255)
                 aGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCAttttt                                              <  1:529809/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:1271196/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:419700/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:2385007/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:2467334/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:2070632/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:2543732/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:696255/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:1875422/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:738978/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:1739750/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:3046847/69‑1 (MQ=255)
                    cGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCt                                                            <  1:1755378/52‑1 (MQ=255)
                     gATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:654742/68‑1 (MQ=255)
                        aCATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCt                                       >  1:1150073/1‑69 (MQ=255)
                        aCATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCt                                       >  1:314911/1‑69 (MQ=255)
                        aCATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCt                                       >  1:570781/1‑69 (MQ=255)
                        aCATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCt                                       >  1:1998395/1‑69 (MQ=255)
                        aCATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCt                                       >  1:756101/1‑69 (MQ=255)
                        aCATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATGTTTTCTGCCt                                       >  1:1305999/1‑69 (MQ=255)
                         cATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCtg                                      >  1:2420977/1‑69 (MQ=255)
                         cATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCtg                                      >  1:315410/1‑69 (MQ=255)
                         cATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCtg                                      >  1:2110077/1‑69 (MQ=255)
                         cATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCtg                                      >  1:426294/1‑69 (MQ=255)
                         cATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCtg                                      >  1:1471540/1‑69 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:3262454/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:3257885/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          >  1:317610/1‑46 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          >  1:862626/1‑46 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:3301342/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:566771/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:3366266/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          >  1:2202244/1‑46 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:2969596/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          >  1:1332993/1‑46 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:3379487/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          >  1:1794863/1‑46 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:266016/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          >  1:710509/1‑46 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:2085484/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:243543/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          <  1:692564/46‑1 (MQ=255)
                            ttGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCa                                                          >  1:424385/1‑46 (MQ=255)
                                ttGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCTGTgcgcg                                >  1:3000473/1‑68 (MQ=255)
                                ttGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCTGTgcgcg                                >  1:821697/1‑68 (MQ=255)
                                ttGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCTGTgcgcg                                >  1:1244271/1‑68 (MQ=255)
                                ttGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCTGTgcgcg                                >  1:2456748/1‑68 (MQ=255)
                                ttGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCTGTgcgcg                                >  1:2369391/1‑68 (MQ=255)
                                      aCTTTTGCCATCTCAACCTGACGAGCTATTTGCTCATTCGCCATTTTTTCt                                           <  1:1290450/51‑1 (MQ=255)
                                                  tcAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCaaaa             >  1:3341586/1‑69 (MQ=255)
                                                  tcAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCaaa              >  1:2719429/1‑68 (MQ=255)
                                                           cGAGCTATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCg    >  1:2706590/1‑69 (MQ=255)
                                                           cGAGCTATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCg    >  1:1176144/1‑69 (MQ=255)
                                                           cGAGCTATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCg    >  1:3092081/1‑69 (MQ=255)
                                                           cGAGCTATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCg    >  1:1745466/1‑69 (MQ=255)
                                                           cGAGCTATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCg    >  1:185703/1‑69 (MQ=255)
                                                           cGAGCTATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCg    >  1:1089114/1‑69 (MQ=255)
                                                           cGAGCTATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCg    >  1:1914511/1‑69 (MQ=255)
                                                           cGAGCTATTTGCTCATTCGCCATTTTTTC‑‑CCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGt  >  1:2186680/1‑69 (MQ=255)
                                                           cGAGCTATTTGCTCATTCGCCATATTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCg    >  1:3300445/1‑69 (MQ=255)
                                                             aGCAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGt  <  1:1485589/69‑1 (MQ=255)
                                                               cAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGt  <  1:683431/67‑1 (MQ=255)
                                                               cAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGt  <  1:603343/67‑1 (MQ=255)
                                                               cAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGt  <  1:1838565/67‑1 (MQ=255)
                                                               cAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGt  <  1:1207072/67‑1 (MQ=255)
                                                               cAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGt  <  1:3380205/67‑1 (MQ=255)
                                                               cAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGt  <  1:993829/67‑1 (MQ=255)
                                                               cAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTAATCAAAATTTTCTGCGGt  <  1:1320926/67‑1 (MQ=255)
                                                                aaTTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGt  <  1:2951563/66‑1 (MQ=255)
                                                                |                                                                 
GCTGCTCCGGCGTTAACAGGCGATACATTTGGTTGCGGACTTTTGCCATCTCAACCTGACGAGCAATTTGCTCATTCGCCATTTTTTCTGCCTGTGCGCGCACAGCGTTTTCATCAAAATTTTCTGCGGT  >  minE/2325701‑2325830

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: