Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,628,292 T→C 63.8% V254V (GTT→GTC livJ → leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,628,2920TC63.8% 47.9 / 43.4 77V254V (GTT→GTClivJleucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (16/12);  new base C (38/11);  total (54/23)
Fisher's exact test for biased strand distribution p-value = 7.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCACGCGCGGCAGGGCTGAAAACTCAGTTTATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAACTACG  >  minE/2628233‑2628356
                                                           |                                                                
cAGGCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCt                                                         <  1:473343/69‑1 (MQ=25)
cAGGCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCt                                                         <  1:329331/69‑1 (MQ=25)
cAGGCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCt                                                         <  1:203066/69‑1 (MQ=25)
cAGGCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCt                                                         <  1:1836526/69‑1 (MQ=25)
cAGGCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCt                                                         <  1:3325578/69‑1 (MQ=25)
cAGGCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCt                                                         <  1:3034606/69‑1 (MQ=25)
cAGGCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCt                                                         <  1:151484/69‑1 (MQ=25)
cAGGCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCt                                                         <  1:1359807/69‑1 (MQ=25)
cAGGCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCt                                                         <  1:2840951/69‑1 (MQ=25)
cAGGCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCt                                                         <  1:927996/69‑1 (MQ=25)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTaag                                                      >  1:3142369/1‑68 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTAAc                                                      >  1:59129/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTAAc                                                      >  1:2864417/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTAAc                                                      >  1:1414155/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTAAc                                                      >  1:1346686/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTAAc                                                      >  1:1765653/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTAAc                                                      >  1:3396672/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTAAc                                                      >  1:1816423/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTAAc                                                      >  1:2092205/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTAAc                                                      >  1:1221242/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCTCGCTGTCTAAc                                                      >  1:2769484/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGTCGCGCTGTCTAAc                                                      >  1:3104365/1‑69 (MQ=21)
   gCACGTGCGGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGTGGCTAACGACTCGCTGTCTAAc                                                      >  1:1191590/1‑69 (MQ=21)
    cACGCGCGGCAGGGCTGAAAACTCAGTTTATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACa                                                     >  1:1904070/1‑69 (MQ=255)
    cACGCGCGGCAGGGCTGAAAACTCAGTTTATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACa                                                     >  1:325848/1‑69 (MQ=255)
    cACGCGCGGCAGGGCTGAAAACTCAGTTTATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACa                                                     >  1:2147406/1‑69 (MQ=255)
    cACGCGCGGCAGGGCTGAAAACTCAGTTTATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACa                                                     >  1:1692091/1‑69 (MQ=255)
          cgGCAGGGCTGAAAACCCAGTTTATGGGGCCGGAAGGCGGGGCTAACGTCTCGCTGTCt                                                         <  1:87247/59‑1 (MQ=21)
                      aaaCTCAGTTTATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGc                                     >  1:1830383/1‑67 (MQ=255)
                      aaaCTCAGTTTATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGc                                     >  1:492646/1‑67 (MQ=255)
                       aaCTCAGTTTATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGa                                  >  1:1591918/1‑69 (MQ=255)
                              tttATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGgc                             <  1:3281494/67‑1 (MQ=255)
                              tttATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGgc                             <  1:3166888/67‑1 (MQ=255)
                              tttATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGgc                             <  1:1792731/67‑1 (MQ=255)
                              tttATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGgc                             <  1:2803986/67‑1 (MQ=255)
                              tttATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTATAACATTGCGGGCGAATCAGCGGAAGGgc                             <  1:2197056/67‑1 (MQ=255)
                                           aaGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGc              <  1:2982148/69‑1 (MQ=255)
                                           aaGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGc              <  1:966135/69‑1 (MQ=255)
                                           aaGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGc              <  1:1299950/69‑1 (MQ=255)
                                           aaGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGc              <  1:2346607/69‑1 (MQ=255)
                                            aGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGc              <  1:2515148/68‑1 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGca                                    >  1:2327530/1‑45 (MQ=37)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:2730266/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:447446/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:1402698/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:570191/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:1783928/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:1382379/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:1930439/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:827712/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:2146829/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:2164269/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:893077/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:3067129/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:913585/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:3028793/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:922711/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:2411219/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:2441131/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:280416/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:2784736/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:2637848/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTCTAACATTGCGGGTGAATCAGCAGAAGG‑‑‑GCTGGTGACCAAGCCga           >  1:1550935/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTATAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:3234008/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTGTATAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:1775779/1‑68 (MQ=255)
                                            aGGCGTGGCTAACGTCTCGCTCTCTAACATTGCGGGTGAATCAGCAGAAGGGCTGCTGGTGACCAAGc              >  1:1177197/1‑68 (MQ=255)
                                                   gCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAAc      >  1:3336569/1‑69 (MQ=255)
                                                   gCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAAc      >  1:3296265/1‑69 (MQ=255)
                                                   gCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAAc      >  1:747179/1‑69 (MQ=255)
                                                   gCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAAc      >  1:2136572/1‑69 (MQ=255)
                                                   gCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAAc      >  1:2290426/1‑69 (MQ=255)
                                                   gCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAAc      >  1:305324/1‑69 (MQ=255)
                                                   gCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAAc      >  1:2411146/1‑69 (MQ=255)
                                                   gCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAAc      >  1:2794430/1‑69 (MQ=255)
                                                   gCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAAc      >  1:2521148/1‑69 (MQ=255)
                                                         gTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAACTACg  <  1:379559/67‑1 (MQ=255)
                                                         gTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAACTACg  <  1:1322212/67‑1 (MQ=255)
                                                           |                                                                
CAGGCACGCGCGGCAGGGCTGAAAACTCAGTTTATGGGGCCGGAAGGTGTGGCTAACGTTTCGCTGTCTAACATTGCGGGCGAATCAGCGGAAGGGCTGCTGGTGACCAAGCCGAAGAACTACG  >  minE/2628233‑2628356

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: