Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,600,711	A→G	13.1%	S405P (TCA→CCA)	ydeV ←	predicted sugar kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,600,711	0	A	G	13.1%	147.6 / 15.4	62	S405P (TCA→CCA)	ydeV	predicted sugar kinase
Reads supporting (aligned to +/- strand):  ref base A (30/23);  new base G (3/5);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

CCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTC  >  W3110S.gb/1600641‑1600774                                                                       |                                                                ccGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGa                                                                 >  1:1418237/1‑71 (MQ=255)     cTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACa                                                             <  1:873466/71‑1 (MQ=255)     cTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACa                                                             <  1:3446550/71‑1 (MQ=255)     cTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTCCTGCAAGTTACACGCTGATACa                                                             <  1:3665245/71‑1 (MQ=255)        cAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAAt                                                           <  1:282225/70‑1 (MQ=255)        cAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATc                                                          <  1:105983/71‑1 (MQ=255)            gACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAAt                                                           >  1:4354884/1‑66 (MQ=255)                 aCGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGGTACAATCGCCGCAttt                                                 >  1:374449/1‑71 (MQ=255)                 aCGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGCTACAATCGCCGCAttt                                                 >  1:3179832/1‑71 (MQ=255)                 aCGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAAt                                                           >  1:2915941/1‑61 (MQ=255)                 aCGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCAttt                                                 >  1:425672/1‑71 (MQ=255)                  cGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCa                                                    <  1:4131247/67‑1 (MQ=255)                    aTGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGGTACAATCGCCGCATTttct                                              <  1:684213/71‑1 (MQ=255)                    aTGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGAt                                                                <  1:3407535/53‑1 (MQ=255)                    aTGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTttct                                              <  1:4158408/71‑1 (MQ=255)                     tGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTttctt                                             <  1:1550674/71‑1 (MQ=255)                       aaGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGGTACAATCGCCGCATTttcttc                                            <  1:2556349/70‑1 (MQ=255)                       aaGGATGAATATTCGAGAAATCAGCAATTTGCTCCAAGTTACACGCTGGTACAATCGCCGCATTttcttc                                            >  1:3597690/1‑70 (MQ=255)                           aTGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAAt                                                           >  1:3676086/1‑51 (MQ=255)                            tGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAgcg                                       <  1:3021532/70‑1 (MQ=255)                              aaTATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCAt                                                   >  1:1682604/1‑56 (MQ=255)                              aaTATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCAt                                                   >  1:1168135/1‑56 (MQ=255)                              aaTATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCa                                     >  1:4574618/1‑70 (MQ=255)                               atatTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACg                                   >  1:1678015/1‑71 (MQ=255)                               atatTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACg                                   >  1:542397/1‑71 (MQ=255)                                tatTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGGTACAATCGCCGCAtt                                                  <  1:1753913/55‑1 (MQ=255)                                tatTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCa                                          <  1:3703610/63‑1 (MQ=255)                                tatTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCTCg                                   >  1:3263784/1‑70 (MQ=255)                                tatTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCAc                                    >  1:4302901/1‑69 (MQ=255)                                tatTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCAGTTTCTTCCAGCGCa                                     >  1:1020537/1‑68 (MQ=255)                                tatTCGAGAAATCAGCAATATGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACg                                   >  1:4297468/1‑70 (MQ=255)                                  ttCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATa                                                               >  1:3757750/1‑40 (MQ=255)                                  ttCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATcg                                                         <  1:3146567/46‑1 (MQ=255)                                  ttCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTttctt                                             >  1:531265/1‑58 (MQ=255)                                  ttCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGaa                                >  1:1417041/1‑71 (MQ=255)                                   tCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGGTACAATCGCCGCATTTTCTTCCAGCGCACGGaa                                >  1:3939690/1‑70 (MQ=255)                                    cGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACACTCGCCGCATTTTCTTCCAGAGCACGGAAc                               <  1:1143604/70‑1 (MQ=255)                                    cGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAAc                               <  1:1015829/70‑1 (MQ=255)                                    cGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAAc                               <  1:4425836/70‑1 (MQ=255)                                     gagaAATCAGCAATTTGCTGCAAGTTACACGCTGGTACa                                                             <  1:2397975/39‑1 (MQ=255)                                     gagaAATCAGCAATTTGCTGCAAGTTACACGCTGATACaa                                                            >  1:2194647/1‑40 (MQ=255)                                     gagaAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACaa                             >  1:1236094/1‑71 (MQ=255)                                     gagaAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACaa                             <  1:2134058/71‑1 (MQ=255)                                     gagaAATCAGCAAATTGCTGCAAGTTACAGGCTGATACaa                                                            >  1:1637018/1‑40 (MQ=38)                                           tCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGc                       >  1:1435699/1‑71 (MQ=255)                                             aGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAAc                               <  1:1101726/61‑1 (MQ=255)                                              gCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTttc                                               >  1:3159538/1‑44 (MQ=255)                                                  tttGCTGCAAGTTACACGCTGATACAATCGCCGCATTttct                                              <  1:3135677/41‑1 (MQ=255)                                                  tttGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTGGCTTTGTTa                <  1:2601434/71‑1 (MQ=255)                                                   ttGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCa                                          >  1:1883409/1‑44 (MQ=255)                                                    tgctgcAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACaa                             >  1:1783174/1‑56 (MQ=255)                                                    tgctgcAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAAt                            >  1:2420681/1‑57 (MQ=255)                                                    tgctgcAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACa              >  1:2118019/1‑71 (MQ=255)                                                     gctgcAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACAt             >  1:3065026/1‑71 (MQ=255)                                                       tgcAAGTTACACGCTGGTACAATCGCCGCATTTTCTTCCAGCGCa                                     <  1:4122892/45‑1 (MQ=255)                                                         cAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTa          <  1:4046848/70‑1 (MQ=255)                                                          aaGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATc        >  1:3587779/1‑71 (MQ=255)                                                           aGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATc        >  1:373321/1‑70 (MQ=255)                                                            gTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTa                <  1:2187018/61‑1 (MQ=255)                                                             ttACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTAc               <  1:2114505/61‑1 (MQ=255)                                                                cacGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGt   >  1:983513/1‑70 (MQ=255)                                                                cacGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTc  <  1:526167/71‑1 (MQ=255)                                                                 acGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCg      >  1:3422311/1‑66 (MQ=255)                                                                       |                                                                CCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTACATTTATCCGGGTC  >  W3110S.gb/1600641‑1600774

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