Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,355,756 G→A 8.9% S159S (TCC→TCT sapC ← predicted antimicrobial peptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,355,7560GA 8.9% 126.0 / 2.6 45S159S (TCC→TCTsapCpredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (16/25);  new base A (2/2);  total (18/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 = 8.09e-01

ATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTAT  >  W3110S.gb/1355687‑1355826
                                                                     |                                                                      
aTGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTgg                                                                        <  1:864607/70‑1 (MQ=255)
   cACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACa                                                                    <  1:3314402/71‑1 (MQ=255)
    aCCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACaa                                                                   <  1:2256321/71‑1 (MQ=255)
    aCCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACa                                                                    <  1:1044130/70‑1 (MQ=255)
    aCCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACa                                                                    <  1:711415/70‑1 (MQ=255)
       aTGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATGGCGTgaga                                                                      >  1:1930072/1‑65 (MQ=255)
       aTGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATGGCGTGAGACAAACt                                                                >  1:777629/1‑71 (MQ=255)
       aTGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACt                                                                <  1:923855/71‑1 (MQ=255)
       aTGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACt                                                                <  1:3126739/71‑1 (MQ=255)
          cTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCgg                                                             >  1:2904114/1‑71 (MQ=255)
             tAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTc                                                               <  1:1498566/66‑1 (MQ=255)
               aaTCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCgg                                                        >  1:841247/1‑71 (MQ=255)
               aaTCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCgg                                                        >  1:475263/1‑71 (MQ=255)
                aTCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGc                                                       <  1:170402/71‑1 (MQ=255)
                 tCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGc                                                       >  1:3160975/1‑70 (MQ=255)
                   gAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATGGCGTGAGACAAACTc                                                               <  1:1163952/60‑1 (MQ=255)
                   gAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATGGCGTGAGACAAACTc                                                               <  1:1630742/60‑1 (MQ=255)
                   gAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCaaa                                                    >  1:1302168/1‑71 (MQ=255)
                             tACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTc                                                               <  1:1835138/50‑1 (MQ=255)
                               tgcgGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCgg                                                             >  1:2933303/2‑50 (MQ=255)
                               cgcgGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAAc                                                                 >  1:2988024/1‑46 (MQ=255)
                               cgcgGCAGCAATGCCAGCCAGACAGAAAACATAGCGTGGGACAAACTCGGTCCGGCaaa                                                    >  1:207468/1‑59 (MQ=255)
                                gcgGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATa                                       <  1:668995/71‑1 (MQ=255)
                                gcgGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATa                                       <  1:3477222/71‑1 (MQ=255)
                                  ggCAGCAATGCCAGCCAGACAGCAAACATGGCGTGAGACAAACTCGGTCCGGCAAACGCTACaa                                            >  1:2891542/1‑64 (MQ=255)
                                   gcagcaATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACt                                                                <  1:1362574/43‑1 (MQ=255)
                                   gcagcaATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACt                                                                <  1:3471013/43‑1 (MQ=255)
                                   gcagcaATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAAtt                                    <  1:3756111/71‑1 (MQ=255)
                                     agcaATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCg                                                              >  1:1326411/1‑43 (MQ=255)
                                     agcaATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATaa                                      >  1:2826319/1‑67 (MQ=255)
                                       caATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCa                                                      >  1:2524217/1‑49 (MQ=255)
                                          tGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCaa                               >  1:2004506/1‑69 (MQ=255)
                                           gccagccaGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTAc                                              <  1:661838/53‑1 (MQ=255)
                                            ccagccagACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAAcagc                           <  1:3678112/71‑1 (MQ=255)
                                                ccagACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCg                       <  1:1552356/71‑1 (MQ=255)
                                                 cagACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAAca                             >  1:3181833/1‑64 (MQ=255)
                                                 cagACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAAca                             >  1:1747729/1‑64 (MQ=255)
                                                 cagACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCg                       <  1:1556865/70‑1 (MQ=255)
                                                  agACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGaa                     <  1:1607551/71‑1 (MQ=255)
                                                     cAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGAAgg                   >  1:1500965/1‑70 (MQ=255)
                                                     cAGCAAACATAGCGTGGGACAAACTCGGTACGGCAAACGCTACAACGATAATTGCCAACAGCAGCGAAggg                  <  1:2859245/71‑1 (MQ=255)
                                                       gCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGAAGGGa                 <  1:21265/70‑1 (MQ=255)
                                                       gCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGAAGGGAt                >  1:866718/1‑71 (MQ=255)
                                                          aaCTTAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGaa                     <  1:1663271/63‑1 (MQ=255)
                                                          aaCATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGaa                     <  1:75083/63‑1 (MQ=255)
                                                          aaCATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGAACGGATCGc             <  1:1739556/71‑1 (MQ=255)
                                                             aTAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCg                       <  1:2815347/58‑1 (MQ=255)
                                                                     ggACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTAt  >  1:3820197/1‑71 (MQ=255)
                                                                     |                                                                      
ATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTAT  >  W3110S.gb/1355687‑1355826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: