Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 740,531 T→C F201F (TTT→TTC phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb740,5310TC77.4% 114.1 / 29.4 62F201F (TTT→TTCphrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (12/2);  new base C (48/0);  total (60/2)
Fisher's exact test for biased strand distribution p-value = 4.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

ATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGC  >  W3110S.gb/740500‑740587
                               |                                                        
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCGGGTGGAAGAAAAAGCGGCGa                                    >  1:2733968/1‑54 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:2305946/1‑61 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:2107174/1‑61 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:2000188/1‑61 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:1906841/1‑61 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:2874385/1‑61 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:2192729/1‑61 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:57444/1‑61 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:1413537/1‑61 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:1397233/1‑61 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:1337470/1‑61 (MQ=255)
aTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCa                             >  1:1112301/1‑61 (MQ=255)
           tCTTTCTATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1977880/1‑55 (MQ=39)
           tCTTTCGATACTGCGCATTTCCCGTTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:53155/1‑55 (MQ=39)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAg                                          >  1:403132/1‑37 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGAt                                   >  1:605498/1‑44 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2655277/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2477407/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2722110/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2648094/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2643428/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2722648/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2552925/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2780278/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2373958/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2817205/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2863494/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2871385/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:298459/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:395014/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:537567/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:588118/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:621916/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:659573/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:797842/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:892384/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1944263/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1165596/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1188060/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1122766/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1455494/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1583914/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1663051/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1697059/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1848490/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1875898/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1884064/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2458769/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1977035/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2044028/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2147645/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2151146/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2171145/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2263061/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1040082/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2403788/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:2426722/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1077629/1‑55 (MQ=255)
           tCTTTCGATACTGCGCATTTCCCGGCGGAAGaaaaa                                           >  1:276573/1‑36 (MQ=37)
           tCTTTCGACACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTAc                        >  1:1906202/1‑55 (MQ=39)
                          cATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGc  <  1:1347110/62‑1 (MQ=255)
                          cATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGc  <  1:2445056/62‑1 (MQ=255)
                               |                                                        
ATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGC  >  W3110S.gb/740500‑740587

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: