Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,296,088 C→T 48.0% E661E (GAG→GAA katG ← catalase/hydroperoxidase HPI(I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,296,0880CT48.0% 3.8 / 66.2 77E661E (GAG→GAAkatGcatalase/hydroperoxidase HPI(I)
Reads supporting (aligned to +/- strand):  ref base C (16/24);  new base T (19/18);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-02

GGCTGGCCGTAAATTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGT  >  minE/2296036‑2296149
                                                    |                                                             
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:243021/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:2844590/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:185885/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:347096/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:1758063/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:2357155/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:1669907/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:3366343/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:2403989/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:1394125/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:748161/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:748327/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:2491526/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:25868/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:2682495/69‑1 (MQ=255)
ggctggctGTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTCCGATTCgtcgg                                               <  1:1453774/69‑1 (MQ=255)
agctggctgTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:1851481/68‑1 (MQ=255)
 gctggctgTGTATTTCACTTCGCCAGTTTCGCGGTCACGGCCTTCGAACAGTTCTTTCGATTCgtcgg                                               <  1:500323/68‑1 (MQ=255)
        gTAAATTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTc                                       <  1:2455323/69‑1 (MQ=255)
        gTAAATTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTc                                       <  1:841872/69‑1 (MQ=255)
          aaaTTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCa                                     <  1:3193976/69‑1 (MQ=255)
           aaTTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCAc                                    >  1:2905066/1‑69 (MQ=255)
           aaTTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCAc                                    >  1:2685723/1‑69 (MQ=255)
           aaTTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCAc                                    >  1:2882285/1‑69 (MQ=255)
           aaTTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCAc                                    >  1:1703803/1‑69 (MQ=255)
           aaTTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCAc                                    >  1:1577549/1‑69 (MQ=255)
           aaTTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCAc                                    >  1:649113/1‑69 (MQ=255)
           aaTTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCAc                                    >  1:699056/1‑69 (MQ=255)
           aaTTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCAc                                    >  1:1343683/1‑69 (MQ=255)
           aaTATCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTc                                       >  1:2836257/1‑66 (MQ=255)
                   ttCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACg                            <  1:1145158/69‑1 (MQ=255)
                   ttCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACg                            <  1:3034661/69‑1 (MQ=255)
                   ttCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACg                            <  1:903574/69‑1 (MQ=255)
                       ccGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCa                          <  1:501727/67‑1 (MQ=255)
                       ccGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCa                          <  1:280381/67‑1 (MQ=255)
                       ccGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCa                          <  1:2979934/67‑1 (MQ=255)
                            ttCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCa                    >  1:2683138/1‑68 (MQ=255)
                            ttCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCa                    >  1:2686994/1‑68 (MQ=255)
                             tcgcggtcacggCCTTCGACCAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTATCTCATATCCa                    >  1:1313494/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACTCATATGCa                    >  1:447977/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACTCATATCCa                    >  1:2163641/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACTCATATCCa                    >  1:3180481/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACTCATATCCa                    >  1:2528658/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACTCATATCCa                    >  1:3368796/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACTCATATCCa                    >  1:2080417/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACTCATATCCa                    >  1:471745/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACTCATATCCa                    >  1:1205162/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACTCATATCCa                    >  1:1001910/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCa                    >  1:3189737/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCa                    >  1:1024511/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCa                    >  1:288540/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCa                    >  1:132886/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCa                    >  1:2670251/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCa                    >  1:2577386/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCa                    >  1:2401837/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCa                    >  1:2229925/4‑67 (MQ=255)
                             tcgcggtcacggCCTTCGAACAGTTCTTTCGATTCGTCGGACGCTTTCCACTCGTAACGCATATCCa                    >  1:298009/4‑67 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:989272/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:963460/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:1274872/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:1451223/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:1838977/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:1947280/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:2033042/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:2120422/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:2306004/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:2639859/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:3163485/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:3078595/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:2982339/67‑1 (MQ=255)
                                               aaCAGCTCTTTCGATTCGTCGATCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGt  <  1:2761512/67‑1 (MQ=255)
                                                aCAGCTCTTTCGATTCGGCGGTCGCTTTCCACTCGTAACGCa                          <  1:3278593/42‑1 (MQ=255)
                                                   gCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCa                 >  1:416745/1‑48 (MQ=255)
                                                   gCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCa                 >  1:576221/1‑48 (MQ=255)
                                                   gCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCa                 >  1:1532995/1‑48 (MQ=255)
                                                   gCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCa                 >  1:1344696/1‑48 (MQ=255)
                                                   gCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCa                 >  1:119679/1‑48 (MQ=255)
                                                    |                                                             
GGCTGGCCGTAAATTTCACTTCGCCGGTTTCACGGTCACGGCCTTCGAACAGCTCTTTCGATTCGTCGGTCGCTTTCCACTCGTAACGCATATCCAGCAAGTTCACGAAGAAGT  >  minE/2296036‑2296149

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: