Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 330,966 T→G 68.3% intergenic (+246/‑629) betT → / → yahA choline transporter of high affinity/predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb330,9660TG68.3% 91.5 / 24.5 44intergenic (+246/‑629)betT/yahAcholine transporter of high affinity/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (6/8);  new base G (13/17);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

AGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTG  >  W3110S.gb/330903‑331032
                                                               |                                                                  
aGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAggt                                                                  >  1:2301714/1‑66 (MQ=255)
  tttAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCg                                                            >  1:756400/1‑70 (MQ=255)
   ttAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGt                                                           >  1:1764117/1‑70 (MQ=255)
    tAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTa                                                          <  1:1944454/70‑1 (MQ=255)
     aGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTaa                                                         <  1:2600578/70‑1 (MQ=255)
         ttATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATc                                                             >  1:1634585/1‑62 (MQ=255)
           aTTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAAc                                                   >  1:2708560/1‑70 (MQ=255)
           aTTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTAAATc                                                      <  1:1800089/67‑1 (MQ=255)
            tttAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTAAATc                                                      <  1:1500617/66‑1 (MQ=255)
              tAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTaaa                                                        >  1:2621362/1‑62 (MQ=255)
              tAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTAAATCAACTGa                                                <  1:2399607/70‑1 (MQ=255)
                 tGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATc                                                      <  1:1024384/61‑1 (MQ=255)
                 tGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTAAATCAACTGAAAtt                                            <  1:717959/71‑1 (MQ=255)
                   tAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATc                                                      >  1:840056/1‑59 (MQ=255)
                   tAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTAAATCAACTGAAATTa                                           <  1:2504859/70‑1 (MQ=255)
                       tatTGCTGATCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTAAATCAACTGAAAtt                                            <  1:1179254/65‑1 (MQ=255)
                        atTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCaa                                     >  1:2756241/1‑71 (MQ=255)
                               aTCATTTGAAATGACGCATTATTCATGAGAAAggtg                                                                 <  1:2746024/36‑1 (MQ=255)
                               aTCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTAAATCAACTGAAAtt                                            <  1:800914/57‑1 (MQ=255)
                                tCATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGt                             <  1:254696/71‑1 (MQ=255)
                                 cATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGt                             >  1:2314381/1‑70 (MQ=255)
                                 cATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTAAATCaa                                                    >  1:633187/1‑47 (MQ=255)
                                 cATTTGAAATGACGCATTATTCATGAGAAAGGTGTATCGTAAATCAACTGaaa                                              >  1:2345187/1‑53 (MQ=255)
                                   tttGAAATGACGCATTATTCATGAGAAATGTGTATCGTa                                                          <  1:327490/39‑1 (MQ=255)
                                       aaaTGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTaa                      <  1:424798/71‑1 (MQ=255)
                                       aaaTGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTaa                      <  1:408102/71‑1 (MQ=255)
                                        aaTGACGCATTATTCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAg                     <  1:1691304/71‑1 (MQ=255)
                                            aCGCATTATTCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGttt                 <  1:605403/71‑1 (MQ=255)
                                                 ttattCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGtt                            <  1:1906227/55‑1 (MQ=255)
                                                 ttattCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGtt                            <  1:2661663/55‑1 (MQ=255)
                                                 ttattCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTAt                          <  1:279496/57‑1 (MQ=255)
                                                 ttattCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGttt                 <  1:956653/66‑1 (MQ=255)
                                                   attCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTAt                          >  1:1339251/1‑55 (MQ=255)
                                                   attCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGt                             <  1:2410610/52‑1 (MQ=255)
                                                   attCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGTTTTAATTATc          >  1:1906556/1‑71 (MQ=255)
                                                   attCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTaa               >  1:206885/1‑66 (MQ=255)
                                                   attCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTaa               >  1:2074421/1‑66 (MQ=255)
                                                   attCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATc          >  1:1851106/1‑71 (MQ=255)
                                                   attCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATc          >  1:475461/1‑71 (MQ=255)
                                                   attCATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATc          >  1:839394/1‑71 (MQ=255)
                                                      cATGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCtg        <  1:1846907/70‑1 (MQ=255)
                                                       aTGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCtgtg      <  1:2360634/71‑1 (MQ=255)
                                                       aTGAGAAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCtgtg      <  1:808680/71‑1 (MQ=255)
                                                          agaAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCtgtgtgt   >  1:991285/1‑71 (MQ=255)
                                                          agaAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCtgtgtgt   >  1:456765/1‑71 (MQ=255)
                                                           gaAAGGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCtgtgtgtg  <  1:2620755/71‑1 (MQ=255)
                                                               tgtgTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGttt                 <  1:2657572/52‑1 (MQ=255)
                                                               |                                                                  
AGTTTAGTCTTATTTAATGTAAATATTGCTGATCATTTGAAATGACGCATTATTCATGAGAAATGTGTATCGTAAATCAACTGAAATTAACGCAACCATTTGTTATTTAAGGTTTAATTATCTGTGTGTG  >  W3110S.gb/330903‑331032

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: