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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb237,7260AG14.9% 113.6 / 12.8 47N131S (AAC→AGC) yafTpredicted aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (16/24);  new base G (2/5);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGA  >  W3110S.gb/237660‑237793
                                                                  |                                                                   
caTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 <  1:283958/71‑1 (MQ=255)
caTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 <  1:2806515/71‑1 (MQ=255)
caTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 <  1:1844808/71‑1 (MQ=255)
 aTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAgg                                                                <  1:2804044/71‑1 (MQ=255)
 aTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 >  1:916376/1‑70 (MQ=255)
 aTTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 >  1:1209769/1‑70 (MQ=255)
   tctcTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGa                                                               <  1:1160097/70‑1 (MQ=255)
         gTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATGTTGc                                                         >  1:2428102/1‑70 (MQ=255)
         gTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGcc                                                        <  1:816345/71‑1 (MQ=255)
         gTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGcc                                                        <  1:2297572/71‑1 (MQ=255)
         gTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGc                                                         <  1:1965124/70‑1 (MQ=255)
              ccGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCa                                                       >  1:2207219/1‑67 (MQ=255)
              ccGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTg                                                    >  1:2779688/1‑70 (MQ=255)
                    aGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGaa                                             <  1:1545747/71‑1 (MQ=255)
                      gCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAg                                                      <  1:957867/60‑1 (MQ=255)
                      gCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAata                                           >  1:1045636/1‑71 (MQ=255)
                       cAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAg                                                      >  1:123968/1‑59 (MQ=255)
                       cAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAatat                                          >  1:2734779/1‑71 (MQ=255)
                        aaaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAatat                                          >  1:1180324/1‑70 (MQ=255)
                        aaaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATg                                         <  1:2504538/71‑1 (MQ=255)
                         aaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATGTTGCCAGTGCTGAGaa                                             <  1:2581505/66‑1 (MQ=255)
                         aaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATGTTGCCAGTGCTGAGaa                                             <  1:1859061/66‑1 (MQ=255)
                         aaaGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGa                                        <  1:1926319/71‑1 (MQ=255)
                             ccTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:1278605/71‑1 (MQ=255)
                               tctGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAg                                                                 <  1:320647/40‑1 (MQ=255)
                                 tGACTTGTAACAAAAACAAAGATAAAAATGAGAGCGAGGATGTTGCCAGTGCTGAGAATATGAAc                                      >  1:643104/1‑65 (MQ=255)
                                 tGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGtt                                                           >  1:1325297/1‑44 (MQ=255)
                                 tGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATg                                >  1:2150688/1‑71 (MQ=255)
                                 tGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATg                                >  1:2801193/1‑71 (MQ=255)
                                 tGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATg                                >  1:1368582/1‑71 (MQ=255)
                                     ttGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATg                                         <  1:65545/58‑1 (MQ=255)
                                       gTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAgg                            <  1:2182565/69‑1 (MQ=255)
                                       gTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAgg                            <  1:193325/69‑1 (MQ=255)
                                       gTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGACAGTg                                                    >  1:540107/1‑45 (MQ=255)
                                          acaaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGc                           <  1:1276810/67‑1 (MQ=255)
                                           caaaaacaaaGATAAAAATGAGAGCGAGGATGTTGCCAGTGCTGAGAATATg                                         <  1:738076/52‑1 (MQ=255)
                                           caaaaacaaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGaa                                       >  1:1972718/1‑54 (MQ=255)
                                                 caaaGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGa                                  >  1:2054235/1‑53 (MQ=255)
                                                       tAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATg                                         <  1:1325750/40‑1 (MQ=255)
                                                         aaaaTGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAAc                                      >  1:240954/1‑41 (MQ=255)
                                                          aaaTGAGAGCGAGGATGTTGCCACTGCTGAGAatatt                                         <  1:1930755/37‑2 (MQ=37)
                                                           aaTGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:2658736/41‑1 (MQ=255)
                                                           aaTGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:2537600/41‑1 (MQ=255)
                                                             tGAGAGCGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAg    <  1:1433266/71‑1 (MQ=255)
                                                             tGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:883633/39‑1 (MQ=255)
                                                             tGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTg                                    <  1:2060600/39‑1 (MQ=255)
                                                               agaACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGa  <  1:1312353/71‑1 (MQ=255)
                                                                 aaCGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAgg                            <  1:30972/43‑1 (MQ=255)
                                                                  |                                                                   
CATTCTCTGGTATCCCGGACAGGCAAAAGCCTCTGACTTGTAACAAAAACAAAGATAAAAATGAGAACGAGGATGTTGCCAGTGCTGAGAATATGAACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGA  >  W3110S.gb/237660‑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: