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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb237,7260AG15.1% 130.9 / 16.1 53N131S (AAC→AGC) yafTpredicted aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (15/30);  new base G (3/5);  total (18/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

TCCACATTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGA  >  W3110S.gb/237656‑237793
                                                                      |                                                                   
tCCACATTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGaa                                                                     >  1:2254214/1‑71 (MQ=255)
tCCACATTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGaa                                                                     >  1:2823377/1‑71 (MQ=255)
  cacaTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAAc                                                                    <  1:1671658/70‑1 (MQ=255)
  cacaTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAAc                                                                    >  1:2302119/1‑70 (MQ=255)
      ttCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGa                                                               >  1:378731/1‑71 (MQ=255)
        ctctGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATg                                                             >  1:3202735/1‑71 (MQ=255)
             gTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGcc                                                        <  1:2168022/71‑1 (MQ=255)
             gTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGcc                                                        <  1:3688723/71‑1 (MQ=255)
               aTCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAg                                                      >  1:1650108/1‑71 (MQ=255)
                tCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGa                                                               >  1:1794973/1‑61 (MQ=255)
                tCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGa                                                               <  1:3911819/61‑1 (MQ=255)
                  ccGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGaa                                                                     >  1:1914132/1‑53 (MQ=255)
                     gACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGa                                                               <  1:3269600/56‑1 (MQ=255)
                           cAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 <  1:2812726/48‑1 (MQ=255)
                           cAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 <  1:3898115/48‑1 (MQ=255)
                            aaaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATg                                         <  1:773735/71‑1 (MQ=255)
                            aaaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATg                                         <  1:3057957/71‑1 (MQ=255)
                             aaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAg                                                                 <  1:2468712/46‑1 (MQ=255)
                             aaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAg                                                                 <  1:2548850/46‑1 (MQ=255)
                             aaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAg                                                                 <  1:3060550/46‑1 (MQ=255)
                             aaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGa                                        <  1:1315300/71‑1 (MQ=255)
                                gCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATg                                                             >  1:1888453/1‑47 (MQ=255)
                                gCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACt                                     >  1:3432119/1‑71 (MQ=255)
                                 ccTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:2509125/71‑1 (MQ=255)
                                   tctGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAg                                                                 <  1:1102464/40‑1 (MQ=255)
                                   tctGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAg                                                                 <  1:2081320/40‑1 (MQ=255)
                                   tctGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGGTGAGaa                                             <  1:4865/60‑1 (MQ=255)
                                   tctGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGaa                                             <  1:2733830/60‑1 (MQ=255)
                                   tctGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGaa                                             <  1:10099/60‑1 (MQ=255)
                                   tctGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGaa                                             <  1:4062042/60‑1 (MQ=255)
                                   tctGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGaa                                             <  1:2172397/60‑1 (MQ=255)
                                   tctGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGaa                                             <  1:2279721/60‑1 (MQ=255)
                                   tctGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGaa                                             <  1:2377353/60‑1 (MQ=255)
                                   tctGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGaa                                             <  1:2918660/60‑1 (MQ=255)
                                     tGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 >  1:1152407/1‑38 (MQ=255)
                                     tGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATg                                >  1:3970633/1‑71 (MQ=255)
                                       aCTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATGTTGCCAGTGCTGAGAATATGAACTGGa                                  >  1:360891/1‑67 (MQ=255)
                                       aCTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATg                                                             >  1:418927/1‑40 (MQ=255)
                                       aCTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:4055585/65‑1 (MQ=255)
                                       aCTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:112929/65‑1 (MQ=255)
                                       aCTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:2113808/65‑1 (MQ=255)
                                       aCTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:1260321/65‑1 (MQ=255)
                                           gTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCa                          >  1:3811652/1‑71 (MQ=255)
                                              acaaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTg                       >  1:3138446/1‑71 (MQ=255)
                                              acaaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTg                       >  1:308638/1‑71 (MQ=255)
                                              acaaaaacaaaGATAAAAATGACAACGAGGATGTTGca                                                        >  1:2654417/1‑37 (MQ=39)
                                                aaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCa                              <  1:2283303/62‑1 (MQ=255)
                                                aaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCa                              <  1:696829/62‑1 (MQ=255)
                                                  aaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATg                                         <  1:802690/49‑1 (MQ=255)
                                                  aaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATg                                         <  1:845163/49‑1 (MQ=255)
                                                   aacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCa                              >  1:317780/1‑59 (MQ=255)
                                                         gATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGa                                  >  1:788476/1‑49 (MQ=255)
                                                         gATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGa                                  >  1:801163/1‑49 (MQ=255)
                                                                 tGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAg    <  1:2532883/71‑1 (MQ=255)
                                                                   agaACGAGGATGTTGCCAGAGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGa  <  1:3305234/71‑1 (MQ=255)
                                                                    gaACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:2285944/36‑1 (MQ=255)
                                                                      |                                                                   
TCCACATTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGA  >  W3110S.gb/237656‑237793

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: