Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,355,750 A→G 10.6% A161A (GCT→GCC 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,7500AG10.6% 120.1 / 5.7 47A161A (GCT→GCCsapCpredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (16/26);  new base G (3/2);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGAAGGGATCGC  >  W3110S.gb/1355687‑1355815
                                                               |                                                                 
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)
                                                               |                                                                 
ATGCACCATGCTGTAAATCGAACGTACCATACGCGGCAGCAATGCCAGCCAGACAGCAAACATAGCGTGGGACAAACTCGGTCCGGCAAACGCTACAACGATAATTGCCAACAGCAGCGAAGGGATCGC  >  W3110S.gb/1355687‑1355815

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: