Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	237,868	A→G	8.5%	V178V (GTA→GTG)	yafT →	predicted aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	237,868	0	A	G	8.5%	226.5 / 11.8	82	V178V (GTA→GTG)	yafT	predicted aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (29/46);  new base G (4/3);  total (33/49)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

AGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTG  >  W3110S.gb/237798‑237937                                                                       |                                                                      aGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTa                                                                       <  1:1164092/71‑1 (MQ=255)  gAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTAt                                                                      <  1:3392195/71‑1 (MQ=255)       cGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTcc                                                                 <  1:1314427/71‑1 (MQ=255)           tattGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTcc                                                                 <  1:2654111/67‑1 (MQ=255)                tttACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATaaa                                                        <  1:828877/71‑1 (MQ=255)                   aCCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATaaaaa                                                      <  1:3588533/70‑1 (MQ=255)                    ccAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAAt                                                     <  1:783512/70‑1 (MQ=255)                     cAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGa                                                   >  1:1154546/1‑71 (MQ=255)                     cAGGATATGTTGCAGACAGGAAAATATGACTATGCGCTGAAATCAACAGTATGGTcc                                                                 <  1:637513/57‑1 (MQ=255)                      aGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTcc                                                                 <  1:1575858/56‑1 (MQ=255)                      aGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGAc                                                  <  1:3211238/71‑1 (MQ=255)                       ggATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGAc                                                  >  1:1286244/1‑70 (MQ=255)                        gATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATaaaaa                                                      >  1:854330/1‑65 (MQ=255)                        gATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACa                                                 >  1:241376/1‑70 (MQ=255)                         atatGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACaaa                                               >  1:5010177/1‑71 (MQ=255)                          tatGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTGTGGTCCGACTATAAAAATGACATAc                                              >  1:3551418/1‑71 (MQ=255)                          tatGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTGTGGTCCGACTATAAAAATGACAAAc                                              >  1:626220/1‑71 (MQ=255)                            tGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTGTGGTCCGACTATAAAAATGACAAACt                                             <  1:1481463/70‑1 (MQ=255)                            tGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTGTGGTCCGACTATAAAAATGACAAACt                                             <  1:2184710/70‑1 (MQ=255)                            tGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTcc                                                                 <  1:4866253/50‑1 (MQ=255)                            tGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTcc                                                                 <  1:191635/50‑1 (MQ=255)                            tGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACt                                             <  1:667677/70‑1 (MQ=255)                            tGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTc                                            <  1:5120159/71‑1 (MQ=255)                            tGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTc                                            <  1:3311637/71‑1 (MQ=255)                             gTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATaaa                                                        <  1:2515679/58‑1 (MQ=255)                             gTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACaaa                                               >  1:257798/1‑67 (MQ=255)                              ttGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTGTGGTCCGACTATaaa                                                        >  1:2762672/1‑57 (MQ=255)                              ttGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCAc                                          >  1:2740558/1‑71 (MQ=255)                              ttGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCAc                                          >  1:3596294/1‑71 (MQ=255)                              ttGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCAc                                          >  1:675771/1‑71 (MQ=255)                                gCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTGTGGTCCGACTATAAAAATGACAAACTCACt                                         >  1:2021303/1‑70 (MQ=255)                                gCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATaaa                                                        <  1:2783655/55‑1 (MQ=255)                                gCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACt                                         <  1:401566/70‑1 (MQ=255)                                gCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACt                                         <  1:1664007/70‑1 (MQ=255)                                   gACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTcc                                                                 <  1:2875390/43‑1 (MQ=255)                                    aCAGGAAAATATGACTCTGCGCTGAAATCAACAGTGTGGTCCGACTATAAAAATGACAAACTCACTGACGc                                    <  1:1925619/71‑1 (MQ=255)                                     cAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGa                                                               >  1:4348547/1‑43 (MQ=255)                                     cAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACg                                     >  1:3570926/1‑69 (MQ=255)                                     cAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGc                                    >  1:1446367/1‑70 (MQ=255)                                     cAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCt                                   <  1:3679368/71‑1 (MQ=255)                                     cAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCt                                   <  1:4903621/71‑1 (MQ=255)                                     cAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCt                                   <  1:194101/71‑1 (MQ=255)                                       ggAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATaaaa                                                       <  1:4452514/49‑1 (MQ=255)                                        gAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGAc                                                  <  1:1879384/53‑1 (MQ=255)                                        gAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACt                                             <  1:3613043/58‑1 (MQ=255)                                         aaaaTATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCt                                   <  1:4170634/67‑1 (MQ=255)                                           aaTATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTAt                                 >  1:4973605/1‑67 (MQ=255)                                           aaTATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTcc                             <  1:3275050/71‑1 (MQ=255)                                              atGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACaaa                                               >  1:1028512/1‑50 (MQ=255)                                               tGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCa                                           >  1:3699421/1‑53 (MQ=255)                                               tGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCt                           >  1:599677/1‑69 (MQ=255)                                               tGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCt                           >  1:2236732/1‑69 (MQ=255)                                               tGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCt                           >  1:179882/1‑69 (MQ=255)                                               tGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTgc                         >  1:4731653/1‑71 (MQ=255)                                               tGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTg                          <  1:2879956/70‑1 (MQ=255)                                               tGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTg                          <  1:3063268/70‑1 (MQ=255)                                               tGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTg                          <  1:3337637/70‑1 (MQ=255)                                               tGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTg                          <  1:3897701/70‑1 (MQ=255)                                                gACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTgcg                        >  1:3449923/1‑71 (MQ=255)                                                 aCTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTgcgc                       <  1:3925161/71‑1 (MQ=255)                                                  ctctGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCt                      <  1:4079097/71‑1 (MQ=255)                                                   tctGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACt                                         <  1:2117461/51‑1 (MQ=255)                                                    ctGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACaaa                                               >  1:3842390/1‑44 (MQ=255)                                                    ctGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACaaa                                               >  1:2324398/1‑44 (MQ=255)                                                     tGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTgcgc                       >  1:2509982/1‑67 (MQ=255)                                                     tGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTAc                    >  1:561044/1‑70 (MQ=255)                                                        gcTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTgg                <  1:1887704/71‑1 (MQ=255)                                                          tGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACt                                             <  1:3312298/40‑1 (MQ=255)                                                          tGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCt                                   <  1:1090328/50‑1 (MQ=255)                                                          tGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTa              >  1:4740555/1‑71 (MQ=255)                                                           gAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCt                                   <  1:3455174/49‑1 (MQ=255)                                                           gAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTAc             >  1:316385/1‑71 (MQ=255)                                                           gAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTAc             >  1:2511234/1‑71 (MQ=255)                                                            aaaTCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACg            >  1:3791200/1‑71 (MQ=255)                                                            aaaTCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACg            <  1:1105899/71‑1 (MQ=255)                                                               tCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGttt         <  1:2559336/71‑1 (MQ=255)                                                               tCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGttt         <  1:943663/71‑1 (MQ=255)                                                                 aaCAGTATGGTCCGACTATAAAAATGACAAACTCACTGa                                       >  1:2860254/1‑39 (MQ=255)                                                                   cAGTATTGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACg     <  1:966008/71‑1 (MQ=255)                                                                     gTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCt   <  1:3336467/71‑1 (MQ=255)                                                                     gTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCt   <  1:1676638/71‑1 (MQ=255)                                                                     gTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCt   <  1:1158149/71‑1 (MQ=255)                                                                      tATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCtg  <  1:2156772/71‑1 (MQ=255)                                                                       aTGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTcc                             <  1:3289588/43‑1 (MQ=255)                                                                       |                                                                      AGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTG  >  W3110S.gb/237798‑237937

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