Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,602,757 T→A 31.3% L260F (TTA→TTT trkH ← potassium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,602,7570TA31.3% 82.4 / 57.3 64L260F (TTA→TTTtrkHpotassium transporter
Reads supporting (aligned to +/- strand):  ref base T (21/23);  new base A (8/12);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

GCGTAAACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTAATAG  >  W3110S.gb/3602687‑3602827
                                                                      |                                                                      
gCGTAAACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACtt                                                                        <  1:1094108/71‑1 (MQ=255)
 cGTAAACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACtt                                                                        <  1:311678/70‑1 (MQ=255)
   tAAACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACtt                                                                        <  1:3285115/68‑1 (MQ=255)
   tAAACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACtt                                                                        <  1:1161651/68‑1 (MQ=255)
     aaCTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACa                                                                    >  1:1073069/1‑70 (MQ=255)
     aaCTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAg                                                                   >  1:415674/1‑71 (MQ=255)
      tcTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACa                                                                    <  1:2829112/68‑1 (MQ=255)
          tACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTaa                                                                      >  1:320427/1‑63 (MQ=255)
          tACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTaa                                                                      >  1:3052586/1‑63 (MQ=255)
            cGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACa                                                                    >  1:2746753/1‑63 (MQ=255)
            cGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACa                                                                    >  1:2332657/1‑63 (MQ=255)
            cGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAg                                                            >  1:980109/1‑71 (MQ=255)
            cGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAg                                                            >  1:257588/1‑71 (MQ=255)
            cGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAg                                                            >  1:2031280/1‑71 (MQ=255)
              ccGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGtt                                                          <  1:3261516/71‑2 (MQ=255)
              ccGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:200269/71‑1 (MQ=255)
              ccGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:1348692/71‑1 (MQ=255)
               cGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:1771997/70‑1 (MQ=255)
               cGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:2537062/70‑1 (MQ=255)
                gATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:1126035/69‑1 (MQ=255)
                gATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:4022518/69‑1 (MQ=255)
                   aaaCATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGAc                                                     >  1:2730956/1‑71 (MQ=255)
                       aTGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTa                                                 <  1:811044/71‑1 (MQ=255)
                             aaTTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAg                                                >  1:285543/1‑66 (MQ=255)
                             aaTTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAg                                                >  1:3801609/1‑66 (MQ=255)
                             aaTTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACa                                           <  1:2813700/71‑1 (MQ=255)
                             aaTTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACa                                           <  1:520206/71‑1 (MQ=255)
                                tCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTaaa                                                                      >  1:3089910/1‑41 (MQ=255)
                                tCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTaa                                                                      >  1:1936238/1‑41 (MQ=255)
                                tCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTaa                                                                      >  1:1095420/1‑41 (MQ=255)
                                tCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTa                                                                       >  1:2024202/1‑40 (MQ=255)
                                 ccGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCg                                       >  1:701599/1‑71 (MQ=255)
                                    gATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGgag                                    <  1:1106959/71‑1 (MQ=255)
                                     aTCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGgag                                    <  1:2888/70‑1 (MQ=255)
                                      tCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGAt                                  <  1:795532/71‑1 (MQ=255)
                                           ccAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:3392117/42‑1 (MQ=255)
                                           ccAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:2343858/42‑1 (MQ=255)
                                           ccAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:1877418/42‑1 (MQ=255)
                                           ccAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:3063942/42‑1 (MQ=255)
                                             aaTAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCa                                                        >  1:2871351/1‑42 (MQ=255)
                                             aaTAAACCTTCAGACTACGCCCACTAAACAGTGGAAAGTGCAGACCGTAGTTACAGCCGGAGATcagcag                            >  1:3499311/1‑70 (MQ=255)
                                             aaTAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagcag                            >  1:3625325/1‑70 (MQ=255)
                                             aaTAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagcag                            >  1:419306/1‑70 (MQ=255)
                                             aaTAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcaccag                            >  1:3395442/1‑70 (MQ=255)
                                              aTAAACCTTCAGACTACGCCCACTTAACAGTGaaaa                                                             <  1:3790531/36‑1 (MQ=255)
                                               tAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:2578968/38‑1 (MQ=255)
                                               tAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGaa                         <  1:2478114/71‑1 (MQ=255)
                                               tAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGaa                         <  1:2836058/71‑1 (MQ=255)
                                               tAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGaa                         <  1:364655/71‑1 (MQ=255)
                                               tAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGaa                         <  1:3580994/71‑1 (MQ=255)
                                                 aaCCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAgg                           >  1:2306806/1‑67 (MQ=255)
                                                 aaCCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAgg                           >  1:2500473/1‑67 (MQ=255)
                                                 aaCCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAgg                           >  1:2534000/1‑67 (MQ=255)
                                                     ttCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGaa                         <  1:340653/65‑1 (MQ=255)
                                                        aGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAAt                <  1:3570258/71‑1 (MQ=255)
                                                        aGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAAt                <  1:3659422/71‑1 (MQ=255)
                                                           cTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGAt             <  1:164969/71‑1 (MQ=255)
                                                                 ccACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGtt       <  1:1053760/71‑1 (MQ=255)
                                                                 ccACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGtt       <  1:3222712/71‑1 (MQ=255)
                                                                 ccACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagcag                            <  1:2820260/50‑1 (MQ=255)
                                                                 ccACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagcag                            <  1:4003792/50‑1 (MQ=255)
                                                                 ccACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagcag                            <  1:2371622/50‑1 (MQ=255)
                                                                 ccACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagcag                            <  1:2197538/50‑1 (MQ=255)
                                                                 ccACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagcag                            <  1:1795814/50‑1 (MQ=255)
                                                                  cACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagca                             >  1:2541306/1‑48 (MQ=255)
                                                                    cTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGa              >  1:1256311/1‑61 (MQ=255)
                                                                    cTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGa              >  1:2505972/1‑61 (MQ=255)
                                                                    cTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTaa     >  1:827220/1‑70 (MQ=255)
                                                                    cTTAACAGTGAAAAGTCCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGa              >  1:2724716/1‑61 (MQ=255)
                                                                      tAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagca                             >  1:3368651/1‑44 (MQ=255)
                                                                      tAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagca                             >  1:3060631/1‑44 (MQ=255)
                                                                      tAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagca                             >  1:2124266/1‑44 (MQ=255)
                                                                      tAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTaa     >  1:1734175/1‑68 (MQ=255)
                                                                      tAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTAATAg  >  1:1225101/1‑71 (MQ=255)
                                                                      |                                                                      
GCGTAAACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTAATAG  >  W3110S.gb/3602687‑3602827

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: