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

ATTCTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAACTATTGCGAAA  >  W3110S.gb/3374343‑3374475
                                                                   |                                                                  
aTTCTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATgg                                                                <  1:1498540/71‑1 (MQ=255)
  tCTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATgg                                                                >  1:556425/1‑69 (MQ=255)
  tCTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATgg                                                                >  1:1237365/1‑69 (MQ=255)
   cTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAg                                                             >  1:1704970/1‑71 (MQ=255)
   cTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAg                                                             >  1:1821627/1‑71 (MQ=255)
   cTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAg                                                             >  1:1881078/1‑71 (MQ=255)
       aaGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGAt                                                                  <  1:1720058/62‑1 (MQ=255)
          gTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAGGGtttt                                                       <  1:1104222/70‑1 (MQ=255)
          gTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAGGGtttt                                                       <  1:1154284/70‑1 (MQ=255)
          gTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAGGGtttt                                                       <  1:525477/70‑1 (MQ=255)
                gtgATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAGGGTTTTctct                                                   >  1:2236887/1‑68 (MQ=255)
                 tgATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGtt                                                         <  1:646372/61‑1 (MQ=255)
                 tgATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGtt                                                         <  1:1854863/61‑1 (MQ=255)
                 tgATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTtg                                                <  1:133391/70‑1 (MQ=255)
                   aTTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGtttt                                                       >  1:2406736/1‑61 (MQ=255)
                    ttAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAg                                                             >  1:2409197/1‑54 (MQ=255)
                        catcatCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCtcat                                        <  1:12440/71‑1 (MQ=255)
                         atcatcAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCtcatc                                       <  1:1034964/71‑1 (MQ=255)
                         atcatcAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCtcatc                                       <  1:1232951/71‑1 (MQ=255)
                         atcatcAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCtcatc                                       <  1:100919/71‑1 (MQ=255)
                            atcaGTGACATCCTATCACAGGATTGAAAGTAGGGGaaaa                                                                   >  1:1175607/1‑40 (MQ=255)
                            atcaGTGACATCCTATCACAGGATTGAAAGTAGGGGaaaa                                                                   >  1:2089032/1‑40 (MQ=255)
                            atcaGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAg                                                             >  1:13723/1‑46 (MQ=255)
                            atcaGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTcatcat                                     >  1:1276450/1‑70 (MQ=255)
                               aGTGACATCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTAc                                  <  1:236946/70‑1 (MQ=255)
                               aGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTtc                                                      >  1:351038/1‑50 (MQ=255)
                                   aCATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTa                                   >  1:2043874/1‑65 (MQ=255)
                                     aTCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAAtt                           <  1:1461985/71‑1 (MQ=255)
                                      tCCTATCACAGGATTGAAAGTAGGGGGAAGATGGCAGGGTTTTCTCTTTGAGCCTCATCATTACCAt                               <  1:1549240/67‑1 (MQ=255)
                                      tCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAGGGTTTtctc                                                    <  1:1674874/45‑1 (MQ=255)
                                      tCCTATCACAGGATTGAAAGTAGGGGAAGATGGCAGGGTTTtctc                                                    <  1:2085925/45‑1 (MQ=255)
                                      tCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCtcatc                                       <  1:567380/58‑1 (MQ=255)
                                        cTATCACAGGATTGAAAGTAGGGGAAGATGGCAGGGTTTTCTCtt                                                  >  1:653180/1‑45 (MQ=255)
                                        cTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAAc                        >  1:745702/1‑71 (MQ=255)
                                        cTATCACAGGAGTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTaa                         >  1:1536560/1‑70 (MQ=255)
                                           tCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACgg                      >  1:1634926/1‑70 (MQ=255)
                                             acaGGATTGAAAGTAGGGGAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGaa                    >  1:1736161/1‑70 (MQ=255)
                                             acaGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGaa                    >  1:807123/1‑70 (MQ=255)
                                                 gATTGAAAGTAGGGGAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAAtt               <  1:672081/71‑1 (MQ=255)
                                                 gATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAAt                            <  1:2104527/58‑1 (MQ=255)
                                                  aTTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAAtt               >  1:726846/1‑70 (MQ=255)
                                                   ttGAAAGTAGGGGAAGATGGCAGGGTTTTCTCTTtg                                                >  1:856558/1‑36 (MQ=255)
                                                     gAAAGTAGGGGAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAAc            <  1:2581800/70‑1 (MQ=255)
                                                     gAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAAc            <  1:1474011/70‑1 (MQ=255)
                                                      aaaGTAGGGGAAGATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAACt           <  1:256140/70‑1 (MQ=255)
                                                         gTAGGGGAAGATGGCAGGGTTTTCTCTTTGTGCCtcatc                                       <  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)
                                                                   |                                                                  
ATTCTGTAAGGTCAGTGTGATTAACATCATCAGTGACATCCTATCACAGGATTGAAAGTAGGGGAAAATGGCAGGGTTTTCTCTTTGTGCCTCATCATTACCATAATTAACGGAATAATTAACTATTGCGAAA  >  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: