Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,506,580 C→T 54.8% E253E (GAG→GAA yfbS ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,506,5800CT54.8% 20.7 / 104.7 124E253E (GAG→GAAyfbSpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (37/19);  new base T (57/11);  total (94/30)
Fisher's exact test for biased strand distribution p-value = 3.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 124 total aligned reads displayed.

AACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACG  >  minE/1506514‑1506645
                                                                  |                                                                 
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGGTGACGTTAGCGCCATAACGCTc                                                                 >  1:417315/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATACCGCTc                                                                 >  1:2186837/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:1306921/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:1740732/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:1833306/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:1490620/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:1847352/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:1984645/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:1389183/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:2082504/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:2234645/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:2479736/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:2561837/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:2629179/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:3100386/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:425160/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:667626/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:900909/1‑69 (MQ=255)
aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                                 >  1:101581/1‑69 (MQ=255)
           tAACGCGACGAAAACGCCCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTgag                                                     <  1:2536807/69‑1 (MQ=255)
            aaCGCGACGAAATCGCCGCCAGCGTTCAACACCGTTGACGTTAGCGCCATAACGCTCACGTAATTTgag                                                     <  1:1167207/69‑1 (MQ=255)
            aaCGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTgag                                                     <  1:218324/69‑1 (MQ=255)
                                aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                 <  1:2450911/69‑1 (MQ=255)
                                aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                 <  1:1562906/69‑1 (MQ=255)
                                aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                 <  1:2846184/69‑1 (MQ=255)
                                aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                 <  1:1353016/69‑1 (MQ=255)
                                aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                                 <  1:3325288/69‑1 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACTTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:2404824/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTgaga                                                    >  1:1975993/1‑49 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTg                                      >  1:94951/1‑63 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGCCCAAt                                >  1:1492457/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCaa                                 >  1:1265655/1‑68 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:1664365/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:687173/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:2546606/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:895347/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:2589399/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:891384/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:3018104/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:3059620/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:818378/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:814195/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:3224354/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:757971/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:3355321/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:3429928/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:411442/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:624061/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:246128/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:2403023/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:1902827/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:1996005/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:201513/1‑69 (MQ=255)
                                 gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:22446/1‑69 (MQ=255)
                                 gCGTTCAACACAGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                                >  1:1746644/1‑69 (MQ=255)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:1564483/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:1899741/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:757350/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:578591/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:51494/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:1084989/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:3325019/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:327901/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:2272811/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:2984966/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:2701024/1‑38 (MQ=38)
                                     tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                           >  1:2387348/1‑38 (MQ=38)
                                          acCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGa                       <  1:134659/69‑1 (MQ=255)
                                              aTGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATcc                    <  1:626901/68‑1 (MQ=255)
                                                gACGTTTGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCgg                  >  1:441444/1‑68 (MQ=255)
                                                gACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCggg                 >  1:597353/1‑69 (MQ=255)
                                                gACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCggg                 >  1:2054223/1‑69 (MQ=255)
                                                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa            <  1:2425349/69‑1 (MQ=255)
                                                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa            <  1:526215/69‑1 (MQ=255)
                                                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa            <  1:1138439/69‑1 (MQ=255)
                                                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa            <  1:3380156/69‑1 (MQ=255)
                                                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa            <  1:2928920/69‑1 (MQ=255)
                                                     tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa            <  1:2084769/69‑1 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCg                        >  1:86768/1‑52 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:314851/1‑66 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:1022097/1‑66 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:1000889/1‑66 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:1897469/1‑66 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:1228404/1‑66 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:1482271/1‑66 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:3269526/1‑66 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:2946549/1‑66 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:1377213/1‑66 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:277893/1‑66 (MQ=255)
                                                          ccATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGAATa          >  1:3372213/1‑66 (MQ=255)
                                                           cATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCa      <  1:10652/69‑1 (MQ=255)
                                                           cATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCa      <  1:458030/69‑1 (MQ=255)
                                                           cATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCa      <  1:428290/69‑1 (MQ=255)
                                                            aTAACGCTCACGTAATTTGAGATCATCCAGCCTTTGACCAATCATTGGCGATCCGGGGCGAATAGCCa      <  1:741905/68‑1 (MQ=255)
                                                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGAc   >  1:1483738/1‑68 (MQ=255)
                                                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg  >  1:1218260/1‑69 (MQ=255)
                                                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg  >  1:2418544/1‑69 (MQ=255)
                                                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg  >  1:801023/1‑69 (MQ=255)
                                                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg  >  1:1701211/1‑69 (MQ=255)
                                                               aCGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACg  >  1:2832073/1‑69 (MQ=255)
                                                                  |                                                                 
AACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAATAGCCAGACG  >  minE/1506514‑1506645

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: