Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 740,630 C→T 42.1% G234G (GGC→GGT phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb740,6300CT42.1% 19.6 / 44.2 38G234G (GGC→GGTphrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base T (10/6);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTGCATCGCTTG  >  W3110S.gb/740563‑740696
                                                                   |                                                                  
tACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTAc                                                                  >  1:3841901/1‑70 (MQ=255)
tACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTAc                                                                  >  1:428216/1‑70 (MQ=255)
tACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTAc                                                                  >  1:2598443/1‑70 (MQ=255)
tACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTAc                                                                  >  1:4758029/1‑70 (MQ=255)
     cAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCg                                                            <  1:4509037/71‑1 (MQ=255)
        ttttGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGttt                                                         <  1:4013196/71‑1 (MQ=255)
          ttGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGt                                                       >  1:3359475/1‑71 (MQ=255)
          ttGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGt                                                       >  1:3309815/1‑71 (MQ=255)
          ttGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGt                                                       >  1:2957972/1‑71 (MQ=255)
            gCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTc                                                      <  1:4736556/70‑1 (MQ=255)
                  aCGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGTCGGCCAGc                                               >  1:5072595/1‑71 (MQ=255)
                       gCCGGAGAATATGGGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGt                                                       >  1:4689273/1‑58 (MQ=255)
                        ccGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGTCGGCCAGCCTGGc                                          <  1:4261695/70‑1 (MQ=255)
                        ccGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGTCGGCCAGCCTGGc                                          <  1:4531649/70‑1 (MQ=255)
                        ccGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGTCGGCCAGCCTGGc                                          <  1:1453174/70‑1 (MQ=255)
                           gagaATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAAcg                                      <  1:2999518/71‑1 (MQ=255)
                             gaATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCGAGCCTGGCAAcggg                                    >  1:791686/1‑71 (MQ=255)
                             gaATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAAcggg                                    >  1:1696473/1‑71 (MQ=255)
                                 atGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGTCGGCCAGCCTGGCAAcgggcggg                                >  1:4261243/1‑71 (MQ=255)
                                  tGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGt                               >  1:3121030/1‑71 (MQ=255)
                                     gCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCg                                                            <  1:2657093/39‑1 (MQ=255)
                                     gCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTAt                            <  1:2836494/71‑1 (MQ=255)
                                       aacaacGAGATTTTCCGGCAGTGGAAGGTACCAGc                                                              <  1:1317729/35‑1 (MQ=255)
                                       aacaacGAGATTTTCCGGCAGTGGAAGGTACCAGc                                                              <  1:5143275/35‑1 (MQ=255)
                                       aacaacGAGATTTTCCGGCAGTGGAAGGTACCAGc                                                              <  1:4748408/35‑1 (MQ=255)
                                         caacGAGATTTTCCGGCAGTGGAAGGTACCAGCCg                                                            >  1:4975850/1‑35 (MQ=255)
                                           acGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTAtcgcctc                      <  1:2783353/71‑1 (MQ=255)
                                           acGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTAtcgcctc                      <  1:2155655/71‑1 (MQ=255)
                                            cGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAAc                                       >  1:2780836/1‑53 (MQ=255)
                                             gagaTTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTc                                                      >  1:4986845/1‑37 (MQ=255)
                                             gagaTTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTc                                                      >  1:4120061/1‑37 (MQ=255)
                                             gagaTTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTAtcgcctcgc                    <  1:3439398/71‑1 (MQ=255)
                                             gagaTTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAAAGGGCGGGTTAtcgcctcgc                    <  1:3420951/71‑1 (MQ=255)
                                                    tCCGGCAGTGGAAGGTACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCt             >  1:1358176/1‑71 (MQ=255)
                                                          aGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTGCATc       >  1:4486047/1‑71 (MQ=255)
                                                              gAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTGCATCGCtt   >  1:1502119/1‑71 (MQ=255)
                                                              gAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTGCATCGCt    >  1:1917487/1‑70 (MQ=255)
                                                              gAAGGCACCAGCCGTTTGTGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTGCATCGCttg  >  1:322466/1‑71 (MQ=255)
                                                                   |                                                                  
TACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTGCATCGCTTG  >  W3110S.gb/740563‑740696

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: