Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,355,187 A→G 14.8% V95A (GTG→GCG) ‡ sapB ← antimicrobial peptide transport ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,355,1870AG14.8% 134.4 / 14.2 53V95A (GTG→GCG) ‡sapBantimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (18/27);  new base G (3/5);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

CGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAA  >  NC_000913/1355044‑1355335
                                                                                                                                               |                                                                                                                                                    
cGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACgggg                                                                                                                                                 <  2:293816/149‑1 (MQ=255)
   gCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATg                                                                                                                                              >  2:46724/1‑149 (MQ=255)
    caaccaacCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCGTGGGGATga                                                                                                                                             <  3:139697/149‑2 (MQ=255)
     aaccaaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCGTGGGGATgat                                                                                                                                            <  1:175083/149‑3 (MQ=255)
      accaaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCca                                                                                                                                           >  4:331679/1‑149 (MQ=255)
        caaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCcaac                                                                                                                                         >  2:35635/1‑149 (MQ=255)
        caaccTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCcaac                                                                                                                                         >  2:108667/1‑149 (MQ=255)
                          cAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGaa                                                                                                                       >  4:158142/1‑149 (MQ=255)
                           aaaaCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAAt                                                                                                                                      <  1:332770/133‑1 (MQ=255)
                           aaaaCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAAt                                                                                                                                      >  2:332770/1‑133 (MQ=255)
                            aaaCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAgg                                                                                                                     >  2:14400/1‑149 (MQ=255)
                              aCAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCg                                                                                                                   >  4:17621/1‑149 (MQ=255)
                               taGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGa                                                                                                                  <  3:38318/148‑1 (MQ=255)
                                    tCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATg                                                                                                             >  3:146484/1‑149 (MQ=255)
                                      aaCAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGca                                                                                                           >  4:94836/1‑149 (MQ=255)
                                      aaCAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGca                                                                                                           <  1:306835/149‑1 (MQ=255)
                                       aCAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGcac                                                                                                          >  3:139339/1‑149 (MQ=255)
                                           aaGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGtt                                                                                                      <  2:346737/149‑1 (MQ=255)
                                              cgcCAGCCAGAACTCCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCACTAATCCCCACGTGGATGCCAACTATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCa                                                                                                   <  2:35214/149‑1 (MQ=255)
                                              cgcCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCa                                                                                                   <  2:163178/149‑1 (MQ=255)
                                               gccagccaGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAt                                                                                                  <  3:158142/149‑1 (MQ=255)
                                                ccagccagAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAtc                                                                                                 <  1:59200/149‑1 (MQ=255)
                                                             cGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGa                                                                                    <  3:257230/149‑1 (MQ=255)
                                                              gggATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCGCCGTGGGGATga                                                                                                                                             >  4:343021/1‑90 (MQ=255)
                                                              gggATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCGTGGGGATga                                                                                                                                             <  3:343021/91‑2 (MQ=255)
                                                              gggATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGAc                                                                                   <  2:59393/149‑1 (MQ=255)
                                                                   tGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtt                                                                              <  3:17621/149‑1 (MQ=255)
                                                                     aaaaGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtttt                                                                            <  3:331679/149‑1 (MQ=255)
                                                                     aaaaGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtttt                                                                            <  1:46724/149‑1 (MQ=255)
                                                                         gCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACt                                                                        <  4:146484/149‑1 (MQ=255)
                                                                           ccAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTaa                                                                          >  4:288384/1‑145 (MQ=255)
                                                                           ccAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTaa                                                                          <  3:288384/145‑1 (MQ=255)
                                                                           ccAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGt                                                                      >  2:392648/1‑149 (MQ=255)
                                                                            cAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGa                                                                                                               <  2:160755/107‑1 (MQ=255)
                                                                            cAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGa                                                                                                               >  1:160755/1‑107 (MQ=255)
                                                                            cAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAg                                                                                                        <  4:396571/114‑1 (MQ=255)
                                                                            cAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAg                                                                                                        >  3:396571/1‑114 (MQ=255)
                                                                            cAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGtt                                                                     <  2:255808/149‑1 (MQ=255)
                                                                                              gATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGt                                                   <  2:365330/149‑1 (MQ=255)
                                                                                                  aaGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAAt                                               <  1:392648/149‑1 (MQ=255)
                                                                                                       aTCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCaaaa                                                                                                                            >  1:402662/1‑67 (MQ=255)
                                                                                                       aTCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCaaaa                                                                                                                            <  2:402662/67‑1 (MQ=255)
                                                                                                          cTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGaaa                                       <  3:94836/149‑1 (MQ=255)
                                                                                                                          gTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTg                       <  1:108667/149‑1 (MQ=255)
                                                                                                                          gTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTg                       <  1:35635/149‑1 (MQ=255)
                                                                                                                              tCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGtt                                                  <  1:37564/118‑1 (MQ=255)
                                                                                                                              tCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATc                   <  4:139339/149‑1 (MQ=255)
                                                                                                                              tCCCGGCAATCATCCCCACGGGGAGGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGtt                                                  >  2:37564/1‑118 (MQ=255)
                                                                                                                                  tgttggcatcCCCGTGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCt                                            <  3:56590/114‑1 (MQ=255)
                                                                                                                                  tgttggcatcCCCGTGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCt                                            >  4:56590/7‑120 (MQ=255)
                                                                                                                                    ttggcatcCCCGTGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCg             >  3:251973/5‑149 (MQ=255)
                                                                                                                                         atcCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTaaa        >  1:263680/1‑149 (MQ=255)
                                                                                                                                          tcCCCGCGTGGATGCCAACGATCAGCGCAAAACCGAAGGCGAGGATGCACTGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAAc       <  2:263680/149‑1 (MQ=255)
                                                                                                                                               aCGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGaa  >  3:373951/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTTGACCGTTAATGCTGGAAACACCAAAATCCCAGTGGATCAGGCCGTTAAACCAGAA  >  NC_000913/1355044‑1355335

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