Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,041,067 C→A 9.0% intergenic (+118/+78) ygfZ → / ← yqfA predicted folate‑dependent regulatory protein/predicted oxidoreductase, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,041,0670CA 9.0% 204.8 / 4.8 66intergenic (+118/+78)ygfZ/yqfApredicted folate‑dependent regulatory protein/predicted oxidoreductase, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (24/36);  new base A (2/4);  total (26/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

ATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAA  >  W3110S.gb/3041002‑3041133
                                                                 |                                                                  
aTCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCAc                                                                  <  1:1399987/68‑1 (MQ=255)
aTCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCAc                                                                  <  1:997683/68‑1 (MQ=255)
  cATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCa                                                                   >  1:1507099/1‑65 (MQ=255)
   aTGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGaa                                                            <  1:2396033/71‑1 (MQ=255)
    tGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAAt                                                           <  1:6561896/71‑1 (MQ=255)
    tGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAAt                                                           <  1:7439787/71‑1 (MQ=255)
        tGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATcc                                                         <  1:6353216/69‑1 (MQ=255)
        tGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCgg                                                       >  1:2600312/1‑71 (MQ=255)
            aTTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATcc                                                         <  1:6018907/65‑1 (MQ=255)
            aTTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGt                                                   <  1:5419677/71‑1 (MQ=255)
             tttGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGtt                                                  <  1:661138/71‑1 (MQ=255)
             tttGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGtt                                                  >  1:6758347/1‑71 (MQ=255)
               tGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcc                                                               >  1:5625290/1‑56 (MQ=37)
                 aTGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGc                                              <  1:6216130/71‑1 (MQ=255)
                  tGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCt                                             >  1:7198019/1‑71 (MQ=255)
                  tGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCt                                             <  1:6972597/71‑1 (MQ=255)
                  tGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGAGTTCACGCCGAATCCGGCAGTtgtg                                               >  1:4866994/1‑69 (MQ=255)
                    ccTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCa                                                                   >  1:643948/1‑47 (MQ=255)
                     cTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCa                                                     >  1:2513226/1‑60 (MQ=255)
                     cTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCa                                                     >  1:7731168/1‑60 (MQ=255)
                     cTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCa                                                     >  1:653681/1‑60 (MQ=255)
                        tATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATg                                       >  1:1539361/1‑71 (MQ=255)
                        tATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTAACGCCGAATCCGGCAGTTGTGCTCCGATg                                       >  1:1494849/1‑71 (MQ=255)
                            acaacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTg                                   <  1:6239846/71‑1 (MQ=255)
                            acaacaGGGACGAGGCCCGTAGATCAGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCt                                    <  1:7365896/70‑1 (MQ=255)
                             caacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGcc                                     <  1:5229617/68‑1 (MQ=255)
                             caacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGcc                                     <  1:5909706/68‑1 (MQ=255)
                             caacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGcc                                     <  1:2517497/68‑1 (MQ=255)
                             caacaGGGACGAGGACCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGa                                         >  1:4316130/1‑64 (MQ=255)
                              aacaGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcct                                                              <  1:1932087/42‑2 (MQ=18)
                              aacaGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcct                                                              <  1:1873169/42‑2 (MQ=18)
                              aacaGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcct                                                              <  1:1859930/42‑2 (MQ=18)
                              aacaGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGcct                                                              <  1:5838786/42‑2 (MQ=18)
                              aacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGaa                                                            <  1:2839773/44‑1 (MQ=255)
                                caGGGACGAGGCCCGTAGATCGGATAAGGCGTTCa                                                                   >  1:6849127/1‑35 (MQ=255)
                                caGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAAt                                                           <  1:6754047/43‑1 (MQ=255)
                                caGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAAt                                                           <  1:3116404/43‑1 (MQ=255)
                                caGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAAt                                                           <  1:491414/43‑1 (MQ=255)
                                caGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGc                               >  1:5286024/1‑71 (MQ=255)
                                 aGGGACGAGGCCCGTAGATCGGATAAGGCGTTCAc                                                                  <  1:5069377/35‑1 (MQ=255)
                                  gggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGtt                                                  >  1:4486976/1‑50 (MQ=255)
                                    gACGAGTCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACg                           >  1:1014840/1‑71 (MQ=255)
                                      cGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCt                         >  1:6197303/1‑71 (MQ=255)
                                          gCCCGTCGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgc                      <  1:7732961/70‑1 (MQ=255)
                                          gCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTcc                                           <  1:7079246/49‑1 (MQ=255)
                                            ccGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCt                         <  1:4223146/65‑1 (MQ=255)
                                            ccGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCt                         <  1:2614173/65‑1 (MQ=255)
                                            ccGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGcgcg                   <  1:6229392/71‑1 (MQ=255)
                                            ccGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGcgcg                   >  1:6686788/1‑71 (MQ=255)
                                              gTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTc                 >  1:2554602/1‑71 (MQ=255)
                                              gTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGAGCGTc                 >  1:2493260/1‑71 (MQ=255)
                                               tAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgc                      >  1:2753443/1‑65 (MQ=255)
                                                 gATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCtt               <  1:2160858/70‑1 (MQ=255)
                                                 gATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTa              >  1:3436937/1‑71 (MQ=255)
                                                 gATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTa              >  1:3445948/1‑71 (MQ=255)
                                                     ggATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAg          <  1:5409450/71‑1 (MQ=255)
                                                     ggATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAg          <  1:4108870/71‑1 (MQ=255)
                                                       aTAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGc        <  1:5090092/71‑1 (MQ=255)
                                                           ggCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgcgc                    <  1:2445158/55‑1 (MQ=255)
                                                            gCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTgg                       >  1:3686811/1‑51 (MQ=255)
                                                            gCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAc    >  1:3731967/1‑70 (MQ=255)
                                                            gCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACa   <  1:353770/71‑1 (MQ=255)
                                                             cGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACaa  <  1:7340001/71‑1 (MQ=255)
                                                             cGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACaa  <  1:5180625/71‑1 (MQ=255)
                                                               ttCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTAt             <  1:161036/58‑1 (MQ=255)
                                                                tCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTAt             <  1:7729320/57‑1 (MQ=255)
                                                                 |                                                                  
ATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAA  >  W3110S.gb/3041002‑3041133

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