Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 639,164 G→A 36.5% P63P (CCG→CCA ahpF → alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb639,1640GA36.5% 64.0 / 63.3 63P63P (CCG→CCAahpFalkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (13/27);  new base A (12/11);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGG  >  W3110S.gb/639095‑639226
                                                                     |                                                              
ggCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCg                                                                <  1:226317/70‑1 (MQ=255)
ggCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCg                                                                <  1:327782/70‑1 (MQ=255)
     aaaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCtt                                                            <  1:4100763/69‑1 (MQ=255)
     aaaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCtt                                                            <  1:345181/69‑1 (MQ=255)
      aaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTcc                                                         >  1:4592191/1‑71 (MQ=255)
      aaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTcc                                                         >  1:2008697/1‑71 (MQ=255)
        tCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTg                                                       <  1:1369392/71‑1 (MQ=255)
           cAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATc                                                    <  1:1449240/71‑1 (MQ=255)
           cAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATc                                                    <  1:1096710/71‑1 (MQ=255)
             gAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCAc                                                  >  1:4000960/1‑71 (MQ=255)
                cTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGAt                                                     <  1:4160202/65‑1 (MQ=255)
                cTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATc                                                    <  1:344292/66‑1 (MQ=255)
                 tGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTcc                                                         <  1:4820691/60‑1 (MQ=255)
                 tGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATc                                                    <  1:3159504/65‑1 (MQ=255)
                  gTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAAcc                                             >  1:1754604/1‑71 (MQ=255)
                   tCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAAccc                                            <  1:2948647/71‑1 (MQ=255)
                     aGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCa                                           <  1:2817748/70‑1 (MQ=255)
                     aGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCa                                           <  1:25896/70‑1 (MQ=255)
                     aGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAg                                          >  1:1098800/1‑71 (MQ=255)
                     aGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGAt                                                     <  1:424363/60‑1 (MQ=255)
                       aCAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGt                                        >  1:3625777/1‑71 (MQ=255)
                             tCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAAc                                  <  1:3652551/71‑1 (MQ=255)
                               aCCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAAcc                                 <  1:598613/70‑1 (MQ=255)
                               aCCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAAcc                                 <  1:4986879/70‑1 (MQ=255)
                               aCCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAAcc                                 <  1:3509433/70‑1 (MQ=255)
                               aCCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCa                                <  1:2498242/71‑1 (MQ=255)
                               aCCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCa                                <  1:1398004/71‑1 (MQ=255)
                                ccTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAg                               <  1:4062928/71‑1 (MQ=255)
                                ccTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAg                               <  1:4747360/71‑1 (MQ=255)
                                  tttAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCttt                                                           <  1:2432496/41‑1 (MQ=255)
                                  tttAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCttt                                                           <  1:3789173/41‑1 (MQ=255)
                                  tttAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCtt                                                            <  1:854525/40‑1 (MQ=255)
                                   ttAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCa                                >  1:3419715/1‑67 (MQ=255)
                                    tAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTg                                                       <  1:4751409/43‑1 (MQ=255)
                                     aaaGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGAt                                                     >  1:2200262/1‑44 (MQ=255)
                                     aaaGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGAt                                                     >  1:313809/1‑44 (MQ=255)
                                     aaaGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGAt                                                     >  1:4980615/1‑44 (MQ=255)
                                     aaaGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGAt                                                     >  1:2396407/1‑44 (MQ=255)
                                      aagaagATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGt                                        <  1:887209/56‑1 (MQ=255)
                                            aTAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATc                                                    <  1:2278854/38‑1 (MQ=255)
                                             tAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg                  <  1:4677731/71‑1 (MQ=255)
                                             tAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg                  <  1:159145/71‑1 (MQ=255)
                                             tAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg                  <  1:1451022/71‑1 (MQ=255)
                                              aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAgg                                         >  1:1993215/1‑47 (MQ=255)
                                               aCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGc                 <  1:712272/70‑1 (MQ=255)
                                                 aGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAgg              >  1:1772705/1‑71 (MQ=255)
                                                 aGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAgg              >  1:4459165/1‑71 (MQ=255)
                                                   cTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAgg                              >  1:2672397/1‑53 (MQ=255)
                                                     tGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAAc                                  >  1:4106694/1‑47 (MQ=255)
                                                      gCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTccc         >  1:4637350/1‑71 (MQ=255)
                                                          gTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCt     >  1:2501784/1‑71 (MQ=255)
                                                          gTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCt     >  1:4937157/1‑71 (MQ=255)
                                                          gTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCt     <  1:4858001/71‑1 (MQ=255)
                                                           tGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTg    <  1:5020408/71‑1 (MQ=255)
                                                           tGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg                  <  1:2164873/57‑1 (MQ=255)
                                                            gCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTg    >  1:253258/1‑70 (MQ=255)
                                                             cGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTggg  >  1:4538052/1‑71 (MQ=255)
                                                              gTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAAcc                                 <  1:3468170/39‑1 (MQ=255)
                                                              gTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAAcc                                 <  1:585486/39‑1 (MQ=255)
                                                               tAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAAcc                                 <  1:2100154/38‑1 (MQ=255)
                                                               tAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAAcc                                 >  1:2887507/1‑38 (MQ=255)
                                                               tAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAAcc                                 <  1:522888/38‑1 (MQ=255)
                                                               tAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAAcc                                 >  1:3838652/1‑38 (MQ=255)
                                                                  gCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAg                               >  1:1440633/1‑37 (MQ=255)
                                                                  gCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAg                               >  1:781042/1‑37 (MQ=255)
                                                                     gTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTcccc        <  1:2805820/57‑1 (MQ=255)
                                                                     |                                                              
GGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGG  >  W3110S.gb/639095‑639226

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: