Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 740,342 C→G 7.0% G138G (GGC→GGG phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb740,3420CG 7.0% 162.6 / 2.2 57G138G (GGC→GGGphrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base G (2/2);  total (27/30)
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.36e-01

GAAGTTGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGG  >  W3110S.gb/740274‑740411
                                                                    |                                                                     
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: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)
                                 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:2986453/1‑71 (MQ=255)
                                    ggATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCtt                                 <  1:569343/71‑1 (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:2390418/1‑71 (MQ=255)
                                        ttGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACg                             >  1:1222540/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)
                                                                   gcgcGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAAc    <  1:3520640/69‑1 (MQ=255)
                                                                   gcgcGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAAc    <  1:1175795/69‑1 (MQ=255)
                                                                   gcgcGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACgg  >  1:133318/1‑71 (MQ=255)
                                                                    |                                                                     
GAAGTTGAAAGAGCGCTGCGTAACGTGGTGTGTGAAGGATTTGATGACAGCGTGATCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGG  >  W3110S.gb/740274‑740411

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: