Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	237,726	A→G	8.2%	N131S (AAC→AGC)	yafT →	predicted aminopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	237,726	0	A	G	8.2%	158.1 / 6.3	62	N131S (AAC→AGC)	yafT	predicted aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (31/25);  new base G (3/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

TCCACATTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGA  >  W3110S.gb/237656‑237793                                                                       |                                                                    tCCACATTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGaa                                                                     >  1:4559226/1‑71 (MQ=255)   cacaTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAAc                                                                    >  1:3873360/1‑70 (MQ=255)     caTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 <  1:154967/71‑1 (MQ=255)     caTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 <  1:1643506/71‑1 (MQ=255)       ttCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGa                                                               >  1:2811350/1‑71 (MQ=255)         ctctGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATg                                                             >  1:3275041/1‑71 (MQ=255)         ctctGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATg                                                             >  1:3916930/1‑71 (MQ=255)              gTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTg                                                          <  1:1864802/69‑1 (MQ=255)              gTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGc                                                         <  1:67610/70‑1 (MQ=255)              gTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGc                                                         <  1:2751433/70‑1 (MQ=255)                 tCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAg                                                      <  1:1997036/70‑1 (MQ=255)                   ccGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTg                                                    >  1:495040/1‑70 (MQ=255)                   ccGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGc                                                   >  1:4206598/1‑71 (MQ=255)                    cGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGtt                                                           <  1:2408496/62‑1 (MQ=255)                    cGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCt                                                  >  1:660030/1‑71 (MQ=255)                       aCAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTgag                                               <  1:1832557/71‑1 (MQ=255)                       aCAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTgag                                               <  1:1057073/71‑1 (MQ=255)                       aCAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTga                                                <  1:3763344/70‑1 (MQ=255)                         aGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATGTTGCCAGTGCTGAGaa                                             <  1:2585357/71‑1 (MQ=255)                           gCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAata                                           >  1:2458080/1‑71 (MQ=255)                            cAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATGTTGCCAGTGCTGAGAatat                                          >  1:487922/1‑71 (MQ=255)                            cAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATGTTGCCAGTGCTGAGAatat                                          >  1:2653441/1‑71 (MQ=255)                                aGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATGTTGCCAGTGCTGAGAATATGAAc                                      >  1:1864894/1‑71 (MQ=255)                                aGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAAc                                      >  1:2422874/1‑71 (MQ=255)                                 gCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACt                                     >  1:4070609/1‑71 (MQ=255)                                  ccTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:2519464/71‑1 (MQ=255)                                    tctGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGaa                                             <  1:1273578/60‑1 (MQ=255)                                      tGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGa                                                               >  1:4021133/1‑40 (MQ=255)                                      tGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGa                                                               >  1:2923257/1‑40 (MQ=255)                                      tGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATg                                >  1:3964154/1‑71 (MQ=255)                                        aCTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGAt                                                              >  1:1187817/1‑39 (MQ=255)                                        aCTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATg                                >  1:438489/1‑69 (MQ=255)                                        aCTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCa                              <  1:2216754/71‑1 (MQ=255)                                        aCTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCa                              <  1:530474/71‑1 (MQ=255)                                          ttGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGa                                        <  1:3356738/59‑1 (MQ=255)                                            gTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTg                                                          >  1:1932712/1‑39 (MQ=255)                                            gTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATg                                         <  1:530343/56‑1 (MQ=255)                                            gTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCa                          >  1:1121834/1‑71 (MQ=255)                                               acaaaaacaaaGATAAAAATGAGAGCGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTg                       <  1:2604464/71‑1 (MQ=255)                                               acaaaaacaaaGATAAAAATGAGAACGAGGATGTTGc                                                         <  1:445278/37‑1 (MQ=255)                                               acaaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAata                                           >  1:1091016/1‑51 (MQ=255)                                               acaaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGaa                                       >  1:1556132/1‑55 (MQ=255)                                               acaaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACt                        >  1:2816793/1‑70 (MQ=255)                                               acaaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTg                       >  1:2701843/1‑71 (MQ=255)                                                 aaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTgaga                                              <  1:251006/46‑1 (MQ=255)                                                 aaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAAc                                      >  1:3983850/1‑54 (MQ=255)                                                  aaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGaa                                       >  1:2093920/1‑52 (MQ=255)                                                  aaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCCGGCACTGCGt                    >  1:327098/1‑71 (MQ=255)                                                  aaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGt                    <  1:3438066/71‑1 (MQ=255)                                                   aaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAat                                            >  1:1528622/1‑46 (MQ=255)                                                   aaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATg                                         <  1:1931205/49‑1 (MQ=255)                                                    aacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGaa                                       >  1:812084/1‑50 (MQ=255)                                                     acaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGtttt                 <  1:672534/71‑1 (MQ=255)                                                      caaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCa                          >  1:4093120/1‑61 (MQ=255)                                                      caaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCa                          >  1:1487881/1‑61 (MQ=255)                                                         aGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTgg             >  1:503473/1‑71 (MQ=255)                                                           aTAAAAATGAGACCGAGGATGTTGCCAGTGCTGAGAat                                            >  1:2825848/1‑38 (MQ=255)                                                           aTAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAgg                            <  1:2981045/54‑1 (MQ=255)                                                             aaaaaTGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCa                          >  1:2424467/1‑54 (MQ=255)                                                                  tGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAg    <  1:2590014/71‑1 (MQ=255)                                                                   gagaACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAg    >  1:522866/1‑70 (MQ=255)                                                                    agaACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAgg   <  1:1224365/70‑1 (MQ=255)                                                                    agaACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGa  <  1:4563278/71‑1 (MQ=255)                                                                       |                                                                    TCCACATTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGA  >  W3110S.gb/237656‑237793

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