Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,981,016 C→T 16.2% E195E (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,0160CT16.2% 91.6 / 7.1 55E195E (GAG→GAArobDNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base T (6/3);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

ACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAT  >  minE/2980951‑2981081
                                                                 |                                                                 
aCCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACctct                                                                >  1:461423/1‑69 (MQ=255)
          ccGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTGTtcgtcg                                                       <  1:1049678/68‑1 (MQ=255)
          ccGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTGTtcgtcg                                                       <  1:2835091/68‑1 (MQ=255)
          ccGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTGTtcgtcg                                                       <  1:1817934/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)
                        ccATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACTTCTTgt                                                             >  1:2035147/1‑47 (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:1661205/50‑1 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:225785/50‑1 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:2201105/50‑1 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:1665090/50‑1 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:2857660/50‑1 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:1618062/50‑1 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:3027642/50‑1 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:1165685/50‑1 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTcgtcgt                                                      <  1:840207/50‑1 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGAcg                                   >  1:865011/1‑69 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGAcg                                   >  1:293876/1‑69 (MQ=255)
                             tGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGAcg                                   >  1:782296/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:3084633/69‑1 (MQ=255)
                                                  ggtggtATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTagag              <  1:2124200/69‑1 (MQ=255)
                                                  ggtggtATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTagag              <  1:2216779/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:2293936/69‑1 (MQ=255)
                                                           gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:3355215/69‑1 (MQ=255)
                                                           gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:3267332/69‑1 (MQ=255)
                                                           gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:1597935/69‑1 (MQ=255)
                                                           gAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACcggcgg     <  1:2861557/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:2288283/1‑69 (MQ=255)
                                                              tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:2374075/1‑69 (MQ=255)
                                                              tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:2995992/1‑69 (MQ=255)
                                                              tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:518142/1‑69 (MQ=255)
                                                              tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAt  >  1:1066093/1‑69 (MQ=255)
                                                                 |                                                                 
ACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAAT  >  minE/2980951‑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: