Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 501,987 C→T 43.3% P148P (CCC→CCT modA → molybdate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE501,9870CT43.3% 14.6 / 37.8 53P148P (CCC→CCTmodAmolybdate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (15/15);  new base T (13/10);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

CACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAAC  >  minE/501938‑502051
                                                 |                                                                
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1730154/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:2715456/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:2758511/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:3119218/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:215007/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1762549/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1721745/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1681129/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:445187/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:485498/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1109510/3‑67 (MQ=255)
cgctgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1186386/3‑67 (MQ=255)
cgctgctgAAAGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCgaa                                                 >  1:1395921/3‑67 (MQ=255)
cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:1595590/1‑67 (MQ=255)
cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:1101447/1‑67 (MQ=255)
cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:1699066/1‑67 (MQ=255)
cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:2756262/1‑67 (MQ=255)
cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:604149/1‑67 (MQ=255)
cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:1778036/1‑67 (MQ=255)
cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:761543/1‑67 (MQ=255)
cACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCaaa                                                 >  1:3196363/1‑67 (MQ=255)
 aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga                                              >  1:150431/1‑69 (MQ=255)
 aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga                                              >  1:2077611/1‑69 (MQ=255)
 aCTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAaaga                                              >  1:1885744/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:1105023/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:1205917/69‑1 (MQ=255)
   tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:1232763/69‑1 (MQ=255)
   tgctgAATGGCGGTCGCCTCGCGGTTGGTGATCCGGAACATGTCCCTGCTGGCATTTATGCGAaagaag                                            <  1:2347321/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:1167438/69‑1 (MQ=255)
            gCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCaaa                                   <  1:2247919/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:2361496/69‑1 (MQ=255)
                                  tCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGct             <  1:2744315/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)
                                                 |                                                                
CACTGCTGAATGGCGGTCGCCTGGCGGTTGGCGATCCGGAACATGTTCCCGCTGGCATTTATGCAAAAGAAGCACTGCAAAAACTGGGCGCATGGGATACGCTCTCTCCGAAAC  >  minE/501938‑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: