Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,426,674 T→C 83.8% D393D (GAT→GAC trkG → potassium transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,426,6740TC83.8% 91.3 / 10.3 37D393D (GAT→GACtrkGpotassium transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base C (11/20);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TAAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTATACTGG  >  W3110S.gb/1426610‑1426741
                                                                |                                                                   
tAAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                                <  1:3005085/70‑1 (MQ=255)
tAAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCgtgt                                                                <  1:2462310/70‑1 (MQ=255)
  aaTCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAAt                                                             <  1:2032960/71‑1 (MQ=255)
  aaTCAGCCTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGAAAGATAGTTACAGACCGTGTAAt                                                             <  1:1048721/71‑1 (MQ=255)
      aGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAg                                                          <  1:724580/70‑1 (MQ=255)
        cTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCgtgt                                                                <  1:1640301/62‑1 (MQ=255)
        cTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCgtgt                                                                <  1:696308/62‑1 (MQ=255)
        cTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCgtgt                                                                <  1:3403519/62‑1 (MQ=255)
            ctcATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGt                                                    <  1:1054012/70‑1 (MQ=255)
               aTCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTa                                                   <  1:3062956/68‑1 (MQ=255)
               aTCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTa                                                   <  1:767113/68‑1 (MQ=255)
               aTCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTa                                                   <  1:548423/68‑1 (MQ=255)
               aTCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTa                                                   <  1:3688355/68‑1 (MQ=255)
               aTCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTa                                                   <  1:2764229/68‑1 (MQ=255)
               aTCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTa                                                   <  1:2091710/68‑1 (MQ=255)
                  ccAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGt                                             >  1:3091345/1‑71 (MQ=255)
                      agCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttct                                         <  1:2112862/71‑1 (MQ=255)
                       gCGTTGTTGAGTGTAAATATAGGAGGGAAGATAGTTACAGACCgtg                                                                 >  1:761106/1‑46 (MQ=255)
                               gAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttcttttttctt                                 <  1:567262/70‑1 (MQ=255)
                               gAGTGTAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttcttttttctt                                 <  1:3264448/70‑1 (MQ=255)
                                    tAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATgg                                                <  1:2775825/50‑1 (MQ=255)
                                    tAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATgg                                                <  1:3478115/50‑1 (MQ=255)
                                    tAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGtttcttttttctt                                 <  1:3082732/65‑1 (MQ=255)
                                    tAAATGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATAct                           >  1:1619470/1‑71 (MQ=255)
                                     aaaTGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATActc                          >  1:1025968/1‑71 (MQ=255)
                                     aaaTGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATAct                           >  1:3703593/1‑70 (MQ=255)
                                      aaTGTAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATActc                          >  1:879724/1‑70 (MQ=255)
                                          tAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCtt                      >  1:3310776/1‑70 (MQ=255)
                                          tAGGAGGGAAGATAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCtt                      >  1:198769/1‑70 (MQ=255)
                                                     aTAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttttt          <  1:2394174/71‑1 (MQ=255)
                                                     aTAGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGttttt          >  1:2095686/1‑71 (MQ=255)
                                                       aGTTACAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTta        >  1:3243164/1‑71 (MQ=255)
                                                           aCAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTATACt    >  1:3658021/1‑71 (MQ=255)
                                                           aCAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTATACt    >  1:2518092/1‑71 (MQ=255)
                                                           aCAGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTATACt    >  1:1340011/1‑71 (MQ=255)
                                                             aGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTATACTgg  <  1:3038515/71‑1 (MQ=255)
                                                             aGACCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTATACTgg  <  1:899078/71‑1 (MQ=255)
                                                                |                                                                   
TAAATCAGCTTTCTCATCCCAGAGCGTTGTTGAGTGTAAATGTAGGAGGGAAGATAGTTACAGATCGTGTAATGAGGTCTGTATGGAGTTTCTTTTTTCTTTATACTCTCTTCACGGTGTTTTTTATACTGG  >  W3110S.gb/1426610‑1426741

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: