Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,704,876 C→T 57.1% M164I (ATG→ATA trkD ← potassium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,704,8760CT57.1% 3.1 / 24.0 21M164I (ATG→ATAtrkDpotassium transporter
Reads supporting (aligned to +/- strand):  ref base C (5/4);  new base T (7/5);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAACCCCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGG  >  W3110S.gb/3704811‑3704933
                                                                 |                                                         
cGTAACCCCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACtata                                                        <  1:996380/69‑1 (MQ=255)
cGTAACCCCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACtata                                                        <  1:1538098/69‑1 (MQ=255)
 gTAACCCCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGcc                                                     <  1:1724470/71‑1 (MQ=255)
          gTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGt                                            <  1:1264188/71‑1 (MQ=255)
           tCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTg                                                  <  1:1132189/64‑1 (MQ=255)
                ccAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGa                                      <  1:953260/71‑1 (MQ=255)
                    aaTCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATc                                   <  1:315359/70‑1 (MQ=255)
                      tCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCat                                 >  1:862085/1‑70 (MQ=255)
                      tCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCat                                 >  1:511747/1‑70 (MQ=255)
                          aaaCCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCataaa                              >  1:405230/1‑69 (MQ=255)
                          aaaCCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCataaata                            <  1:1298394/71‑1 (MQ=255)
                                    aGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGt                     >  1:952619/1‑68 (MQ=255)
                                       aTGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCa                                  <  1:941009/52‑1 (MQ=255)
                                             ggCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCataa                               >  1:1336017/1‑49 (MQ=255)
                                             ggCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAAt         >  1:1675299/1‑71 (MQ=255)
                                                   aaCAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTga      >  1:1635169/1‑68 (MQ=255)
                                                   aaCAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTga      >  1:1054571/1‑68 (MQ=255)
                                                   aaCAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAAt         <  1:1198336/65‑1 (MQ=255)
                                                    aCAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                       >  1:988874/1‑50 (MQ=255)
                                                    aCAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAgg  >  1:63037/1‑71 (MQ=255)
                                                    aCAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAgg  >  1:639246/1‑71 (MQ=255)
                                                     cAGCTTACCGACTATAGCCGTGCCATGTTTTTGaa                                     >  1:417825/1‑35 (MQ=255)
                                                        cTTACCGACCATAGCCGTGCCATGTTTTTGAATCataaat                             <  1:197464/40‑1 (MQ=255)
                                                                 |                                                         
CGTAACCCCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGG  >  W3110S.gb/3704811‑3704933

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