Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 452,584 T→G 26.8% P183P (CCT→CCG phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE452,5840TG26.8% 64.2 / 31.1 56P183P (CCT→CCGphrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (16/25);  new base G (8/7);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAA  >  minE/452529‑452651
                                                       |                                                                   
gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                                           >  1:521613/1‑66 (MQ=255)
gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                                           >  1:2542580/1‑66 (MQ=255)
gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                                           >  1:2519245/1‑66 (MQ=255)
gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                                           >  1:2660126/1‑66 (MQ=255)
gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                                           >  1:612876/1‑66 (MQ=255)
gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                                           >  1:2031065/1‑66 (MQ=255)
gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                                           >  1:2737384/1‑66 (MQ=255)
    ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                                                    >  1:316003/1‑69 (MQ=255)
    ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                                                    >  1:2482986/1‑69 (MQ=255)
    ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                                                    >  1:798798/1‑69 (MQ=255)
    ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                                                    >  1:2247162/1‑69 (MQ=255)
    ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                                                    >  1:2243116/1‑69 (MQ=255)
    ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                                                    >  1:2131850/1‑69 (MQ=255)
    ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                                                    >  1:1023626/1‑69 (MQ=255)
    ccGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTg                                                    >  1:171342/1‑69 (MQ=255)
           gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGATCCCGCGCCATCCATTACGCTGAATTATc                                             <  1:347568/69‑1 (MQ=255)
           gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc                                             <  1:861349/69‑1 (MQ=255)
           gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc                                             <  1:3322340/69‑1 (MQ=255)
           gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc                                             <  1:1967639/69‑1 (MQ=255)
           gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc                                             <  1:3016770/69‑1 (MQ=255)
           gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc                                             <  1:2411340/69‑1 (MQ=255)
           gCGTCGCTGCACCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCGCGCCATCCATTACGCTGAATTATc                                             <  1:1356270/69‑1 (MQ=255)
            cGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATcc                                            >  1:2829707/1‑69 (MQ=255)
                   gcgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg                                     <  1:3282142/69‑1 (MQ=255)
                   gcgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg                                     <  1:1255017/69‑1 (MQ=255)
                   gcgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg                                     <  1:1499769/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt                                    <  1:2214036/67‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt                                    <  1:3371760/67‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt                                    <  1:2937722/67‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt                                    <  1:206218/67‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt                                    <  1:1754501/67‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt                                    <  1:3172490/67‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:584911/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:596540/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:869025/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:3330189/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:3240963/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:912557/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:1749968/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:2110332/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:2207286/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:2327249/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:2655848/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:2431938/69‑1 (MQ=255)
                      ccAAAAGTTCGTAGTAGCGGATCGATAGAGCCCT‑‑CCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:3069582/67‑1 (MQ=255)
                        aaaaGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:743626/67‑1 (MQ=255)
                                           tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCAtttcc               >  1:2363561/1‑65 (MQ=255)
                                           tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTccc              >  1:2850948/1‑68 (MQ=255)
                                           tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg             >  1:3107005/1‑69 (MQ=255)
                                           tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg             >  1:2632027/1‑69 (MQ=255)
                                           tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg             >  1:2788698/1‑69 (MQ=255)
                                           tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg             >  1:332464/1‑69 (MQ=255)
                                           tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg             >  1:294648/1‑69 (MQ=255)
                                           tCGATAGAGCCCGCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCg             >  1:3109170/1‑69 (MQ=255)
                                                  agCCCTCGCCATACATTACGCTGAATTATCCTCGTCAGTCt                                  <  1:2529100/41‑1 (MQ=255)
                                                      cTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGaaaa  <  1:820525/69‑1 (MQ=255)
                                                       |                                                                   
GATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAA  >  minE/452529‑452651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: