Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,613,998 A→G 48.7% I126I (ATT→ATC eutH ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,613,9980AG48.7% ‑5.8 / 31.3 37I126I (ATT→ATCeutHpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (17/2);  new base G (8/10);  total (25/12)
Fisher's exact test for biased strand distribution p-value = 5.14e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TCGATAATGCCGAGCGCCACCGGAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCGAGGAAGAAGCCGCC  >  minE/1613936‑1614066
                                                              |                                                                    
tCGATAATGCCGAGCGCCACCGGAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTaa                                                                 >  1:73989/1‑68 (MQ=255)
tCGATAATGCCGAGCGCCACCGGAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTaa                                                                 >  1:1185699/1‑68 (MQ=255)
tCGATAATGCCGAGCGCCACCGGAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTaa                                                                 >  1:3251834/1‑68 (MQ=255)
tCGATAATGCCGAGCGCCACCGGAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTaa                                                                 >  1:402642/1‑68 (MQ=255)
 cGATAATGCCGAGCGCCACCGGAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTaa                                                                 >  1:650463/1‑67 (MQ=255)
 cGATAATGCCGAGCGCCACCGGAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTaa                                                                 >  1:2253997/1‑67 (MQ=255)
         ccGAGCGCCACCGGAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTAACCCAGAGtat                                                       >  1:3203276/1‑69 (MQ=255)
                   ctggAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCAcgct                                              >  1:2163460/3‑67 (MQ=255)
                   ctggAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCAcgct                                              >  1:2671353/3‑67 (MQ=255)
                   ctggAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCAcgct                                              >  1:2947008/3‑67 (MQ=255)
                   ctggAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCAcgct                                              >  1:3249737/3‑67 (MQ=255)
                   ctggAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCAcgct                                              >  1:2249370/3‑67 (MQ=255)
                   ctggAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCAcgct                                              >  1:2229887/3‑67 (MQ=255)
                   ctggAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCAcgct                                              >  1:3383850/3‑67 (MQ=255)
                   ctggAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCAcgct                                              >  1:3387428/3‑67 (MQ=255)
                           gAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTcg                                     >  1:459975/1‑69 (MQ=255)
                           gAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTcg                                     >  1:1294065/1‑69 (MQ=255)
                           gAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTcg                                     >  1:285834/1‑69 (MQ=255)
                           gAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTc                                      >  1:2157967/1‑68 (MQ=255)
                           gAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTc                                      >  1:1180877/1‑68 (MQ=255)
                                       gTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCACGCTGCTACGTCACCGCCCGCCAGc                         <  1:2518671/69‑1 (MQ=255)
                                       gTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCACGCTGCTACGTCACCGCCCGCCAGc                         <  1:3247162/69‑1 (MQ=255)
                                       gTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCACGCTGCTACGTCACCGCCCGCCAGc                         <  1:558411/69‑1 (MQ=255)
                                       gTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCACGCTGCTACGTCACCGCCCGCCAGc                         <  1:3073621/69‑1 (MQ=255)
                                       gTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCACGCTGCTACGTCACCGCCCGCCAGc                         <  1:571910/69‑1 (MQ=255)
                                       gTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCACGCTGCTACGTCACCGCCCGCCAGc                         <  1:2174726/69‑1 (MQ=255)
                                       gTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCACGCTGCTACGTCACCGCCCGCCAGc                         <  1:1040811/69‑1 (MQ=255)
                                       gTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCACGCTGCTACGTCACCGCCCGCCAGc                         <  1:1209605/69‑1 (MQ=255)
                                       gTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCACGCTGCTACGTCACCGCCCGCCAGc                         <  1:2877644/69‑1 (MQ=255)
                                       gTTGGCCCCATCATCGACCCGAGGATTAACCCAGAGTATAGCCACGCTGCTACGTCACCGCCCGCCAGc                         <  1:2558550/69‑1 (MQ=255)
                                                 tcatcGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCg                >  1:2077449/1‑68 (MQ=255)
                                                   atcGACCCGAGAATTAACCCAGAGTATAGCCACGCGGc                                            <  1:1077037/38‑1 (MQ=255)
                                                   atcGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGctc                       >  1:2307923/1‑59 (MQ=255)
                                                   atcGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCGa               >  1:2425823/1‑67 (MQ=255)
                                                   atcGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCGa               >  1:2379163/1‑67 (MQ=255)
                                                   atcGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCAGCCAGCTCTTTGGCGa               >  1:1021468/1‑67 (MQ=255)
                                                   atcCACCCGAGAATTAACCCAGAGTATAGCCACGCGGc                                            <  1:3367693/38‑1 (MQ=255)
                                                              aaTTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCGAGGAAGAAgccgcc  <  1:2315247/69‑1 (MQ=255)
                                                              aaTTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCGAGGAAGAAgccgcc  <  1:309726/69‑1 (MQ=255)
                                                              aaTTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCGAGGAAGAAgccgcc  <  1:1958711/69‑1 (MQ=255)
                                                              aaTTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCGAGGAAGAAgccgcc  <  1:1846565/69‑1 (MQ=255)
                                                              aaTTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCGAGGAAGAAgccgcc  <  1:1285021/69‑1 (MQ=255)
                                                              aaTTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCGAGGAAGAAgccgcc  <  1:675968/69‑1 (MQ=255)
                                                              |                                                                    
TCGATAATGCCGAGCGCCACCGGAATGGAAAACACAATCGTTGGCCCCATCATCGACCCGAGAATTAACCCAGAGTATAGCCACGCGGCTACGTCGCCGCCCGCCAGCTCTTTGGCGAGGAAGAAGCCGCC  >  minE/1613936‑1614066

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: