Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,357,334 T→G 30.2% T596T (ACA→ACC yehB ← predicted outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,357,3340TG30.2% 63.7 / 41.1 63T596T (ACA→ACCyehBpredicted outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (35/9);  new base G (10/9);  total (45/18)
Fisher's exact test for biased strand distribution p-value = 3.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TTCATTTCCCTGATGTTGATGACTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACATA  >  minE/1357271‑1357400
                                                               |                                                                  
ttCATTTCCCTGATACTGATAACTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTc                                                                <  1:288081/68‑1 (MQ=255)
ttCATTTCCCTGATACTGATAACTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTc                                                                <  1:2359727/68‑1 (MQ=255)
ttCATTTCCCTGATACTGATAACTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTc                                                                <  1:2590596/68‑1 (MQ=255)
ttCATTTCCCTGATACTGATAACTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTc                                                                <  1:1815913/68‑1 (MQ=255)
ttCATTTCCCTGATACTGATAACTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTc                                                                <  1:1693352/68‑1 (MQ=255)
ttCATTTCCCTGATACTGATAACTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTc                                                                <  1:3241475/68‑1 (MQ=255)
ttCATTTCCCTGATACTGATAACTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTc                                                                <  1:3132825/68‑1 (MQ=255)
ttCATTTCCCTGATACTGATAACTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTc                                                                <  1:1253966/68‑1 (MQ=255)
ttCATTTCCCTGATACTGATAACTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTc                                                                <  1:3219998/68‑1 (MQ=255)
     ttCCCTGATGTTGATGACTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGAt                                                           <  1:2454205/68‑1 (MQ=255)
         cTGATGTTGATGACTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGAt                                                           <  1:2308505/64‑1 (MQ=255)
                 gatgaCTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGtatta                                               >  1:902569/1‑68 (MQ=255)
                 gatgaCTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGtatta                                               >  1:1990382/1‑68 (MQ=255)
                 gatgaCTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGtatta                                               >  1:512350/1‑68 (MQ=255)
                 gatgaCTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGtat                                                 <  1:1169715/66‑1 (MQ=255)
                 gatgaCTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTACTGATAATCCCGtattat                                              >  1:323929/1‑69 (MQ=255)
                 gatgaCTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTACTGATAATCCCGtattat                                              >  1:319813/1‑69 (MQ=255)
                     aCTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGa                                          >  1:1579032/1‑69 (MQ=255)
                     aCTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGa                                          >  1:1478538/1‑69 (MQ=255)
                     aCTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGa                                          >  1:393037/1‑69 (MQ=255)
                      cTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTCCTGa                                                            >  1:1550407/1‑50 (MQ=255)
                      cTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTCCTGa                                                            >  1:3185353/1‑50 (MQ=255)
                      cTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTCCTGa                                                            >  1:2634376/1‑50 (MQ=255)
                      cTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTCCTGa                                                            >  1:1059712/1‑50 (MQ=255)
                      cTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTCCTGa                                                            >  1:492450/1‑50 (MQ=255)
                      cTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTCCTGa                                                            >  1:302800/1‑50 (MQ=255)
                      cTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTCCTGa                                                            >  1:398085/1‑50 (MQ=255)
                      cTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTCCTGa                                                            >  1:2170915/1‑50 (MQ=255)
                      cTCAGGTTGACCCCATAGTTAAACTGATCCCGGCTTCCAACGGTTCCTGa                                                            >  1:1409225/1‑50 (MQ=255)
                      cTCAGGTTGACCCCATAGTTAAACTGATCCCAGCTTCCAACGGTTCCTGa                                                            >  1:1003170/1‑50 (MQ=255)
                           gTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCa                                      <  1:701398/67‑1 (MQ=255)
                           gTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCa                                      <  1:471382/67‑1 (MQ=255)
                           gTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCa                                      <  1:367024/67‑1 (MQ=255)
                           gTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCa                                      <  1:1352351/67‑1 (MQ=255)
                           gTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCa                                      <  1:2961393/67‑1 (MQ=255)
                              gACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAAccc                                 <  1:1328473/69‑1 (MQ=255)
                                    aTAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGatc                           >  1:2083462/1‑69 (MQ=255)
                                    aTAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGatc                           >  1:31258/1‑69 (MQ=255)
                                    aTAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGatc                           >  1:537410/1‑69 (MQ=255)
                                    aTAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGatc                           >  1:2201617/1‑69 (MQ=255)
                                       aTTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGatcat                         >  1:1074891/1‑68 (MQ=255)
                                                         cccTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTa        >  1:640107/1‑67 (MQ=255)
                                                         cccTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTa        >  1:1422133/1‑67 (MQ=255)
                                                         cccTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTa        >  1:3245179/1‑67 (MQ=255)
                                                         cccTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTa        >  1:2762344/1‑67 (MQ=255)
                                                         cccTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTa        >  1:1806605/1‑67 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACata  >  1:1005054/1‑69 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACata  >  1:1149555/1‑69 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACat   >  1:3082103/1‑68 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACat   >  1:2855769/1‑68 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACat   >  1:3311876/1‑68 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACat   >  1:2100623/1‑68 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACat   >  1:1039415/1‑68 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACa    >  1:2966391/1‑67 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACa    >  1:1623427/1‑67 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACa    >  1:1573178/1‑67 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACa    >  1:1480487/1‑67 (MQ=255)
                                                             aCTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACa    >  1:862111/1‑67 (MQ=255)
                                                              cTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGAtcatca                       >  1:3189053/1‑47 (MQ=255)
                                                              cTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGAtcatca                       >  1:2523020/1‑47 (MQ=255)
                                                              cTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGAtcatca                       >  1:1826524/1‑47 (MQ=255)
                                                              cTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACata  >  1:2616201/1‑68 (MQ=255)
                                                              cTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACata  >  1:2082254/1‑68 (MQ=255)
                                                               |                                                                  
TTCATTTCCCTGATGTTGATGACTCAGGTTGACACCATAATTGAACTGATCCCGACTCCCTACTGTTCCTGATAATCCCGTATTATTTGAGGCAAACCCCTGATCATCAAACGTCGTTGAGTTAGACATA  >  minE/1357271‑1357400

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: