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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,602,757	0	T	A	8.1%	161.3 / 6.9	62	L260F (TTA→TTT)	trkH	potassium transporter
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base A (3/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

TAAACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTAATAG  >  W3110S.gb/3602690‑3602827                                                                    |                                                                       tAAACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTaa                                                                      >  1:826861/1‑70 (MQ=255)     cTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGt                                                                  >  1:120219/1‑70 (MQ=255)      tGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTg                                                                 >  1:2109457/1‑70 (MQ=255)      tGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTg                                                                 >  1:277578/1‑70 (MQ=255)        tACTCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaaa                                                              <  1:2668330/71‑1 (MQ=255)        tACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaaa                                                              <  1:623309/71‑1 (MQ=255)        tACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGaa                                                               <  1:873052/70‑1 (MQ=255)        tACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGaa                                                               >  1:1004038/1‑70 (MQ=255)          cGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAg                                                            >  1:546351/1‑71 (MQ=255)            ccGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTg                                                          <  1:937118/71‑1 (MQ=255)              gATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCa                                                        <  1:431/71‑1 (MQ=255)               aTAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAg                                                            <  1:136651/66‑1 (MQ=255)                 aaaCATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGAc                                                     >  1:1083739/1‑71 (MQ=255)                 aaaCATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGAc                                                     >  1:265971/1‑71 (MQ=255)                    cATGCTAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGt                                                  <  1:863629/71‑1 (MQ=255)                    cATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCg                                                   >  1:1751388/1‑70 (MQ=255)                    cATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCg                                                   >  1:1032438/1‑70 (MQ=255)                    cATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTAAACAGTGAAAAGTGCAGACCGt                                                  <  1:1785228/71‑1 (MQ=255)                    cATGCGAAATTCCGGATCGCGCCAATAAACCTTAAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCg                                                   >  1:1256490/1‑70 (MQ=255)                      tGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACtt                                                                        <  1:2510600/47‑1 (MQ=255)                      tGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACtt                                                                        <  1:1722550/47‑1 (MQ=255)                        cGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTaa                                                                      >  1:440709/1‑47 (MQ=255)                        cGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGt                                               >  1:972642/1‑70 (MQ=255)                        cGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGt                                               >  1:1295970/1‑70 (MQ=255)                               ccGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAg                                                >  1:307900/1‑62 (MQ=255)                                  gATCGCGCCAATAAACCTTCAGACTACGCCCACtt                                                                        <  1:700393/35‑1 (MQ=255)                                  gATCGCGCCAATAAACCTTCAGACTACGCCCACtt                                                                        <  1:1640993/35‑1 (MQ=255)                                    tCGCGCCAATAAACCTTCAGACTACGCCCACTTAACa                                                                    <  1:1991077/37‑1 (MQ=255)                                     ctcgcCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGAt                                  >  1:1213285/3‑70 (MQ=255)                                     cgcgCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCgg                                      >  1:2227379/1‑66 (MQ=255)                                     cgcgCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGAt                                  >  1:1240201/1‑70 (MQ=255)                                     cgcgCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGAt                                  >  1:2640839/1‑70 (MQ=255)                                     cgcgCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGAt                                  >  1:78127/1‑70 (MQ=255)                                     cgcgCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGAt                                  >  1:658085/1‑70 (MQ=255)                                             tAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGaa                         <  1:1865011/71‑1 (MQ=255)                                             tAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGaa                         <  1:1507204/71‑1 (MQ=255)                                              aaaCCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagca                             >  1:1696693/1‑66 (MQ=255)                                               aaCCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAgg                           >  1:2208677/1‑67 (MQ=255)                                                aCCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGa                       <  1:703828/70‑1 (MQ=255)                                                 ccTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGa                       >  1:1803642/1‑69 (MQ=255)                                                 ccTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGa                       >  1:2762762/1‑69 (MQ=255)                                                 ccTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGa                       >  1:2341140/1‑69 (MQ=255)                                                  cTTCAGACTACGCCCACTTAACAGTGAAAAGTGCa                                                        >  1:784053/1‑35 (MQ=255)                                                  cTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcag                               <  1:1744382/60‑1 (MQ=255)                                                    tCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATcagca                             >  1:1930361/1‑60 (MQ=255)                                                     cAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGc                   <  1:2120803/69‑1 (MQ=255)                                                     cAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGc                   <  1:44758/69‑1 (MQ=255)                                                     cAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGc                   <  1:2400958/69‑1 (MQ=255)                                                     cAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGc                   <  1:2245111/69‑1 (MQ=255)                                                      aGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTa                                                 <  1:1286650/38‑1 (MQ=255)                                                      aGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAAt                <  1:718845/71‑1 (MQ=255)                                                       gACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATg               <  1:543046/71‑1 (MQ=255)                                                       gACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATg               <  1:2248999/71‑1 (MQ=255)                                                        aCTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTa                                             <  1:1035319/40‑1 (MQ=255)                                                        aCTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTa                                             <  1:2234102/40‑1 (MQ=255)                                                         cTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGga                                     >  1:2144733/1‑47 (MQ=255)                                                         cTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGAt             <  1:2141159/71‑1 (MQ=255)                                                          tACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGa              >  1:1085604/1‑69 (MQ=255)                                                           aCGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCa                  <  1:2221399/64‑1 (MQ=255)                                                               ccACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGtt       <  1:1497282/71‑1 (MQ=255)                                                               ccACTAAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGtt       <  1:2646751/71‑1 (MQ=255)                                                                cACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGtt       >  1:110517/1‑70 (MQ=255)                                                                 aCTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTaa     <  1:1926560/71‑1 (MQ=255)                                                                  cTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGgag                                    <  1:1364138/39‑1 (MQ=255)                                                                  cTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGgag                                    <  1:2185077/39‑1 (MQ=255)                                                                  cTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGga                                     <  1:1499032/38‑1 (MQ=255)                                                                   ttAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGgag                                    <  1:1201252/38‑1 (MQ=255)                                                                   ttAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGa                       <  1:1112687/51‑1 (MQ=255)                                                                    tAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGa                       >  1:27539/1‑50 (MQ=255)                                                                    tAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGa                       >  1:331889/1‑50 (MQ=255)                                                                    tAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTAATAg  >  1:1468546/1‑71 (MQ=255)                                                                    tAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTAATAg  >  1:1822832/1‑71 (MQ=255)                                                                    |                                                                       TAAACTGTACGCCGATAAACATGCGAAATTCCGGATCGCGCCAATAAACCTTCAGACTACGCCCACTTAACAGTGAAAAGTGCAGACCGTAGTTACAGCCGGAGATCAGCAGGAAGATAGCAATGATAGTGTTAATAG  >  W3110S.gb/3602690‑3602827

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: ‑