Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,033,537 C→G 8.4% 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.4% 141.0 / 3.2 47A26A (GCC→GCG) PP_2646hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/23);  new base G (2/2);  total (22/25)
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.06e-01

GGCTACTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGG  >  NC_002947/3033391‑3033685
                                                                                                                                                  |                                                                                                                                                     
ggCTACTGCGTGCTCGCCCAGGCCATCAAGACCCTCAATATCCGCCCTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCt                                                                                                                                                     <  8:309747/149‑1 (MQ=255)
 gCTACTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTg                                                                                                                                                    <  4:276683/149‑1 (MQ=255)
 gCTACTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTg                                                                                                                                                    >  3:276683/1‑149 (MQ=255)
   tACTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCG‑CTGGCGCTGgac                                                                                                                                                 >  7:243792/1‑149 (MQ=255)
     cTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGgac                                                                                                                                                 <  1:23803/148‑1 (MQ=255)
        cGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGATCGATGCAAGGAGAGCATACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGgacgac                                                                                                                                              <  8:379044/148‑1 (MQ=255)
               gCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagc                                                                                                                                       <  5:413223/148‑1 (MQ=255)
                     gCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcgcc‑agcg                                                                                                                               <  7:343215/149‑1 (MQ=255)
                          gAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGa                                                                                                                          >  6:262148/1‑149 (MQ=255)
                          gAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGa                                                                                                                          >  5:23690/1‑149 (MQ=255)
                              aCCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCg                                                                                                                      <  7:79399/149‑1 (MQ=255)
                                  tctATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGgtg                                                                                                                  >  3:168625/1‑149 (MQ=255)
                                      tatCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGa                                                                                                              <  5:159400/149‑1 (MQ=255)
                                         ccgccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATcc                                                                                                           >  3:97497/1‑149 (MQ=255)
                                          cgccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCgg                                                                                                                              >  7:165545/1‑129 (MQ=255)
                                          cgccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCgg                                                                                                                              <  8:165545/129‑1 (MQ=255)
                                          cgccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGATCGGCCG‑CTGGCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGgtgtcggcgt                                                                                                           <  8:68146/148‑7 (MQ=17)
                                           gccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCAc                                                                                                         >  2:59424/1‑149 (MQ=255)
                                            ccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcg                                                                                                                                      >  8:221913/1‑120 (MQ=255)
                                            ccgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTccagcg                                                                                                                                      <  7:221913/120‑1 (MQ=255)
                                             cgTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACgg                                                                                                       <  3:46468/149‑1 (MQ=255)
                                                  cGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGc                                                                                                   <  8:130007/148‑1 (MQ=255)
                                                    gcgaGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGc                                                                                                >  4:213739/4‑149 (MQ=255)
                                                    gCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGc                                                                                                <  4:97497/149‑1 (MQ=255)
                                                     ccAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGc                                                                                                                             >  4:182250/1‑119 (MQ=255)
                                                     ccAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGc                                                                                                                             <  3:182250/119‑1 (MQ=255)
                                                     ccAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACgg                                                                                                       >  1:119072/1‑141 (MQ=255)
                                                     ccAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACgg                                                                                                       <  2:119072/141‑1 (MQ=255)
                                                          ggTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGttt                                                                                          >  7:135215/1‑149 (MQ=255)
                                                                     aTGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTgggg                                                                                >  7:231025/1‑148 (MQ=255)
                                                                        cAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAAt                                                                            <  1:59424/149‑1 (MQ=255)
                                                                         aaGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                           >  6:130902/1‑149 (MQ=255)
                                                                         aaGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATc                                                                           <  5:262148/149‑1 (MQ=255)
                                                                         aaGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑ACCGGCCGATCCGGGTGTTGATACACGGTGGCAGCCAGTTTCTATGTGGCGaaaac                                                                           >  6:86456/1‑147 (MQ=255)
                                                                           ggAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCg                                                                          >  5:138276/1‑148 (MQ=255)
                                                                                       aTGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGgcgcg                                                             >  6:94390/1‑149 (MQ=255)
                                                                                               aCTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCAt                                                      <  4:168625/148‑1 (MQ=255)
                                                                                                 tGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGgc                                                   >  8:114162/1‑149 (MQ=255)
                                                                                                           aGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTgg                                         >  2:285782/1‑149 (MQ=255)
                                                                                                              gcgcAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGc                                      <  5:94390/149‑1 (MQ=255)
                                                                                                              gcgcAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGGGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGc                                      <  7:114162/149‑1 (MQ=255)
                                                                                                                   aCTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTAc                                 <  3:213739/149‑1 (MQ=255)
                                                                                                                            acacCCGGATCGGCCG‑CTGGCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAAc                        >  4:122436/1‑149 (MQ=18)
                                                                                                                             cacCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCTGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGTTCTGGTGCGTTACCCGCCAAACt                       <  5:86456/149‑1 (MQ=255)
                                                                                                                             cacCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACt                       <  8:135215/149‑1 (MQ=255)
                                                                                                                                          ccgctggCGCTGGACGTCGTCCAGGGCCTGGCGGCAGCGCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCCGCACg           <  3:122436/144‑1 (MQ=25)
                                                                                                                                            ccAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGGGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGAAGCTGCACGGc         >  1:334008/1‑148 (MQ=255)
                                                                                                                                                 cccTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAg   <  5:130902/149‑1 (MQ=255)
                                                                                                                                                 cccTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTAGATCAACGTTGTCACCCAATAAATAAAAAAAAAAATCGAAAGCTGGGCGCGCGAGCATGGCGCAGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAg   >  4:268276/1‑149 (MQ=255)
                                                                                                                                                  ccTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAgg  >  1:167170/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCTACTGCGTGGTCGCCCAGGCCATGAAGACCCTCTATATCCGCCGTTTCGGCCAGTGGTTCTGAGCGATGCAAGGAGAGCAGACCATGCGAGTACTGATATGTGCAGGGCGCAACTACGCCGACACCCGGCGCTGCCGCCAGGCCCTGGACGACGTCCAGCGCC‑AGCGGCCGATCCGGGTGTTGATCCACGGTGGCAGCCAGTTTCTTGGTGGGGAAATCGAAAGCTGGGCGCGCGAGCATGGCGCCGATCTGGTGCGTTACCCGCCAAACTGGCAGCTGCACGGCAAGCAGG  >  NC_002947/3033391‑3033685

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