Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 150,058 G→A 12.1% P298S (CCG→TCG)  PP_0142 ← ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947150,0580GA12.1% 78.1 / 3.4 33P298S (CCG→TCG) PP_0142ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (14/15);  new base A (2/2);  total (16/17)
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.75e-01

TCGGCGCTGCCGCTGACCTTGAAGCCTTCCAGGCAGCCGATGGCGGCGATCAGGAAGGCAAAGAACGGTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATCAGCAGGGCCAGTACCCGCGGCACCACCAGCAGCTCCATGGGGTTGAGGCCCAGCGTG  >  NC_002947/149919‑150203
                                                                                                                                           |                                                                                                                                                 
tCGGCGCTGCCGCTGACCTTGAAGCCTTCCAGGCAGCCGATGGCGGCGAGCAGGAAGGCAAAGAACGGTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATcc                                                                                                                                         >  2:751965/1‑150 (MQ=255)
         ccggttaCCTTGTAGCCTTCCAGTCAGCCGTTGGCGGCGATCAGGAAGTCAAAGAACGGTGCCTTGGCTAGGCCTACTCGGAAGTGCTGCACGCCAATGTCGCTTTCCAGCAGGGACAGTAACAGTGCCGGCGAAATATCCAGTGTCAgg                                                                                                                                <  2:71391/144‑1 (MQ=255)
          cGCTGACCTTGAAGCCTTCCAGGCAGCCGATGGCGGCGATCAGGAAGGCAAAGAACGGTGCCTTCCTTCTGTCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAggg                                                                                                                               <  1:49227/150‑1 (MQ=255)
           gCTGACCTTGAAGCCTTCCAGGCAGCCGATGGCGGCGATCAGGAAGGCAAAGAACGGTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAggca                                                                                                                              <  1:751965/150‑3 (MQ=255)
              gACCTTGAAGCCTTCCAGGCAGCCGATGGCGGCGATCAGGAAGGCAAAGAACGGTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCa                                                                                                                           <  2:670845/150‑1 (MQ=255)
                 cTTGAAGCCTTCCAGGCAGCCGATGGCGGCGATCAGGAAGGCAAAGAACGGTGCCTTGGCGAGTGCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGGCGGCGAAATATCCAAATGCAGGGCTCATac                                                                                                                        >  1:362858/1‑150 (MQ=255)
                 cTTGAAGCCTTCCAGGCAGCCGATGGCGGCGATCAGGAAGGCAAAGAACGGTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATac                                                                                                                        >  2:297490/1‑150 (MQ=255)
                      aGCCTTCCAGGCAGCCGATGGCGGCGATCAGGAAGGCAAAGAACGGTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCAccg                                                                                                                   >  1:217430/1‑150 (MQ=255)
                         cTTCCAGGCAGCCGATGGCGGCGATCAGGAAGGCAAAGAACGGTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCAccgccc                                                                                                                >  2:591266/1‑150 (MQ=255)
                                         ggcggcGATCAGGAAGGCAAAGAACGGTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCa                                                                                                <  2:217430/150‑1 (MQ=255)
                                                           aaaGAACGGTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGt                                                                              >  1:139345/1‑150 (MQ=255)
                                                                   gTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACg                                                                      <  2:139345/150‑1 (MQ=255)
                                                                   gTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACg                                                                      >  1:123352/1‑150 (MQ=255)
                                                                   gTGCATTGTCCAGGCCTCGCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACCTGGCCGGCGAAATCTCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACg                                                                      <  2:362858/150‑1 (MQ=255)
                                                                       cTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAg                                                                  >  2:518242/1‑150 (MQ=255)
                                                                             cAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAAt                                                            <  2:510998/150‑1 (MQ=255)
                                                                                 ccTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATcagc                                                        <  2:123352/150‑1 (MQ=255)
                                                                                  cTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGACGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGc                                                                                                    <  1:342500/105‑1 (MQ=255)
                                                                                  cTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGACGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGc                                                                                                    >  2:342500/1‑105 (MQ=255)
                                                                                  cTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGACGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGc                                                                                                    <  1:342450/105‑1 (MQ=255)
                                                                                  cTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGACGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGc                                                                                                    >  2:342450/1‑105 (MQ=255)
                                                                                      cAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGCGCGAAGACCACCGCCCCGCTAATGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATCAGCAGGGc                                                   <  1:297490/150‑1 (MQ=255)
                                                                                       aGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCaa                                                                                                          <  2:767236/94‑1 (MQ=255)
                                                                                       aGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCaa                                                                                                          >  1:767236/1‑94 (MQ=255)
                                                                                        ggAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGAt                                                                                                      >  2:231810/1‑97 (MQ=255)
                                                                                        ggAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGAt                                                                                                      <  1:231810/97‑1 (MQ=255)
                                                                                            gTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGTCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATCAGCAGGGCCAGTAc                                             <  2:52169/150‑1 (MQ=255)
                                                                                            gTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATCAGCAGGGCCAGTAc                                             <  1:591266/150‑1 (MQ=255)
                                                                                                       cAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAg                                                                                                                                       >  2:339860/1‑49 (MQ=255)
                                                                                                       cAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAg                                                                                                                                       <  1:339860/49‑1 (MQ=255)
                                                                                                                     agcagcGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATCAGCAGGGCCAGTACCCGCGGCACCACCAGCAGCTCCATg                    >  1:202767/1‑150 (MQ=255)
                                                                                                                                gAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATCAGCAGGGCCAGTACCCGCGGCACCACCAGCAGCTCCATGGGGTTGAgatc         >  2:181783/1‑147 (MQ=255)
                                                                                                                                       gCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATCAGCAGGGCCAGTACCCGCGGCACCACCAGCAGCTCCATGGGGTTGAGGCCCAGCGTg  >  1:694396/1‑150 (MQ=255)
                                                                                                                                          ggCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCAGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATCAGCAGGGCCAGTACCCGCGGCACCACCAGCAGCTCCATggg                  <  1:566398/131‑1 (MQ=255)
                                                                                                                                          ggCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCAGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATCAGCAGGGCCAGTACCCGCGGCACCACCAGCAGCTCCATggg                  >  2:566398/1‑131 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
TCGGCGCTGCCGCTGACCTTGAAGCCTTCCAGGCAGCCGATGGCGGCGATCAGGAAGGCAAAGAACGGTGCCTTGGCCAGGCCTACCAGGAAGTGCTGCACGCCAATGTCGCTTTGCAGCAGCGACAGGAACATGGCCGGCGAAATATCCAGTGTCAGGGCGCATACCACCGCCCCGCCAACGATGCCGCAGATCATGGCGACAAAGGTGAGCAACGGCAGCGAAATCAGCAGGGCCAGTACCCGCGGCACCACCAGCAGCTCCATGGGGTTGAGGCCCAGCGTG  >  NC_002947/149919‑150203

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