Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	58,662	G→A	28.1%	V63V (GTG→GTA)	yabP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	58,662	0	G	A	28.1%	94.0 / 39.7	57	V63V (GTG→GTA)	yabP	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/18);  new base A (6/10);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGAGACGTTACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAATGAGA  >  W3110S.gb/58593‑58726                                                                      |                                                                 ggAGACGTTACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTg                                                                  >  1:1947477/1‑70 (MQ=255)     aCGTTACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATaa                                                             <  1:351783/71‑1 (MQ=255)     aCGTTACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATaa                                                             >  1:1153073/1‑71 (MQ=255)        ttACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCa                                                          <  1:4465207/71‑1 (MQ=255)          aCGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAg                                                         <  1:4123920/70‑1 (MQ=255)              cATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTgg                                                    >  1:5091940/1‑71 (MQ=255)               aTCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTGGc                                                   <  1:968407/71‑1 (MQ=255)                tCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGc                                                   >  1:1520424/1‑70 (MQ=255)                  tGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCaa                                                 >  1:960555/1‑70 (MQ=255)                    ttCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACAt                                              >  1:337963/1‑71 (MQ=255)                      cTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATc                                             >  1:2273949/1‑70 (MQ=255)                      cTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATc                                             >  1:2004723/1‑70 (MQ=255)                      cTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATc                                             <  1:2502693/70‑1 (MQ=255)                      cTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATc                                             <  1:2101777/70‑1 (MQ=255)                       tGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATc                                             <  1:4413043/69‑1 (MQ=255)                            aGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCt                                      >  1:2300881/1‑71 (MQ=255)                            aGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCt                                      >  1:4026062/1‑71 (MQ=255)                            aGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCt                                      >  1:198892/1‑71 (MQ=255)                              aTAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAAt                                                               <  1:4881340/44‑1 (MQ=255)                              aTAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAAt                                                               <  1:5001897/44‑1 (MQ=255)                              aTAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCt                                      >  1:139164/1‑69 (MQ=255)                               tAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCt                                      >  1:3907159/1‑68 (MQ=255)                               tAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAg                                   <  1:1752932/71‑1 (MQ=255)                               tAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAg                                   <  1:3606948/71‑1 (MQ=255)                                aaaCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGa                                                        >  1:2470186/1‑49 (MQ=255)                                aaaCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCaa                                                 >  1:4796462/1‑56 (MQ=255)                                 aaCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATaa                                                             <  1:1889874/43‑1 (MQ=255)                                 aaCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTg                                                     <  1:504239/51‑1 (MQ=255)                                  aCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCaa                                           >  1:1839757/1‑60 (MQ=255)                                  aCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAAt                                <  1:1599564/71‑1 (MQ=255)                                  aCAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTAAGAAt                                <  1:1037681/71‑1 (MQ=255)                                        tGAGGCTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACttt                           >  1:937446/1‑70 (MQ=255)                                             cTTATAAACTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTggg                     >  1:55375/1‑71 (MQ=255)                                               tataAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACt                        <  1:440707/66‑1 (MQ=255)                                                ataAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCa                                            >  1:2570286/1‑45 (MQ=255)                                                ataAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGc                                        >  1:3585695/1‑49 (MQ=255)                                                ataAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGaa                                 >  1:4649493/1‑56 (MQ=255)                                                ataAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGaa                                 >  1:4131662/1‑56 (MQ=255)                                                 taAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAg                                   <  1:1876079/53‑1 (MQ=255)                                                  aaaCTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTAt                >  1:3525497/1‑71 (MQ=255)                                                   aaCTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATc               <  1:3485213/71‑1 (MQ=255)                                                   aaCTTCTGTTCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATc               <  1:1291644/71‑1 (MQ=255)                                                   aaCTACTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTAAGAATACTTTACTg                       >  1:3026120/1‑63 (MQ=255)                                                    aCTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGaa                                 >  1:3344603/1‑52 (MQ=255)                                                    aCTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAAt                                >  1:1462082/1‑53 (MQ=255)                                                           ttCAATTTCGTGAATAACCAGACTGTCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACgag         <  1:4141216/69‑1 (MQ=255)                                                           ttCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAg                                   <  1:4327178/43‑1 (MQ=255)                                                           ttCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAg                                   <  1:3603387/43‑1 (MQ=255)                                                           ttCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAg                                   <  1:2529506/43‑1 (MQ=255)                                                           ttCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACgag         <  1:4002949/69‑1 (MQ=255)                                                           ttCAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTggg                     <  1:981588/57‑1 (MQ=255)                                                            tCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAAt      >  1:3327266/1‑71 (MQ=255)                                                             cAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACt                        >  1:4729302/1‑52 (MQ=255)                                                             cAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACt                        >  1:3472698/1‑52 (MQ=255)                                                             cAATTTCGTAAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAAt      <  1:2665142/70‑1 (MQ=255)                                                                tttCGTGAATAACCAGACTGGCAACATCAACGCCTCAGaa                                 <  1:3547420/40‑1 (MQ=255)                                                                tttCGTAAATAACCAGACTGGCAACATCAACGCCTCAGaa                                 <  1:2743869/40‑1 (MQ=255)                                                                  tCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAATgaga  >  1:3100369/1‑69 (MQ=255)                                                                      |                                                                 GGAGACGTTACGCCATCTGTTCTGGTCAGATAAACAGACTGAGGCTTATAAACTTCTGTTCAATTTCGTGAATAACCAGACTGGCAACATCAACGCCTCAGAATACTTTACTGGGGCTATCAACGAGAATGAGA  >  W3110S.gb/58593‑58726

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