Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,910,142 C→T 58.9% V18I (GTT→ATT)  ygdB ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,910,1420CT58.9% 24.7 / 52.1 73V18I (GTT→ATT) ygdBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (3/27);  new base T (30/13);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 2.80e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGAATCCTG  >  minE/1910087‑1910206
                                                       |                                                                
gagGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAgg                                                     <  1:270717/69‑1 (MQ=255)
gagGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAgg                                                     <  1:1803969/69‑1 (MQ=255)
gagGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAgg                                                     <  1:3010259/69‑1 (MQ=255)
gagGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAgg                                                     <  1:1178540/69‑1 (MQ=255)
gagGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAgg                                                     <  1:2046845/69‑1 (MQ=255)
gagGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAgg                                                     <  1:1237226/69‑1 (MQ=255)
gagGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAgg                                                     <  1:2341575/69‑1 (MQ=255)
gagGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAgg                                                     <  1:423446/69‑1 (MQ=255)
 agGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAgg                                                     <  1:573321/68‑1 (MQ=255)
 agGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAgg                                                     <  1:2088097/68‑1 (MQ=255)
     taaaCTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCagcagc                                                           <  1:2439285/57‑1 (MQ=255)
     aaaaCTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCagcagc                                                           <  1:2366447/58‑1 (MQ=255)
     aaaaCTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCagcagc                                                           <  1:2451954/58‑1 (MQ=255)
     aaaaCTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCagcagc                                                           <  1:3423948/58‑1 (MQ=255)
     aaaaCTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCagcagc                                                           <  1:1082157/58‑1 (MQ=255)
         cTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGcc                                            <  1:567392/69‑1 (MQ=255)
         cTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGcc                                            <  1:2309184/69‑1 (MQ=255)
         cTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGcc                                            <  1:781117/69‑1 (MQ=255)
         cTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGcc                                            <  1:2453432/69‑1 (MQ=255)
         cTGCGATCCTGTTGGCTCATTCCTTGTAATAGAAAACTACCCAAAACCAGCAGCATCAGGACCAGAGcc                                            <  1:3080747/69‑1 (MQ=255)
                 cTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGa                                    <  1:1895284/69‑1 (MQ=255)
                 cTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGa                                    <  1:285289/69‑1 (MQ=255)
                 cTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGa                                    <  1:974340/69‑1 (MQ=255)
                 cTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGa                                    <  1:3323529/69‑1 (MQ=255)
                                     aTAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTca                <  1:2615731/69‑1 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCa                                           >  1:2508616/1‑41 (MQ=38)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:3247994/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:3286130/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:513600/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:617617/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:697386/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:77756/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:819088/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:860302/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:886244/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:906125/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:2834870/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:2712466/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:100498/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:2705720/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:982501/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:2606566/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:1044570/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:1307774/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:1311690/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:1433007/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:1484398/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:1689757/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:2066093/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:2067123/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:213011/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:2140217/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:2593714/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:2547373/1‑69 (MQ=255)
                                      tAGCAAGCTACCCAAAATCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcac               >  1:2306171/1‑69 (MQ=255)
                                      tAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTca                <  1:1994340/68‑1 (MQ=255)
                                       aGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcaca              >  1:2159154/1‑69 (MQ=255)
                                       aGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcaca              >  1:1180895/1‑69 (MQ=255)
                                       aGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTcaca              >  1:768120/1‑69 (MQ=255)
                                          aaaCTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAgg            <  1:2490382/68‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:2847182/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:1105111/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:1319136/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:2790496/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:1539851/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:3134533/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:1985715/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:2026489/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:355527/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:2283285/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:303310/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGGATCCTg  <  1:217901/69‑1 (MQ=255)
                                                   aaaaTCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACACGTTGGATCCTg  <  1:2701123/69‑1 (MQ=255)
                                                       |                                                                
GAGGCAAAACTGCGATCCTGTTGGCTCATTCCTTGTAATAGCAAACTACCCAAAACCAGCAGCATCAGGACCAGAGCCAGTGACGAAACACCCTTTTCGCGGTTCACAGGTTGAATCCTG  >  minE/1910087‑1910206

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: ‑