Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,678,948 C→T 9.2% A57A (GCG→GCA)  PP_3243 ← GNAT family acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,678,9480CT 9.2% 153.1 / 4.5 54A57A (GCG→GCA) PP_3243GNAT family acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/24);  new base T (3/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

TCGACCACTTGCTGCAGCAACCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCTCGACGTCCGGCAGCCACGTGATGTCT  >  NC_002947/3678811‑3679092
                                                                                                                                         |                                                                                                                                                
tCGACCACTTGCTGCAGCAACCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGTGCCACCCAAc                                                                                                                                        >  1:86280/1‑148 (MQ=255)
    ccACTTGCTGCAGCAACCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTAccc                                                                                                                                   >  7:131774/1‑149 (MQ=255)
        ttGCTGCAGCAACCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCAccc                                                                                                                                           <  1:193910/137‑1 (MQ=255)
        ttGCTGCAGCAACCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCAccc                                                                                                                                           <  1:337267/137‑1 (MQ=255)
        ttGCTGCAGCAACCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCAccc                                                                                                                                           >  2:193910/1‑137 (MQ=255)
        ttGCTGCAGCAACCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGAGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCAccc                                                                                                                                           >  2:337267/1‑137 (MQ=255)
           ctgcAGCAACCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGcc                                                                                                                             <  2:263033/148‑1 (MQ=255)
                  aaCCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTcgacga                                                                                                                      >  8:29186/1‑148 (MQ=255)
                   aCCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGTGCCACCCAACTACCCCCTGCCTCGACGaaa                                                                                                                    <  2:86280/149‑1 (MQ=255)
                      ggCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGcc                                                                                                                 <  4:109667/149‑1 (MQ=255)
                      ggCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGcc                                                                                                                 <  7:29186/149‑1 (MQ=255)
                           ccGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCtgtg                                                                                                            <  1:143827/149‑1 (MQ=255)
                                  cTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCcgcg                                                                                                     >  6:325439/1‑149 (MQ=255)
                                  cTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCcgcg                                                                                                     >  3:264928/1‑149 (MQ=255)
                                   ttGACCTTGGGCCTCCTGGCGCACAGATATTTCATGAAGGTGCAATGCCTGATCAGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACAACCCCCAGACTCGACGAAAGCCCTGTGCGCcgcg                                                                                                     >  8:161931/1‑148 (MQ=255)
                                    tGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCgcgc                                                                                                    >  5:322044/1‑148 (MQ=255)
                                     gCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCaa                                                                                                  >  7:281842/1‑149 (MQ=255)
                                            ggCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGa                                                                                           <  8:15697/149‑1 (MQ=255)
                                                  ctGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCg                                                                                      <  5:325439/148‑1 (MQ=255)
                                                    ggCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCt                                                                                   >  7:336552/1‑149 (MQ=255)
                                                        cACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTgcg                                                                               >  4:122563/1‑149 (MQ=255)
                                                         aCCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGTCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTACTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTgcgc                                                                              <  7:161931/149‑1 (MQ=255)
                                                           cGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCgcca                                                                            <  8:281842/149‑1 (MQ=255)
                                                              tatTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTgcg                                                                                                                                                 <  8:177499/77‑1 (MQ=255)
                                                              tatTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTgcg                                                                                                                                                 >  7:177499/1‑77 (MQ=255)
                                                                   cATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCgcca                                                                            >  2:157643/1‑141 (MQ=255)
                                                                   cATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCgcca                                                                            <  1:157643/141‑1 (MQ=255)
                                                                   cATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCACAGTCGCGTCCTGTGCCACCCAACTACCCACTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAgccaggc                                                                     >  6:121990/1‑148 (MQ=255)
                                                                          gTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATcc                                                                                               <  3:227819/115‑1 (MQ=255)
                                                                          gTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATcc                                                                                               >  4:227819/1‑115 (MQ=255)
                                                                          gTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCg                                                             <  8:336552/149‑1 (MQ=255)
                                                                            gCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCg                                                             >  3:16643/1‑147 (MQ=255)
                                                                            gCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCg                                                             <  4:16643/147‑1 (MQ=255)
                                                                             cAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGc                                                           <  6:322044/148‑1 (MQ=255)
                                                                             cAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCAACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACAAACCCCTGCCTCGACGAAAGCCCTGACCGCCGCGCAATCCATGACCTCGCCTTGCACCAGCCAGGCCAGTGCCGGc                                                           >  1:109732/1‑148 (MQ=255)
                                                                                  ccTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTgc                                                                                <  4:85707/122‑1 (MQ=255)
                                                                                  ccTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTgc                                                                                >  3:85707/1‑122 (MQ=255)
                                                                                            cAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCt                                           <  8:131774/149‑1 (MQ=255)
                                                                                                       aGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACaa                                                 >  3:288113/1‑132 (MQ=255)
                                                                                                       aGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACaa                                                 <  4:288113/132‑1 (MQ=255)
                                                                                                       aGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAg                                >  7:209098/1‑149 (MQ=255)
                                                                                                        ggCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAgc                               >  7:217413/1‑149 (MQ=255)
                                                                                                           aaaaCCCACAATGCGCCCAGTCGCGTCCTGTGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCt                            >  5:305838/1‑149 (MQ=255)
                                                                                                           aaaaCCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCt                            >  1:325629/1‑149 (MQ=255)
                                                                                                            aaaCCCACAATGCGCCCAGTCGCGTCCTGTGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCTc                           <  6:305838/149‑1 (MQ=255)
                                                                                                              aCCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGTCGATAGCCCTGTTCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAATAGCCTGGGCTGCCGAGCGCTCGa                         <  2:109732/149‑1 (MQ=255)
                                                                                                                ccACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCTCGACg                       >  5:81962/1‑149 (MQ=255)
                                                                                                                   caATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCTCGACGTcc                    <  8:217413/149‑1 (MQ=255)
                                                                                                                         gcCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCTCGACGTCCg                   <  4:207310/144‑1 (MQ=255)
                                                                                                                         gcCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCTCGACGTCCg                   >  3:207310/1‑144 (MQ=255)
                                                                                                                          cccAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTg                                          >  6:221176/1‑120 (MQ=255)
                                                                                                                          cccAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTg                                          <  5:221176/120‑1 (MQ=255)
                                                                                                                           ccAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAgccaggc                                                                     <  4:254884/92‑1 (MQ=255)
                                                                                                                           ccAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAgccaggc                                                                     >  3:254884/1‑92 (MQ=255)
                                                                                                                           ccAGTCGCGTCCTGCGCCACCCAACTA‑CCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCTCGACGTCCGGCAGCCAc           >  3:323323/1‑149 (MQ=255)
                                                                                                                                      cTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCTCGACGTCCGGCAGCCACGTGATGTCt  <  8:209098/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                
TCGACCACTTGCTGCAGCAACCGGCGCCCGACACCTTGCCCTTGGGCCTCCTGGCGCACCGATATTTCATGAAGGTGCAATGCCTGATCCGCCAGCCGGGCGCAGGCAAAACCCACAATGCGCCCAGTCGCGTCCTGCGCCACCCAACTACCCCCTGCCTCGACGAAAGCCCTGTGCGCCGCGCAATCCATGACGTCGCCTTGCGCCAGCCAGGCCAGTGCCGGCCAACTGACAAAAGCCTGGGCTGCCGAGCGCTCGACGTCCGGCAGCCACGTGATGTCT  >  NC_002947/3678811‑3679092

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: ‑