Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,374,409 A→G 56.3% intergenic (‑65/‑315) nanR ← / → dcuD DNA‑binding transcriptional dual regulator/predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,374,4090AG56.3% 20.2 / 59.7 48intergenic (‑65/‑315)nanR/dcuDDNA‑binding transcriptional dual regulator/predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (10/11);  new base G (11/16);  total (21/27)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

ATTCTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAACTATTGCGAAA  >  W3110S.gb/3374343‑3374475
                                                                   |                                                                  
aTTCTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATgg                                                                <  1:1498540/71‑1 (MQ=255)
  tCTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATgg                                                                >  1:556425/1‑69 (MQ=255)
  tCTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATgg                                                                >  1:1237365/1‑69 (MQ=255)
   cTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAg                                                             >  1:1704970/1‑71 (MQ=255)
   cTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAg                                                             >  1:1821627/1‑71 (MQ=255)
   cTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAg                                                             >  1:1881078/1‑71 (MQ=255)
       aaGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGAt                                                                  <  1:1720058/62‑1 (MQ=255)
          gTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAGGGtttt                                                       <  1:1104222/70‑1 (MQ=255)
          gTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAGGGtttt                                                       <  1:1154284/70‑1 (MQ=255)
          gTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAGGGtttt                                                       <  1:525477/70‑1 (MQ=255)
                gtgATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAGGGTTTTctct                                                   >  1:2236887/1‑68 (MQ=255)
                 tgATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGtt                                                         <  1:646372/61‑1 (MQ=255)
                 tgATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGtt                                                         <  1:1854863/61‑1 (MQ=255)
                 tgATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTtg                                                <  1:133391/70‑1 (MQ=255)
                   aTTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGtttt                                                       >  1:2406736/1‑61 (MQ=255)
                    ttAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAg                                                             >  1:2409197/1‑54 (MQ=255)
                        catcatCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCtcat                                        <  1:12440/71‑1 (MQ=255)
                         atcatcAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCtcatc                                       <  1:1034964/71‑1 (MQ=255)
                         atcatcAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCtcatc                                       <  1:1232951/71‑1 (MQ=255)
                         atcatcAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCtcatc                                       <  1:100919/71‑1 (MQ=255)
                            atcaGTGACATCCTATCACAGGATTGAAAGTA‑GGGGaaaa                                                                   >  1:1175607/1‑40 (MQ=255)
                            atcaGTGACATCCTATCACAGGATTGAAAGTA‑GGGGaaaa                                                                   >  1:2089032/1‑40 (MQ=255)
                            atcaGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAg                                                             >  1:13723/1‑46 (MQ=255)
                            atcaGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTcatcat                                     >  1:1276450/1‑70 (MQ=255)
                               aGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTAc                                  <  1:236946/70‑1 (MQ=255)
                               aGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTtc                                                      >  1:351038/1‑50 (MQ=255)
                                   aCATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTa                                   >  1:2043874/1‑65 (MQ=255)
                                     aTCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAAtt                           <  1:1461985/71‑1 (MQ=255)
                                      tCCTATCACAGGATTGAAAGTAGGGGGAAGATGGCAGGGTTTTCTCTTTGAGCCTCATCATTACCAt                               <  1:1549240/67‑1 (MQ=255)
                                      tCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAGGGTTTtctc                                                    <  1:1674874/45‑1 (MQ=255)
                                      tCCTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAGGGTTTtctc                                                    <  1:2085925/45‑1 (MQ=255)
                                      tCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCtcatc                                       <  1:567380/58‑1 (MQ=255)
                                        cTATCACAGGATTGAAAGTA‑GGGGAAGATGGCAGGGTTTTCTCtt                                                  >  1:653180/1‑45 (MQ=255)
                                        cTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAAc                        >  1:745702/1‑71 (MQ=255)
                                        cTATCACAGGAGTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTaa                         >  1:1536560/1‑70 (MQ=255)
                                           tCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACgg                      >  1:1634926/1‑70 (MQ=255)
                                             acaGGATTGAAAGTA‑GGGGAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGaa                    >  1:1736161/1‑70 (MQ=255)
                                             acaGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGaa                    >  1:807123/1‑70 (MQ=255)
                                                 gATTGAAAGTA‑GGGGAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAAtt               <  1:672081/71‑1 (MQ=255)
                                                 gATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAAt                            <  1:2104527/58‑1 (MQ=255)
                                                  aTTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAAtt               >  1:726846/1‑70 (MQ=255)
                                                   ttGAAAGTA‑GGGGAAGATGGCAGGGTTTTCTCTTtg                                                >  1:856558/1‑36 (MQ=255)
                                                     gAAAGTA‑GGGGAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAAc            <  1:2581800/70‑1 (MQ=255)
                                                     gAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAAc            <  1:1474011/70‑1 (MQ=255)
                                                      aaaGTA‑GGGGAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAACt           <  1:256140/70‑1 (MQ=255)
                                                         gTA‑GGGGAAGATGGCAGGGTTTTCTCTTTGTGCCtcatc                                       <  1:1602351/39‑1 (MQ=255)
                                                               ggAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAACTATTGCGaaa  <  1:508328/71‑1 (MQ=255)
                                                               ggAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAACTATTGCGaaa  <  1:2581163/71‑1 (MQ=255)
                                                               ggAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAACTATTGCGaaa  <  1:2044254/71‑1 (MQ=255)
                                                                gAAGATTGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAACTATTGCGaaa  >  1:2257253/1‑70 (MQ=255)
                                                                   aaTGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAAtt               <  1:1683250/54‑1 (MQ=255)
                                                                   |                                                                  
ATTCTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTA‑GGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAACTATTGCGAAA  >  W3110S.gb/3374343‑3374475

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