Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 718,857 T→C 9.6% T354A (ACC→GCC)  PP_0613 ← amidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947718,8570TC 9.6% 109.0 / 2.7 42T354A (ACC→GCC) PP_0613amidase family protein
Reads supporting (aligned to +/- strand):  ref base T (21/16);  new base C (2/2);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTGAAGGTTTCCAGCGGGTC  >  NC_002947/718720‑719001
                                                                                                                                         |                                                                                                                                                
cATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGcc                                                                                                                                         >  1:237138/1‑147 (MQ=255)
 aTCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAg                                                                                                                                      <  6:386577/149‑1 (MQ=255)
    cccGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTTGGCGAATACACCCAGGc                                                                                                                                    >  8:292292/1‑148 (MQ=255)
    cccGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTTGGCGAATACACCCAGGc                                                                                                                                    <  7:292292/148‑1 (MQ=255)
     ccGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTTGGCGAATACACCCAGGcgc                                                                                                                                  >  7:246721/1‑147 (MQ=255)
        aCCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTTGGCGAATACACCCAGGcgcttg                                                                                                                               <  8:246721/149‑6 (MQ=255)
         cccAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCt                                                                                                                              <  7:313367/149‑1 (MQ=255)
               cggcACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCAc                                                                                                                        <  2:237138/149‑1 (MQ=255)
                          gtCCGGCTTCGAACGCCCCCAACGGCAGCATCGGCAATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGc                                                                                                              >  1:255267/1‑148 (MQ=255)
                          gtCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGc                                                                                                              >  7:316699/1‑148 (MQ=255)
                           tCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCATGTTCGGCa                                                                                                             >  3:37390/1‑148 (MQ=255)
                             cGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATc                                                                                                          >  5:261136/1‑149 (MQ=255)
                               gCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCa                                                                                                        <  5:21163/149‑1 (MQ=255)
                                      aCGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGc                                                                                                 <  7:6759/149‑1 (MQ=255)
                                          caccaACGGCAGCATCGGCGATACCCGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCg                                                                                             >  3:331246/1‑149 (MQ=255)
                                             caACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGt                                                                                          <  8:33903/149‑1 (MQ=255)
                                                cGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCaa                                                                                       <  2:66949/149‑1 (MQ=255)
                                                cGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCaa                                                                                       <  8:316699/149‑1 (MQ=255)
                                                cGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAAATCTCCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCaa                                                                                       <  7:65416/149‑1 (MQ=255)
                                                        tCGGCGCTACCAGCACGTCATATCGCTGATGAAAGGCATTCATGGTCGCAATCAATTCTGCGCGTGCCTGAAGAGCCTGGGGGTAGTCGCCAAGGCTGATCTGTGCAGGTGGTTCGGACATCCAGCGCAGCGCCGGGTCCAggacctct                                                                               >  6:109841/1‑142 (MQ=255)
                                                        tCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCt                                                                               <  6:261136/149‑1 (MQ=255)
                                                             gATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTg                                                                          <  2:171911/149‑1 (MQ=255)
                                                                 ccAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCa                                                                      >  3:113452/1‑149 (MQ=255)
                                                                    gCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCgg                                                                                            >  6:263661/1‑124 (MQ=255)
                                                                            tatCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAgg                                                            >  1:145684/1‑148 (MQ=255)
                                                                                                  tCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGcgccgcgccg                                     >  8:79769/1‑149 (MQ=255)
                                                                                                      cGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCa                                  >  5:3989/1‑148 (MQ=255)
                                                                                                         aatcaatTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCa                                                                                    <  3:122369/95‑1 (MQ=255)
                                                                                                         aatcaatTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCa                                                                                    >  4:122369/1‑95 (MQ=255)
                                                                                                         aatcaatTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCCCCGCGCCGGCAAAc                               >  3:28170/1‑148 (MQ=255)
                                                                                                          atcaatTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACca                             >  1:8957/1‑149 (MQ=255)
                                                                                                           tcaatTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACcac                            >  6:35690/1‑149 (MQ=255)
                                                                                                              atTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAgg                         <  6:3989/149‑1 (MQ=255)
                                                                                                               ttCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAggg                        >  5:140610/1‑149 (MQ=255)
                                                                                                                 cTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGgtg                      >  7:218600/1‑149 (MQ=255)
                                                                                                                  tGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGgtgt                     >  3:170153/1‑149 (MQ=255)
                                                                                                                   gCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGtt                    >  3:141570/1‑149 (MQ=255)
                                                                                                                    cccGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGgtg                      <  2:138621/146‑1 (MQ=255)
                                                                                                                      cGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTg                   >  1:95059/1‑147 (MQ=255)
                                                                                                                      cGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTg                   >  3:251241/1‑147 (MQ=255)
                                                                                                                              aaGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTGAAGGTTTcc          <  4:28170/148‑1 (MQ=255)
                                                                                                                                       ggTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTGAAGGTTTCCAGCGGGTc  <  2:145684/147‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                
CATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTGAAGGTTTCCAGCGGGTC  >  NC_002947/718720‑719001

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