Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 572,219 A→G 45.7% I176I (ATT→ATC aqpZ ← aquaporin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE572,2190AG45.7% ‑0.6 / 27.0 35I176I (ATT→ATCaqpZaquaporin
Reads supporting (aligned to +/- strand):  ref base A (4/15);  new base G (13/3);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 6.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTG  >  minE/572155‑572285
                                                                |                                                                  
ccTGGAAGATAGCAACCGCGGTGCTGCGCGCCTGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTa                                                                <  1:1079448/69‑1 (MQ=255)
ccTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTa                                                                <  1:3280701/69‑1 (MQ=255)
ccTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTa                                                                <  1:2343845/69‑1 (MQ=255)
              aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTa                                                                <  1:523685/55‑1 (MQ=255)
              aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:1134956/69‑1 (MQ=255)
              aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:657155/69‑1 (MQ=255)
              aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:2321927/69‑1 (MQ=255)
              aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:2249836/69‑1 (MQ=255)
              aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:1749768/69‑1 (MQ=255)
              aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:1565933/69‑1 (MQ=255)
              aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:155337/69‑1 (MQ=255)
               ccGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTAACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:1932223/68‑1 (MQ=255)
                                  gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt                              <  1:1157276/69‑1 (MQ=255)
                                  gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt                              <  1:1428513/69‑1 (MQ=255)
                                  gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt                              <  1:2581928/69‑1 (MQ=255)
                                   ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGAt                              >  1:3208379/1‑68 (MQ=255)
                                   ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                             >  1:333286/1‑69 (MQ=255)
                                              gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCAGCGCAAAAcc                  >  1:1980053/1‑69 (MQ=255)
                                                          aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGaa      <  1:1948636/69‑1 (MQ=255)
                                                          aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGaa      <  1:1435480/69‑1 (MQ=255)
                                                          aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGaa      <  1:2623768/69‑1 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa      >  1:948065/1‑68 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa      >  1:3009856/1‑68 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa      >  1:2063640/1‑68 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa      >  1:1992461/1‑68 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:2804589/1‑69 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:2518255/1‑69 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:1209389/1‑69 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:1186120/1‑69 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:3360277/1‑69 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:444566/1‑69 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:1178516/1‑69 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:857902/1‑69 (MQ=255)
                                                           aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:1480765/1‑69 (MQ=255)
                                                              cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTg  >  1:2918117/1‑69 (MQ=255)
                                                                |                                                                  
CCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTG  >  minE/572155‑572285

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: