Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,021,534 A→G 12.5% N677D (AAC→GAC)  PP_2638 → cellulose synthase operon protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,021,5340AG12.5% 75.2 / 4.0 32N677D (AAC→GAC) PP_2638cellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base G (2/2);  total (15/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 = 4.09e-01

CAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCCGATCCGTTGTTGCGGCGCGACGCAGCCCGCCTGGTAGCGCAGGATGATCCT  >  NC_002947/3021386‑3021683
                                                                                                                                                    |                                                                                                                                                     
cAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCaa                                                                                                                                                      >  1:662646/1‑150 (MQ=255)
                          tGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGgtgt                                                                                                                            >  1:446099/1‑150 (MQ=255)
                                        ccTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGt                                                                                                              <  2:604030/150‑1 (MQ=255)
                                           gATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGaa                                                                                                           <  2:446099/150‑1 (MQ=255)
                                            aTGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAac                                                                                                          <  1:471760/150‑1 (MQ=255)
                                                  gCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATc                                                                                                    >  1:24702/1‑150 (MQ=255)
                                                  gCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATc                                                                                                    >  1:661597/1‑150 (MQ=255)
                                                    cTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAg                                                                                                  >  1:124766/1‑150 (MQ=255)
                                                        gggTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGcggcgg                                                                                                                               >  2:377637/1‑117 (MQ=255)
                                                        gggTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGcggcgg                                                                                                                               <  1:377637/117‑1 (MQ=255)
                                                         ggTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAggcggc                                                                                             >  2:582407/1‑150 (MQ=255)
                                                           tGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCGACGCCCAGCGGCGCCTgg                                                                                                                                     >  2:664555/1‑108 (MQ=255)
                                                           tGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCGACGCCCAGCGGCGCCTgg                                                                                                                                     >  2:664570/1‑108 (MQ=255)
                                                           tGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCGACGCCCAGCGGCGCCTgg                                                                                                                                     <  1:664555/108‑1 (MQ=255)
                                                           tGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCGACGCCCAGCGGCGCCTgg                                                                                                                                     <  1:664570/108‑1 (MQ=255)
                                                             aTGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGATCGCCACTTCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGAGAAAACGATCAGGCGGCCAgg                                                                                         >  1:251053/1‑150 (MQ=255)
                                                                        cGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCg                                                                              <  2:662646/150‑1 (MQ=255)
                                                                           ggCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGccc                                                                                                                                                 >  2:794849/1‑80 (MQ=255)
                                                                           ggCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGccc                                                                                                                                                 <  1:794849/80‑1 (MQ=255)
                                                                             ccGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGa                                                                            >  1:207605/1‑147 (MQ=255)
                                                                             ccGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGa                                                                            <  2:207605/147‑1 (MQ=255)
                                                                                 aGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCcgc                                                                     <  1:85143/150‑1 (MQ=255)
                                                                                     aaCCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCgcag                                                                 <  2:756496/150‑1 (MQ=255)
                                                                                              cagGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGcc                                                                       <  2:280734/135‑1 (MQ=255)
                                                                                              cagGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGcc                                                                       >  1:280734/1‑135 (MQ=255)
                                                                                                ggCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGtgtg                                                                                                                           <  2:477942/81‑1 (MQ=255)
                                                                                                ggCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGtgtg                                                                                                                           >  1:477942/1‑81 (MQ=255)
                                                                                                           tGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCTACCCGATCCgttgt                                           <  2:251053/150‑1 (MQ=255)
                                                                                                           tGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCCGATCCgttgt                                           <  2:661597/150‑1 (MQ=255)
                                                                                                                 gCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCCGATCCGTTGTTgcggc                                     <  2:124766/150‑1 (MQ=255)
                                                                                                                           gCCACAATTTGCCGCTGATGACGGAAATGCTCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCCGATCCGTTGTTGCGGCGCGACGCAGc                           <  1:483/150‑1 (MQ=255)
                                                                                                                                      ccGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCCGATCCGTTGTTGCGGCGCGACGCAGCCCGCCTGGTAg                >  1:303503/1‑150 (MQ=255)
                                                                                                                                                 ggCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCCGATCCGTTGTTGCGGCGCGATGCAGCCCGCCTGGTAGCGCAGgatgat     >  1:294944/1‑150 (MQ=255)
                                                                                                                                                    aaCGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGAGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCCGATCCGTTGTTGCGGCGCGACGCAGCCCGCCTGGTAGCGCAGGATGATCCt  >  2:362937/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCTATTACCAGCGTGTACTGGTCGTGCAGCCTGGGCGCCCTGATGCACGCCTTGGGGTGATGGAAAGCGCCGTGGCCGAAGGCAACCTGCAGCAGGCCCGGCACATGCTGGGCGTGAAGGTGCCACAATTTGCCGCTGATGACGGCAACGCCCAGCGGCGCCTGGCGGCGGTGTGGACCGCGCTGGGTGAACACGATCAGGCGGCCAGGCTGCTTGACCGGATTGCCGCGCAGCAGACGCGACCCGATCCGTTGTTGCGGCGCGACGCAGCCCGCCTGGTAGCGCAGGATGATCCT  >  NC_002947/3021386‑3021683

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