Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,033,537 C→G 8.6% A26A (GCC→GCG PP_2646 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,033,5370CG 8.6% 166.0 / 5.3 57A26A (GCC→GCGPP_2646hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base G (3/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

CTACTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCA  >  NC_002947/3033393‑3033679
                                                                                                                                                |                                                                                                                                                
cTACTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGAGCTGCAGCCAAGCCCTg                                                                                                                                               >  4:340964/1‑148 (MQ=255)
             gCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcg                                                                                                                                 >  2:61691/1‑149 (MQ=255)
              cccAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgc                                                                                                                                <  8:279799/149‑1 (MQ=255)
               ccAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgc                                                                                                                                >  8:84111/1‑148 (MQ=255)
                  ggCCATGAAGACCCTCTATATCCGGCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGTGAGCAGAACATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGTCGCTGCCGCCAGGCCCTGGACGACGTccagcgccagc                                                                                                                           <  8:214911/149‑1 (MQ=255)
                      aTGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCGCTGGACGTCGTccagggcctggcg                                                                                                                          >  4:2804/1‑142 (MQ=255)
                      aTGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCGCTGGACGTCGTccagggcctggcg                                                                                                                          <  3:2804/147‑6 (MQ=255)
                           gACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATcc                                                                                                                  >  3:57646/1‑149 (MQ=255)
                              ccTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCggg                                                                                                               >  6:130650/1‑149 (MQ=255)
                                tctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGgtg                                                                                                             >  6:103980/1‑149 (MQ=255)
                                tctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGATCGGCCGCTGGCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGgtg                                                                                                             >  2:249826/1‑149 (MQ=21)
                                 ctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGgtgt                                                                                                            >  3:185895/1‑149 (MQ=255)
                                  tataTCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGgtg                                                                                                             <  5:103980/147‑1 (MQ=255)
                                   atatCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGtt                                                                                                           >  8:275064/1‑148 (MQ=255)
                                    tatCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGgtgt                                                                                                            <  4:185895/146‑1 (MQ=255)
                                    tatCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGATCGGCCGCTGGCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGgtgtcg                                                                                                          <  1:249826/148‑3 (MQ=18)
                                     atCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGa                                                                                                         >  8:31470/1‑148 (MQ=255)
                                        cgccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCa                                                                                                     >  7:274985/1‑149 (MQ=255)
                                          ccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACg                                                                                                   >  7:120894/1‑149 (MQ=255)
                                          ccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACg                                                                                                   <  4:57646/149‑1 (MQ=255)
                                             tttCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGa                                                                                                         >  3:319723/1‑140 (MQ=255)
                                             tttCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGa                                                                                                         <  4:319723/140‑1 (MQ=255)
                                                    cAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCa                                                                                         <  3:340964/149‑1 (MQ=255)
                                                            cTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCtt                                                                                  <  1:61691/148‑1 (MQ=255)
                                                              gAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCa                                                                                             <  5:261533/135‑1 (MQ=255)
                                                              gAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCa                                                                                             >  6:261533/1‑135 (MQ=255)
                                                               aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                 <  7:31470/146‑1 (MQ=255)
                                                               aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                 <  7:275064/146‑1 (MQ=255)
                                                               aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                 <  8:120894/146‑1 (MQ=255)
                                                               aGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtg                                                                                 <  5:32423/146‑1 (MQ=255)
                                                                      cAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGaaa                                                                        <  5:130650/148‑1 (MQ=255)
                                                                       aaGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                      >  6:302082/1‑149 (MQ=255)
                                                                       aaGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                      >  7:106809/1‑149 (MQ=255)
                                                                            gagCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAg                                                                 >  6:344931/1‑149 (MQ=255)
                                                                               cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtggtgg                                                                             <  6:142638/134‑1 (MQ=255)
                                                                               cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTtggtgg                                                                             >  5:142638/1‑134 (MQ=255)
                                                                               cAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTg                                                              >  7:235679/1‑149 (MQ=255)
                                                                                aGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTgg                                                             <  5:302082/149‑1 (MQ=255)
                                                                                  aCCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTggg                                                            >  1:189082/1‑148 (MQ=255)
                                                                                  aCCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTggg                                                            <  2:189082/148‑1 (MQ=255)
                                                                                    cATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGgcgc                                                        >  8:278771/1‑149 (MQ=255)
                                                                                     aTGCGAGTACTGATATGTGCAGGGCGCAACTACGACGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGgcgcg                                                       >  7:300484/1‑149 (MQ=255)
                                                                                      tGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAg                                                                 >  1:93929/1‑139 (MQ=255)
                                                                                      tGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAg                                                                 <  2:93929/139‑1 (MQ=255)
                                                                                                 atatGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGgcgc                                           <  5:344931/149‑1 (MQ=255)
                                                                                                      tgCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATc                                      <  4:287302/149‑1 (MQ=255)
                                                                                                        cAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTg                                    <  7:278771/149‑1 (MQ=255)
                                                                                                        cAGGGCGCAACTACGACGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCt                                     <  8:300484/148‑1 (MQ=255)
                                                                                                           ggCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTg                                 <  8:274985/149‑1 (MQ=255)
                                                                                                           ggCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTg                                 >  2:270349/1‑149 (MQ=255)
                                                                                                                     cGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGc                       >  7:302210/1‑149 (MQ=255)
                                                                                                                          acacCCGGATCGGCCGCTGGCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAAc                  >  2:59089/1‑149 (MQ=18)
                                                                                                                              ccGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCAGGGGGGCAGCCAGTATCTTGGTgggg                                                                           >  4:268049/1‑89 (MQ=255)
                                                                                                                                ggCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATc                                      <  1:254796/123‑1 (MQ=255)
                                                                                                                                ggCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTAGATCAACGTTGTCACCCAATAAATTAAAAAAAAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATc                                      >  2:254796/1‑123 (MQ=255)
                                                                                                                                          ccAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGACGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGc   <  1:270349/149‑1 (MQ=255)
                                                                                                                                          ccAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCa  <  8:106809/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CTACTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCCAGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCA  >  NC_002947/3033393‑3033679

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