Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 740,354 C→T 7.4% T142T (ACC→ACT phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb740,3540CT 7.4% 157.3 / 2.3 54T142T (ACC→ACTphrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (23/27);  new base T (2/2);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

TGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGG  >  W3110S.gb/740290‑740422
                                                                |                                                                    
tGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCGTGGGCGGTGATGACTggc                                                                   >  1:3826547/1‑67 (MQ=255)
tGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCGGGGGCGGTGATGACTggc                                                                   >  1:1420604/1‑67 (MQ=255)
  cGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAAt                                                                <  1:2144321/69‑1 (MQ=255)
  cGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAAt                                                                <  1:1974475/69‑1 (MQ=255)
      aCGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTa                                                                  <  1:3002777/63‑1 (MQ=255)
      aCGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACga                                                           <  1:112297/70‑1 (MQ=255)
       cGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGTTAATCAc                                                             >  1:3154050/1‑67 (MQ=255)
       cGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGTTAATCAc                                                             >  1:3149627/1‑67 (MQ=255)
       cGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCAc                                                             >  1:3557174/1‑67 (MQ=255)
       cGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCAc                                                             >  1:3086921/1‑67 (MQ=255)
                tgAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACa                                                   <  1:380179/68‑1 (MQ=255)
                  aaGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTc                                              >  1:1759305/1‑71 (MQ=255)
                  aaGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATCACCGGTAATCACGAGATGTACAAAGTc                                              <  1:1603020/71‑1 (MQ=255)
                  aaGGATTTGATGACAGCGTGATCCTGCCGCCGGGGGCGGTGATGACTGGCAATCatg                                                            <  1:1813783/57‑3 (MQ=37)
                  aaGGATTTGATGACAGCGTGATCCTGCCGCCGCGGGCGGTGATGACTGGCAATCatg                                                            <  1:304951/57‑3 (MQ=255)
                   aGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCt                                             >  1:2044585/1‑71 (MQ=255)
                    ggttttGATGACAGCGTGCTCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACTAAGTCtt                                            >  1:1215023/4‑71 (MQ=255)
                    ggATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCtt                                            <  1:569343/71‑1 (MQ=255)
                    ggATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCtt                                            >  1:2986453/1‑71 (MQ=255)
                      aTTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCtt                                            <  1:1915733/69‑1 (MQ=255)
                       tttGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTa                                          >  1:1145700/1‑70 (MQ=255)
                       tttGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTa                                          >  1:1141119/1‑70 (MQ=255)
                        ttGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACg                                        >  1:1222540/1‑71 (MQ=255)
                        ttGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACg                                        >  1:2390418/1‑71 (MQ=255)
                         tgatgaCAGCGTGATCCTTCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGc                                       >  1:3537802/1‑71 (MQ=255)
                         tgatgaCAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGc                                       >  1:2205957/1‑71 (MQ=255)
                             gaCAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCttt                                   >  1:1913154/1‑71 (MQ=255)
                                aGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTaa                                <  1:2131144/71‑1 (MQ=255)
                                 gCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAg                               <  1:3119577/71‑1 (MQ=255)
                                 gCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAg                               <  1:2385090/71‑1 (MQ=255)
                                  cGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGa                              <  1:1969111/71‑1 (MQ=255)
                                   gTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGaa                             >  1:2024887/1‑71 (MQ=255)
                                      aTCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACg                                                            <  1:1570745/37‑1 (MQ=255)
                                      aTCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACa                                                   <  1:3036205/46‑1 (MQ=255)
                                      aTCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACa                                                   <  1:468215/46‑1 (MQ=255)
                                       tCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACg                                                            <  1:3232452/36‑1 (MQ=255)
                                        ccTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCt                        >  1:1373180/1‑71 (MQ=255)
                                         cTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGc                          <  1:843578/68‑1 (MQ=255)
                                          tGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGt                                               <  1:535084/46‑1 (MQ=255)
                                              gccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGa                  <  1:3339362/71‑1 (MQ=255)
                                                 tGGCGCGGTGATGACCTGTAATCACGAGATGTACAAAGTCTTTa                                          >  1:1227035/1‑44 (MQ=255)
                                                 tGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTa                                          >  1:80085/1‑44 (MQ=255)
                                                 tGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTa                                          >  1:3673319/1‑44 (MQ=255)
                                                 tGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAAc               >  1:768015/1‑71 (MQ=255)
                                                   gcgcGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAAc               <  1:3520640/69‑1 (MQ=255)
                                                   gcgcGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAAc               <  1:1175795/69‑1 (MQ=255)
                                                   gcgcGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACgg             >  1:133318/1‑71 (MQ=255)
                                                         tgatgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAg                               <  1:2178851/47‑1 (MQ=255)
                                                         tgatgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAg                               <  1:585089/47‑1 (MQ=255)
                                                         tgatgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCgg       >  1:3587410/1‑71 (MQ=255)
                                                           atgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCtgg      <  1:2279439/70‑1 (MQ=255)
                                                           atgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCggg      <  1:1047590/70‑1 (MQ=255)
                                                           atgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCggg      <  1:2723328/70‑1 (MQ=255)
                                                              aCCTGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAgg  <  1:1828719/71‑1 (MQ=255)
                                                               ccGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAgg  <  1:1881810/70‑1 (MQ=255)
                                                                cGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTaa                                >  1:2173142/1‑39 (MQ=255)
                                                                |                                                                    
TGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGG  >  W3110S.gb/740290‑740422

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: