Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 501,984 T→C 46.0% V147V (GTT→GTC modA → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE501,9840TC46.0% 2.0 / 40.3 50V147V (GTT→GTCmodAmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (12/15);  new base C (13/10);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAAC  >  minE/501916‑502051
                                                                    |                                                                   
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     >  1:2970461/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     <  1:2088448/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     >  1:1912252/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     >  1:2619027/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     <  1:3420608/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     >  1:1220561/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     <  1:121057/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     >  1:32705/1‑69 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     <  1:1079247/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     <  1:1075403/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     <  1:115067/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     <  1:3113121/69‑1 (MQ=255)
gACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGtt                                                                     >  1:1193379/1‑69 (MQ=255)
              cTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGAGCCGGAACATGtt                                                                     <  1:2994035/55‑1 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1681129/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:215007/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:485498/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:3119218/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:445187/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1762549/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1730154/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1721745/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:2715456/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1109510/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:2758511/3‑67 (MQ=255)
                      cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1186386/3‑67 (MQ=255)
                      cgctgctgAAAGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1395921/3‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:2756262/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:3196363/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:1101447/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:761543/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:1595590/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:1699066/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:1778036/1‑67 (MQ=255)
                      cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:604149/1‑67 (MQ=255)
                       aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga                                              >  1:2077611/1‑69 (MQ=255)
                       aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga                                              >  1:1885744/1‑69 (MQ=255)
                       aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga                                              >  1:150431/1‑69 (MQ=255)
                       aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga                                              >  1:1622763/1‑69 (MQ=255)
                         tgctgAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaagaag                                            <  1:281908/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaagaag                                            <  1:2418655/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaagaac                                            <  1:379301/69‑2 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTTCTGGCATTTATGCGAaagaag                                            <  1:2386204/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:2568018/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:2532615/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:2451247/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:2347321/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:1232763/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:1205917/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:1105023/69‑1 (MQ=255)
                         tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:1976817/69‑1 (MQ=255)
                         tgcagAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:261820/69‑1 (MQ=255)
                                  gctgTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCaaa                                   <  1:3384187/66‑1 (MQ=255)
                                  gCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCaaa                                   <  1:915316/69‑1 (MQ=255)
                                  gCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCaaa                                   <  1:2247919/69‑1 (MQ=255)
                                  gCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCaaa                                   <  1:1167438/69‑1 (MQ=255)
                                                  tGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGa                   <  1:3024011/69‑1 (MQ=255)
                                                        tCCGGAACATGTTCCCGCTGGCATTTATGCAAAATAAGCACTGCAAAAACTGGGCGCATGGGATACGct             <  1:1479826/69‑1 (MQ=255)
                                                        tCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGct             <  1:3119461/69‑1 (MQ=255)
                                                        tCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGct             <  1:3199469/69‑1 (MQ=255)
                                                        tCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGct             <  1:2744315/69‑1 (MQ=255)
                                                        tCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGct             <  1:2361496/69‑1 (MQ=255)
                                                        tCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGct             <  1:911842/69‑1 (MQ=255)
                                                         ccGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGct             <  1:3373481/68‑1 (MQ=255)
                                                                   ttCCCGCTGGCATTTATGCAAATGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAAc  >  1:2774709/1‑69 (MQ=255)
                                                                   ttCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAAc  >  1:2117815/1‑69 (MQ=255)
                                                                   ttCCCGCTGGCATTTATACAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAAc  >  1:2630307/1‑69 (MQ=255)
                                                                    |                                                                   
GACAGCAAAACCAACTGGACTTCACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAAC  >  minE/501916‑502051

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: