Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,844,803 G→A 64.9% T94T (ACC→ACT cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,844,8030GA64.9% 37.1 / 35.1 64T94T (ACC→ACTcysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base G (14/8);  new base A (27/14);  total (42/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-02

CCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTT  >  minE/1844737‑1844864
                                                                  |                                                             
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:2922765/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:2118874/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:2169801/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:2302887/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:242076/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:2471966/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:2671011/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:27478/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:2887897/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:1555895/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:1547381/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:1464000/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:381834/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:1194362/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:690198/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:851247/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:1113933/1‑69 (MQ=255)
cccTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGt                                                             >  1:1006267/1‑69 (MQ=255)
            ggCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTc                                                 >  1:165683/1‑69 (MQ=255)
            ggCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTc                                                 >  1:531122/1‑69 (MQ=255)
                  aataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATTGt                                                             >  1:2238042/1‑51 (MQ=255)
                  aataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGt                                                             >  1:3156448/1‑51 (MQ=255)
                  aataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGt                                                             >  1:3239518/1‑51 (MQ=255)
                  aataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGt                                                             >  1:353037/1‑51 (MQ=255)
                  aataaATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGt                                                             >  1:410587/1‑51 (MQ=255)
                        tttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                                     <  1:2901993/69‑1 (MQ=255)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                                     <  1:989999/68‑1 (MQ=255)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                                     <  1:1070175/68‑1 (MQ=255)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                                     <  1:1524828/68‑1 (MQ=255)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                                     <  1:2109447/68‑1 (MQ=255)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:3123852/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:1164686/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:739054/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:698926/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:65939/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:1321044/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:3365525/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:3062076/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:1571744/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:2025240/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:1865437/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:2549087/48‑1 (MQ=39)
                         ttACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:2510057/48‑1 (MQ=39)
                         ttACGCTTCTAGGTAGAGAAATAGCGGTAGGCCACATCAATAGTGATg                                                         <  1:1560726/48‑1 (MQ=255)
                           aCGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                                     <  1:398088/66‑1 (MQ=255)
                                              tAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGcc               <  1:1221712/69‑1 (MQ=255)
                                               aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                  >  1:2194153/1‑65 (MQ=255)
                                               aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                  >  1:950700/1‑65 (MQ=255)
                                               aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                  >  1:2250083/1‑65 (MQ=255)
                                               aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                  >  1:1874136/1‑65 (MQ=255)
                                               aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                  >  1:1178773/1‑65 (MQ=255)
                                               aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                  >  1:1874058/1‑65 (MQ=255)
                                               aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                  >  1:2482173/1‑65 (MQ=255)
                                               aGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCaga                  >  1:1369663/1‑65 (MQ=255)
                                                 cGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTg                                     <  1:1537350/44‑1 (MQ=255)
                                                           catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt  >  1:5442/3‑69 (MQ=255)
                                                           catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt  >  1:1174150/3‑69 (MQ=255)
                                                           catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt  >  1:2955002/3‑69 (MQ=255)
                                                           catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt  >  1:2063146/3‑69 (MQ=255)
                                                           catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt  >  1:906691/3‑69 (MQ=255)
                                                           catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt  >  1:1835966/3‑69 (MQ=255)
                                                           catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCtttt  >  1:2761820/3‑69 (MQ=255)
                                                           catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCttt   >  1:1212396/3‑68 (MQ=255)
                                                           catCAATAGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGTGCCAGATCCAGCttt   >  1:3222408/3‑68 (MQ=255)
                                                                  |                                                             
CCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTT  >  minE/1844737‑1844864

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: