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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,704,876	0	C	T	36.0%	52.2 / 48.3	50	M164I (ATG→ATA)	trkD	potassium transporter
Reads supporting (aligned to +/- strand):  ref base C (19/13);  new base T (7/11);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

ACCCCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGGGGAACTATCCA  >  W3110S.gb/3704815‑3704944                                                              |                                                                     aCCCCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGt                                                              <  1:1309718/70‑1 (MQ=255) aCCCCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGt                                                              <  1:701848/70‑1 (MQ=255)  ccccAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGc                                                            <  1:769879/71‑1 (MQ=255)     cAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCAt                                                         >  1:441054/1‑71 (MQ=255)       gTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGt                                                              >  1:2299755/1‑64 (MQ=255)             ccAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGttttt                                                   <  1:2849489/69‑1 (MQ=255)             ccAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTg                                                  <  1:31985/70‑1 (MQ=255)               aGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAAt                                               >  1:1646489/1‑71 (MQ=255)                      aaaaCCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCataaa                                         >  1:4456162/1‑70 (MQ=255)                           cAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAAtaataa                                   <  1:497032/71‑1 (MQ=255)                            aGGTCAGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                                  <  1:4807634/71‑1 (MQ=255)                            aGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                                  <  1:2917990/71‑1 (MQ=255)                               tCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTg                               <  1:4263775/71‑1 (MQ=255)                                 aGCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGtt                                                      <  1:537326/46‑1 (MQ=255)                                  gCATGATCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTg                                                  <  1:940755/49‑1 (MQ=255)                                    aTGATCGGCGCAAACAGCTTACCGACCATAGCCGTGc                                                            <  1:135946/37‑1 (MQ=255)                                    aTGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGaa                           >  1:2431531/1‑70 (MQ=255)                                      gATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAataata                                    <  1:1127528/59‑1 (MQ=255)                                       aTCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAaataat                                     >  1:4014752/1‑57 (MQ=255)                                       aTCGGCGCAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGa                        <  1:557662/70‑1 (MQ=255)                                       aTCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGa                        >  1:4897657/1‑70 (MQ=255)                                       aTCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGa                        >  1:3996980/1‑70 (MQ=255)                                        tCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCataaa                                         >  1:2244433/1‑52 (MQ=255)                                        tCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCataaa                                         >  1:3827857/1‑52 (MQ=255)                                           gcgcAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAAtt                   <  1:3139143/71‑1 (MQ=255)                                            cgcAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTg                                                  <  1:2603875/39‑1 (MQ=255)                                            cgcAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAAtt                   >  1:5103776/1‑70 (MQ=255)                                            cgcAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAAt                    >  1:1233664/1‑69 (MQ=255)                                             gcAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAAtt                   >  1:279473/1‑69 (MQ=255)                                             gcAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTga                 <  1:1232295/71‑1 (MQ=255)                                              cAAACAGCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTgag                <  1:971771/71‑1 (MQ=255)                                              cAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAAtt                   >  1:1676351/1‑68 (MQ=255)                                              cAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAAtt                   >  1:3005648/1‑68 (MQ=255)                                              cAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAAtt                   >  1:4278368/1‑68 (MQ=255)                                                 aCAGCTGACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTgaga               >  1:655599/1‑69 (MQ=255)                                                  cAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                                  >  1:134964/1‑49 (MQ=255)                                                  cAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTga                              >  1:4765163/1‑53 (MQ=255)                                                    gCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGt                                <  1:1335636/49‑1 (MQ=255)                                                    gCTTACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGt                                <  1:1889006/49‑1 (MQ=255)                                                     cTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGGGGaa         >  1:4454569/1‑71 (MQ=255)                                                     cTTAACGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATaaa                                  >  1:2412989/1‑46 (MQ=255)                                                       tACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTgag                             >  1:4680588/1‑49 (MQ=255)                                                       tACCGACTATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGGGGAAc        >  1:3963645/1‑70 (MQ=255)                                                        aCCGACTATAGCCGTGCCATGTTTTTGAATCataa                                          >  1:1325319/1‑35 (MQ=255)                                                        aCCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGGGGAACt       <  1:3612124/70‑1 (MQ=255)                                                        aCCGACCATAGCCGTGCCATGTTCTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGGGGAACt       <  1:4767728/70‑1 (MQ=255)                                                           gACCATAGGCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATaa                     >  1:3335608/1‑53 (MQ=255)                                                           gACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGGGGAACTATcc   <  1:1972636/71‑1 (MQ=255)                                                            aCCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGGGGaa         >  1:4572669/1‑64 (MQ=255)                                                            aCCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGGGGAACTATCCa  <  1:815823/71‑1 (MQ=255)                                                              |                                                                     ACCCCAGTCCTGCCAGAATCAAAAACCAGGTCAGCATGATCGGCGCAAACAGCTTACCGACCATAGCCGTGCCATGTTTTTGAATCATAAATAATAAAGTGAGAACGATAATTGAGAGGGGAACTATCCA  >  W3110S.gb/3704815‑3704944

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