Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,981,022 C→T 14.6% E193E (GAG→GAA rob ← DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,981,0220CT14.6% 93.5 / 5.3 54E193E (GAG→GAArobDNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base T (5/3);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAT  >  minE/2980961‑2981081
                                                             |                                                           
ccGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTGTtcgtcg                                                       <  1:1049678/68‑1 (MQ=255)
ccGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTGTtcgtcg                                                       <  1:1817934/68‑1 (MQ=255)
ccGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTGTtcgtcg                                                       <  1:2835091/68‑1 (MQ=255)
     aGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTa                                                 <  1:2405859/69‑1 (MQ=255)
           tAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGa                                           >  1:2755874/1‑69 (MQ=255)
              ccATCTGCCTGATCGTGGGCTAACGCGGTGGTATAGAATACTTCTTGTTCGTCGTCTTTATCCTGACt                                         >  1:410369/1‑68 (MQ=255)
              ccATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTGTTCGTCGTCTTTATCCTGAc                                          >  1:3417363/1‑67 (MQ=255)
              ccATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTGTTCGTCGTCTTTATCCTGACt                                         >  1:100146/1‑68 (MQ=255)
              ccATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTGTTCGTCGTCTTTATCCTGACt                                         >  1:2479734/1‑68 (MQ=255)
              ccATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTGTTCGTCGTCTTTATCCTGACt                                         >  1:251903/1‑68 (MQ=255)
               cATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCg                                       <  1:2588809/69‑1 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:2857660/50‑1 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:1618062/50‑1 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:1661205/50‑1 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:1665090/50‑1 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:3027642/50‑1 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:1165685/50‑1 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:2201105/50‑1 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:840207/50‑1 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:225785/50‑1 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGAcg                                   >  1:293876/1‑69 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGAcg                                   >  1:782296/1‑69 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGAcg                                   >  1:865011/1‑69 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGAcg                                   >  1:2394348/1‑69 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGAc                                    >  1:3174214/1‑68 (MQ=255)
                   tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGAc                                    >  1:1564423/1‑68 (MQ=255)
                                   aaCGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTGATTCAgga                     >  1:1017225/1‑66 (MQ=255)
                                   aaCGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAgg                      >  1:1440780/1‑66 (MQ=255)
                                        ggtggtATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTagag              <  1:2216779/69‑1 (MQ=255)
                                        ggtggtATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTagag              <  1:2124200/69‑1 (MQ=255)
                                        ggtggtATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTagag              <  1:3084633/69‑1 (MQ=255)
                                         gtggtATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGc             <  1:2681892/69‑1 (MQ=255)
                                               taGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTagag              <  1:2315141/62‑1 (MQ=255)
                                               taGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcgg        >  1:2819108/1‑68 (MQ=255)
                                               taGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcgg        >  1:1648647/1‑68 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     >  1:2165934/1‑69 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:1083605/69‑1 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:469987/69‑1 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:1218108/69‑1 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     >  1:1420091/1‑69 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:3355215/69‑1 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:3267332/69‑1 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:2293936/69‑1 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:2861557/69‑1 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:1597935/69‑1 (MQ=255)
                                                 gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:194366/69‑1 (MQ=255)
                                                    tACCTCTTGCTCGTCGTCTTTATCCTGAGTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:3300949/1‑69 (MQ=255)
                                                    tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTa                 >  1:412771/1‑54 (MQ=255)
                                                    tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:3282462/1‑69 (MQ=255)
                                                    tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:2995992/1‑69 (MQ=255)
                                                    tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:2288283/1‑69 (MQ=255)
                                                    tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:518142/1‑69 (MQ=255)
                                                    tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:1066093/1‑69 (MQ=255)
                                                    tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:2374075/1‑69 (MQ=255)
                                                             |                                                           
CCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAT  >  minE/2980961‑2981081

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: