Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,017,036	C→T	6.6%	A562A (GCG→GCA)	ycbZ ←	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,017,036	0	C	T	6.6%	181.1 / 2.1	60	A562A (GCG→GCA)	ycbZ	predicted peptidase
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base T (2/2);  total (29/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 = 8.59e-01

CGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATA  >  W3110S.gb/1016967‑1017103                                                                      |                                                                    cGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTgcg                                                                    >  1:3717739/1‑71 (MQ=255) cGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTgcg                                                                    >  1:1729761/1‑71 (MQ=255) cGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTgcg                                                                    >  1:512034/1‑71 (MQ=255)        aaCCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACg                                                             <  1:2319019/71‑1 (MQ=255)           cAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACg                                                             <  1:1710433/68‑1 (MQ=255)           cAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGt                                                            <  1:1627321/69‑1 (MQ=255)           cAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTc                                                          >  1:3809489/1‑71 (MQ=255)           cAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTc                                                          <  1:2166043/71‑1 (MQ=255)           cAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTc                                                          <  1:2722405/71‑1 (MQ=255)           cAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTc                                                          <  1:333925/71‑1 (MQ=255)           cAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTT‑CTGTTGCGATGCTTGCGCGATACGTTcc                                                         >  1:2445862/1‑71 (MQ=255)            acttCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTc                                                          <  1:2030440/68‑1 (MQ=255)            aGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTcc                                                         <  1:3384692/71‑1 (MQ=255)            aGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTcc                                                         >  1:1639073/1‑71 (MQ=255)            aGTTCAGCCAACGTAATGGCCATGGAAAACGGTGAAGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTcc                                                         >  1:3699520/1‑71 (MQ=255)                cAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGa                                                     >  1:784337/1‑71 (MQ=255)                cAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGa                                                     >  1:284794/1‑71 (MQ=255)                cAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGa                                                     >  1:2874179/1‑71 (MQ=255)                cAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGa                                                     >  1:1198708/1‑71 (MQ=255)                        gTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGc                                             <  1:2596045/71‑1 (MQ=255)                        gTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGc                                             <  1:3326884/71‑1 (MQ=255)                         tAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCa                                            <  1:459032/71‑1 (MQ=255)                            tGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCa                                         <  1:3041989/71‑1 (MQ=255)                            tGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCCTGCGCGATACGTTCCTGTATGGTTTGCATCa                                         <  1:2357280/71‑1 (MQ=255)                              gCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTg                                              <  1:2054929/64‑1 (MQ=255)                                cATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCa                                            >  1:2288700/1‑64 (MQ=255)                                  tGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGa                                                                 <  1:1666816/41‑1 (MQ=255)                                  tGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGa                                                                 <  1:1595182/41‑1 (MQ=255)                                  tGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGcgtcg                                   <  1:3173927/71‑1 (MQ=255)                                    gAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTc                                                          >  1:2056481/1‑46 (MQ=255)                                    gAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGtt                                 >  1:354422/1‑71 (MQ=255)                                    gAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGtt                                 >  1:1288221/1‑71 (MQ=255)                                     aaaaCGGTGACGTCCTTCCTGTTGCGATGCTTgcgc                                                                   <  1:2971541/36‑1 (MQ=255)                                        aCGGTGACGTCCTTCCTGTTGCGATGCTTGTGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGc                             >  1:1523994/1‑71 (MQ=255)                                        aCGGTGACGTCCTTCCTGTTGCGATGCTTGTGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGc                             >  1:2198450/1‑71 (MQ=255)                                          ggTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCt                           <  1:2791662/71‑1 (MQ=255)                                          ggTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCt                           <  1:2933592/71‑1 (MQ=255)                                          ggTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCt                           <  1:988609/71‑1 (MQ=255)                                            tGACGTCCTTCCTGTTGCGATGCTTGTGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTc                         <  1:616659/71‑1 (MQ=255)                                            tGACGTCCTTCCTGTTGCGATGCTTGTGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTc                         <  1:2084521/71‑1 (MQ=255)                                            tGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCt                           >  1:3422033/1‑69 (MQ=255)                                            tGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTc                         >  1:952760/1‑71 (MQ=255)                                            tGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGATTGGCATTc                         <  1:139397/71‑1 (MQ=255)                                               cGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGcc                      <  1:3062986/71‑1 (MQ=255)                                                  ccttcctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGcc                      >  1:3102884/1‑68 (MQ=255)                                                  ccttcctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGc                       >  1:3092855/1‑67 (MQ=255)                                                  ccttcctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATc                   >  1:2656319/1‑71 (MQ=255)                                                  ccttcctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATc                   >  1:2403853/1‑71 (MQ=255)                                                   cttcctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCCTcc                  >  1:3630273/1‑71 (MQ=255)                                                     tcctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCca                >  1:391585/1‑71 (MQ=255)                                                     tcctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCGTCGCCATCCca                >  1:2160932/1‑71 (MQ=255)                                                      cctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATccc                 <  1:2988006/69‑1 (MQ=255)                                                      cctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCca                >  1:3216193/1‑70 (MQ=255)                                                      cctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCca                >  1:1158353/1‑70 (MQ=255)                                                       ctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGc                                       >  1:1607310/1‑46 (MQ=255)                                                          ttGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACAcc            <  1:2501523/70‑1 (MQ=255)                                                          ttGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACAcc            <  1:1517660/70‑1 (MQ=255)                                                           tGCGATGCTTGCGCGATACGTTACTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCacac             >  1:1247123/1‑68 (MQ=255)                                                            gCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCacac             >  1:1844873/1‑67 (MQ=255)                                                                 gCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGc                             >  1:130527/1‑46 (MQ=255)                                                                  cTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaa    <  1:3188381/70‑1 (MQ=255)                                                                  cTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaa    <  1:1183043/70‑1 (MQ=255)                                                                  cTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaa    <  1:2715558/70‑1 (MQ=255)                                                                    tGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGtt                                 <  1:859072/39‑1 (MQ=255)                                                                     gcgcGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaata  <  1:2085080/69‑1 (MQ=255)                                                                      |                                                                    CGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATA  >  W3110S.gb/1016967‑1017103

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