Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 452,736 C→T 17.8% A234V (GCA→GTA)  phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE452,7360CT17.8% 99.6 / 14.6 45A234V (GCA→GTA) phrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (16/21);  new base T (5/3);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTGCATCGCTTG  >  minE/452671‑452802
                                                                 |                                                                  
cGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAg                                                               >  1:713759/1‑71 (MQ=255)
   cAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGcc                                                             <  1:365885/70‑1 (MQ=255)
   cAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCg                                                            >  1:446346/1‑71 (MQ=255)
   cAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCg                                                            <  1:3550235/71‑1 (MQ=255)
   cAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCg                                                            <  1:3910847/71‑1 (MQ=255)
    aGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGt                                                           <  1:1453546/71‑1 (MQ=255)
    aGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGt                                                           >  1:3588599/1‑71 (MQ=255)
    aGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGt                                                           >  1:3584196/1‑71 (MQ=255)
          gCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTc                                                      <  1:1735815/70‑1 (MQ=255)
           ccAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGTCg                                                     <  1:2811370/70‑1 (MQ=255)
           ccAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCgg                                                    <  1:1401696/71‑1 (MQ=255)
            cAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGTCGGc                                                   <  1:199782/71‑1 (MQ=255)
                aCGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCAc                                                                  >  1:1204319/1‑52 (MQ=255)
                 cGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGcc                                              <  1:1498494/71‑1 (MQ=255)
                 cGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGcc                                              <  1:3329216/71‑1 (MQ=255)
                  ggTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCg                                                     >  1:1486226/1‑63 (MQ=255)
                  ggTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCt                                             <  1:920051/71‑1 (MQ=255)
                   gTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCt                                             >  1:3882247/1‑70 (MQ=255)
                     gCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGTCGGCCAGCCTgg                                           >  1:1407245/1‑70 (MQ=255)
                       cGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGc                                                                    <  1:3065150/43‑1 (MQ=255)
                       cGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGt                                                           <  1:3911745/52‑1 (MQ=255)
                       cGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTg                                                        <  1:2520568/55‑1 (MQ=255)
                         gagaATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAAcg                                      >  1:2998523/1‑71 (MQ=255)
                          agaATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAAcgg                                     >  1:2845571/1‑71 (MQ=255)
                          agaATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAAcgg                                     >  1:989566/1‑71 (MQ=255)
                              tatGAGCAACAACGAGATTTTCCGGCAGTGGAAGGTACCAGc                                                              <  1:2687481/42‑1 (MQ=255)
                                tGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGt                               >  1:1257126/1‑71 (MQ=255)
                                tGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGt                               >  1:858403/1‑71 (MQ=255)
                                 gAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGt                               >  1:89345/1‑70 (MQ=255)
                                   gCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTg                                                        <  1:2331408/43‑1 (MQ=255)
                                   gCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTg                                                        <  1:2258405/43‑1 (MQ=255)
                                   gCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTg                                                        <  1:521284/43‑1 (MQ=255)
                                      acaacGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTAtcgc                         <  1:1155460/71‑1 (MQ=255)
                                       caacGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTAtcgcc                        >  1:2531271/1‑71 (MQ=255)
                                       caacGAGATTTTCCGGCAGTGGAAGGTACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTAtcgcc                        >  1:795254/1‑71 (MQ=255)
                                          cGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTAtcgcctcg                     <  1:583063/71‑1 (MQ=255)
                                          cGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTAtcgc                         >  1:478744/1‑67 (MQ=255)
                                              aTTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAAcggg                                    <  1:1763976/52‑1 (MQ=255)
                                               ttttCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTAtcg                          >  1:3123005/1‑61 (MQ=255)
                                                  tCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGc              <  1:3329512/70‑1 (MQ=255)
                                                    cGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTg           <  1:2730293/71‑1 (MQ=255)
                                                    cGGCAGTGGAAGGCACCAGCCGTTTTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTg           <  1:1685301/70‑1 (MQ=255)
                                                         gTGGAAGGTACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTGCATCg      >  1:2536988/1‑71 (MQ=255)
                                                         gTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTg               <  1:930172/62‑1 (MQ=255)
                                                             aaGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTGCATCGCttg  >  1:244560/1‑71 (MQ=255)
                                                              aGGTACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCtt            >  1:1772084/1‑60 (MQ=255)
                                                                 |                                                                  
CGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGCACCAGCCGTTTGTCGGCCAGCCTGGCAACGGGCGGGTTATCGCCTCGCCAGTGCTTGCATCGCTTG  >  minE/452671‑452802

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: