Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 764,126 T→C 9.7% Q133R (CAG→CGG)  PP_0656 ← putative Amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947764,1260TC 9.7% 104.5 / 3.0 41Q133R (CAG→CGG) PP_0656putative Amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base C (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

GCATGCCCAGCGCCAATGCTGCTTCACGTTGCCCACGGGCCACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAGCGCGGCGGCAAAGCCGAGCAGCACGG  >  NC_002947/763983‑764267
                                                                                                                                               |                                                                                                                                             
gcATGCCCAGCGCCAATGCTGCTTCACGTTGCCCACGGGCCACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACCAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCCCCTGCGGCa                                                                                                                                         >  1:544248/1‑150 (MQ=255)
         gcgcCAATGCTGCTTCACGTTGCCCACGGGCCACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGt                                                                                                                                <  1:937655/150‑1 (MQ=255)
           gcCAATGCTGCTTCACGTTGCCCACGGGCCACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAg                                                                                                                              <  2:55723/150‑1 (MQ=255)
             cAATGCTGCTTCACGTTGCCCACGGGCCACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGAt                                                                                                                            <  1:893395/150‑1 (MQ=255)
                  ctgctTCACGTTGCCCACGGGCCACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATcagca                                                                                                                       <  1:626117/150‑1 (MQ=255)
                         aCGTTGCCCACGGGCCACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTg                                                                                                                <  1:498816/150‑1 (MQ=255)
                            ttGCCCACGGGCCACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGAt                                                                                                             <  1:24950/150‑1 (MQ=255)
                                ccACGGGCCACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATcagc                                                                                                         >  1:833833/1‑150 (MQ=255)
                                       ccACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAgcggcg                                                                                                  >  1:1053339/1‑150 (MQ=255)
                                         aCGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTg                                                                                                >  2:387554/1‑150 (MQ=255)
                                             cGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCAc                                                                                            >  2:953481/1‑150 (MQ=255)
                                                     ccGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAAc                                                                                    <  2:389553/150‑1 (MQ=255)
                                                        gcgcGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGCGTTTGGATCAGCAGCGGCGTGCCACGGAAGAACGAg                                                                                 <  1:641242/150‑1 (MQ=255)
                                                        gcgcGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAg                                                                                 <  2:1076220/150‑1 (MQ=255)
                                                                    tcGCTCAGGAACGCCCCTTAGGCCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTg                                                                     >  2:976145/1‑150 (MQ=255)
                                                                        tCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCca                                                                 >  1:200832/1‑150 (MQ=255)
                                                                                gCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCAcacccacc                                                         <  2:384844/150‑1 (MQ=255)
                                                                                 ccccGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACccacca                                                        <  2:833833/150‑1 (MQ=255)
                                                                                 ccccGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACccacca                                                        <  2:1053339/150‑1 (MQ=255)
                                                                                  cccGTAGTTCAGCGACAAGGCAACGATGCCGGCGCTGACGGCCCCGCGTACCAGGCCCACATACGGCAGCCCCAGGTAGATCAGCAGCACCCGGATCAGCACCGGCCTCCCACGGCAGAACGACGCGCACAAGCTGGCCCCACcccccaa                                                       >  2:842677/1‑149 (MQ=255)
                                                                                   ccGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGc                                                      >  1:718065/1‑150 (MQ=255)
                                                                                       aGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTg                                                                     >  2:999322/1‑131 (MQ=255)
                                                                                       aGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTg                                                                     <  1:999322/131‑1 (MQ=255)
                                                                                              ggACAAGTCAATGATGCCGGCGCTGATGGCGCAGGGTACCAGGCTCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGtt                                           <  2:233000/149‑1 (MQ=255)
                                                                                               gACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTg                                          >  1:151405/1‑150 (MQ=255)
                                                                                                aCAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTgg                                         <  2:227031/150‑1 (MQ=255)
                                                                                                    ggCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATaa                                     >  2:993502/1‑150 (MQ=255)
                                                                                                               cGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAgcg                          >  1:1118251/1‑150 (MQ=255)
                                                                                                               cGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAgcg                          >  1:434859/1‑150 (MQ=255)
                                                                                                                ggCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGCAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAgcgc                         <  2:544248/150‑1 (MQ=255)
                                                                                                                   gcTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAGCgcggc                      <  1:387554/150‑1 (MQ=255)
                                                                                                                     tGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAGCGcggcgg                    >  1:154603/1‑150 (MQ=255)
                                                                                                                               gggTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAGCGCGGCGGCAAAGCCGag          <  2:200832/150‑1 (MQ=255)
                                                                                                                                ggTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTTGGCGAACAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAGCGCGGCGGCAAAGCCGagc         >  1:810621/1‑150 (MQ=255)
                                                                                                                                     cAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAGCGCGGCGGCAAAGCCGAGCAGCAc    >  1:1072317/1‑150 (MQ=255)
                                                                                                                                      aGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGAGCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAGCGCGGCGGCAAAGCCGAGCAGCACg   >  1:821901/1‑150 (MQ=255)
                                                                                                                                       ggCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAGCGCGGCGGCAAAGCCGAGCAGCACgg  >  1:359145/1‑150 (MQ=255)
                                                                                                                                           cACCCGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGAt                                                                                                             <  1:225025/39‑1 (MQ=255)
                                                                                                                                           cACCCGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGAt                                                                                                             >  2:784650/1‑39 (MQ=255)
                                                                                                                                           cACCCGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGAt                                                                                                             <  1:784650/39‑1 (MQ=255)
                                                                                                                                           cACCCGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGAt                                                                                                             >  2:225025/1‑39 (MQ=255)
                                                                                                                                               |                                                                                                                                             
GCATGCCCAGCGCCAATGCTGCTTCACGTTGCCCACGGGCCACGCCGAGGATGCCGGCGCGGAAGATCTCGCTCAGGTAGGCCCCGTAGTTCAGCGACAAGGCAATGATGCCGGCGCTGATGGCGCCGGGTACCAGGCCCACCTGCGGCAGCCCCAGGTAGATCAGCAGAATCTGGATCAGCAGCGGCGTGCCACGGAAGAACGAGGTGTAGAAGCTGGCCACACCCACCAGCACCGCGCTGTTGGATAACCGTGCCAGCGCGGCGGCAAAGCCGAGCAGCACGG  >  NC_002947/763983‑764267

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