Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,477,861	C→T	15.0%	P335L (CCG→CTG)	ydbD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,477,861	0	C	T	15.0%	101.7 / 9.5	40	P335L (CCG→CTG)	ydbD	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (15/19);  new base T (2/4);  total (17/23)
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

GAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATC  >  W3110S.gb/1477791‑1477923                                                                       |                                                               gagaTTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAAcc                                                                >  1:739782/1‑71 (MQ=255)             ttGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGa                                                    <  1:239917/71‑1 (MQ=255)                aCCGATCAGCATTTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAaat                                                 <  1:2145527/71‑1 (MQ=255)                   gATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACTGCTCTCTGTCGAaataat                                              <  1:3556206/71‑1 (MQ=255)                   gATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAaataat                                              <  1:509307/71‑1 (MQ=255)                      cAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGAc                                           >  1:3422886/1‑71 (MQ=255)                      cAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGAc                                           >  1:3258344/1‑71 (MQ=255)                      cAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGAc                                           >  1:2995993/1‑71 (MQ=255)                       aGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAaataa                                               >  1:955200/1‑66 (MQ=255)                       aGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAaataa                                               >  1:85672/1‑66 (MQ=255)                       aGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAaataa                                               >  1:3978024/1‑66 (MQ=255)                         cAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACTGCTCTCTGTCGAAATAATGAcgcg                                        <  1:2345497/71‑1 (MQ=255)                              gggTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAAt                                     >  1:3103935/1‑69 (MQ=255)                              gggTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAAt                                     >  1:2037094/1‑69 (MQ=255)                               ggTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACTGCTCTCTGTCGAAATAATGACGCGAATGAg                                  <  1:3761002/71‑1 (MQ=255)                                 ttACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAgtc                                <  1:2090832/71‑1 (MQ=255)                                  tACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGaa                                                   <  1:3378684/51‑1 (MQ=255)                                    ccAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAgt                                 >  1:3154694/1‑67 (MQ=255)                                     cAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCt                            <  1:2602515/71‑1 (MQ=255)                                     cAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCt                            <  1:1721306/71‑1 (MQ=255)                                      aGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTa                           <  1:1203203/71‑1 (MQ=255)                                      aGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTa                           <  1:500299/71‑1 (MQ=255)                                       gCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGAcgcg                                        <  1:1175950/57‑1 (MQ=255)                                       gCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGAcgcg                                        <  1:1469236/57‑1 (MQ=255)                                       gCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTAc                          >  1:1873379/1‑71 (MQ=255)                                            ttgttCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGAcgcg                                        <  1:1946513/52‑1 (MQ=255)                                            ttgttCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGAcgcg                                        <  1:302430/52‑1 (MQ=255)                                            ttgttCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGAcgcg                                        <  1:893192/52‑1 (MQ=255)                                            ttgttCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGaa                                      >  1:313152/1‑54 (MQ=255)                                                    gggTTGCTCGACCGTGAACTGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCgg             <  1:252270/71‑1 (MQ=255)                                                    gggTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGgtgt                <  1:385873/68‑1 (MQ=255)                                                         gCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCt                        >  1:2956443/1‑55 (MQ=255)                                                          cTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGtt       <  1:3684648/71‑1 (MQ=255)                                                          cTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGtt       <  1:702948/71‑1 (MQ=255)                                                            cGACCGTGAACTGCTCTCTGTCGAAATAATGACGCGAATGAgtcg                               >  1:3628585/1‑45 (MQ=255)                                                            cGACCGTGAACTGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGgt                  >  1:1530667/1‑58 (MQ=255)                                                            cGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTg                       >  1:640877/1‑53 (MQ=255)                                                            cGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTc           >  1:1470532/1‑65 (MQ=255)                                                             gACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGa     <  1:1305644/70‑1 (MQ=255)                                                                   gAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATc  <  1:2612394/67‑1 (MQ=255)                                                                    aaCCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAAt   >  1:3808712/1‑65 (MQ=255)                                                                      ccGCTCTCTGTCGAAATAATGACGCTAATGAGTCGTCt                            >  1:3299014/1‑38 (MQ=255)                                                                      ccGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCg              >  1:2586914/1‑52 (MQ=255)                                                                       |                                                               GAGATTCTGGAGTTGACCGATCAGCAATTGGGTTACCAGCCTGTTGTTCTGGGGTTGCTCGACCGTGAACCGCTCTCTGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATC  >  W3110S.gb/1477791‑1477923

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