Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 740,339 T→G 7.3% P137P (CCT→CCG phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb740,3390TG 7.3% 130.1 / 5.2 55P137P (CCT→CCGphrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (22/29);  new base G (2/2);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

ATGTGGAAGTTGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAAC  >  W3110S.gb/740269‑740409
                                                                      |                                                                      
aTGTGGAAGTTGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCt                                                                        <  1:3214158/71‑1 (MQ=255)
     gAAGTTGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTg                                                                       >  1:562549/1‑67 (MQ=255)
     gAAGTTGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGgcg                                                                    <  1:3226995/70‑1 (MQ=255)
       aGTTGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCgg                                                                 <  1:545677/71‑1 (MQ=255)
        gTTGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGt                                                                <  1:3221457/71‑1 (MQ=255)
        gTTGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGGGATCCTGCCGCCTGGCGCGGt                                                                <  1:3716980/71‑1 (MQ=255)
          tGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGTCGCCTGGCGCGGtga                                                              >  1:3112495/1‑71 (MQ=255)
          tGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGtga                                                              >  1:2989018/1‑71 (MQ=255)
          tGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGtga                                                              >  1:3433089/1‑71 (MQ=255)
            aaaGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGtgatg                                                            <  1:489395/71‑1 (MQ=255)
              agagCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGAc                                                          <  1:34667/71‑1 (MQ=255)
              agagCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGAc                                                          <  1:325043/71‑1 (MQ=255)
              agagCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGAc                                                          <  1:3386631/71‑1 (MQ=255)
              agagCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGAc                                                          <  1:3425575/71‑1 (MQ=255)
                agCGCTGCGTAACGTGGTGTCTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGtga                                                              >  1:2366021/1‑65 (MQ=255)
                     tGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCGTGGGCGGTGATGACTggc                                                      >  1:3826547/1‑67 (MQ=255)
                     tGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCGGGGGCGGTGATGACTggc                                                      >  1:1420604/1‑67 (MQ=255)
                       cGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAAt                                                   <  1:1974475/69‑1 (MQ=255)
                       cGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAAt                                                   <  1:2144321/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)
                                      gAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGtg                                                               <  1:900707/42‑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:1141119/1‑70 (MQ=255)
                                            tttGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTa                             >  1:1145700/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:2385090/71‑1 (MQ=255)
                                                      gCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAg                  <  1:3119577/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:3673319/1‑44 (MQ=255)
                                                                      tGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTa                             >  1:80085/1‑44 (MQ=255)
                                                                      tGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAAc  >  1:768015/1‑71 (MQ=255)
                                                                      |                                                                      
ATGTGGAAGTTGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAAC  >  W3110S.gb/740269‑740409

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