Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,508,790 T→C 12.2% M109V (ATG→GTG) ‡ PP_2203 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,508,7900TC12.2% 96.5 / 5.8 41M109V (ATG→GTG) ‡PP_2203hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/14);  new base C (3/2);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTT  >  NC_002947/2508643‑2508934
                                                                                                                                                   |                                                                                                                                                
gcggcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATg                                                                                                                                                 <  3:154723/149‑1 (MQ=255)
    cACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCg                                                                                                                                             <  3:242631/149‑1 (MQ=255)
       cAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCa                                                                                                                                          >  3:78792/1‑149 (MQ=255)
       cAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGCCCGGCa                                                                                                                                          <  6:59154/149‑1 (MQ=255)
         aaGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGCCCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCAcc                                                                                                                                        >  4:179839/1‑149 (MQ=255)
                    aaaCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCtt                                                                                                                             >  3:246191/1‑149 (MQ=255)
                    aaaCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAgg                                                                                                                                  <  1:194095/143‑1 (MQ=255)
                    aaaCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAgg                                                                                                                                  >  2:194095/1‑143 (MQ=255)
                       cAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAg                                                                                                                          <  7:6001/149‑1 (MQ=255)
                        aGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAgg                                                                                                                         >  2:271601/1‑149 (MQ=255)
                                  ggCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGtt                                                                                                               <  4:78792/149‑1 (MQ=255)
                                   gCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCGGGCACCGCCAGGCCCTTGAGGTCGTTGCGttt                                                                                                              <  2:20180/149‑1 (MQ=255)
                                    cTGTGCTTCTGTAGTGGGCGACCCCGGGGTGGTGGCCATCGTGGGTCACCCCGATGACGTGGTAGTCAGCCTTGTGGTCACCCATCTCCTTGGCGGGCGACCCCTTGGGCTTGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTg                                                                                                             <  2:101904/149‑1 (MQ=255)
                                             agtagtTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAgg                                                                                                    >  5:56206/1‑149 (MQ=255)
                                                              gtgtCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCAGGGGCATGCCGGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTcggccgg                                                                                   <  4:246191/149‑1 (MQ=255)
                                                                     ggCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCTAGCAGGCGCACCTGCTCGGCCGGCACGTgg                                                                            >  4:236501/1‑149 (MQ=255)
                                                                                    ccccGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAAGAGCACCTGCTCGGCCGGCACCTGGCCTTCGACGAACTTg                                                             >  1:256980/1‑149 (MQ=255)
                                                                                     cccGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGCACTTgc                                                            >  2:287806/1‑149 (MQ=255)
                                                                                       cGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgccg                                                          >  3:34843/1‑149 (MQ=255)
                                                                                                 gTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAAc                                                                >  2:67014/1‑133 (MQ=255)
                                                                                                 gTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAAc                                                                <  1:67014/133‑1 (MQ=255)
                                                                                                     gCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGgt                                            >  4:235471/1‑149 (MQ=255)
                                                                                                      cATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGgtg                                           <  1:271601/149‑1 (MQ=255)
                                                                                                                    cccATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAg                             <  4:34843/149‑1 (MQ=255)
                                                                                                                           ccATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgc                      >  3:165219/1‑149 (MQ=255)
                                                                                                                               gCCGGGCATGCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCAACTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCa                  >  8:94315/1‑149 (MQ=255)
                                                                                                                               gCCGGGCATGCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGAAGGCGCAACTGCTCGGCCGGAACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCa                  >  8:89791/1‑149 (MQ=255)
                                                                                                                                ccGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCAc                                                                                               <  5:11426/71‑1 (MQ=255)
                                                                                                                                ccGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCAc                                                                                               >  6:11426/1‑71 (MQ=255)
                                                                                                                                ccGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCTCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCAc                 >  3:60349/1‑149 (MQ=255)
                                                                                                                                ccGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGCCGCAAGGCGTGGCGCCAc                 >  2:247361/1‑149 (MQ=255)
                                                                                                                                  gggcatgCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCaa                              <  5:34953/128‑1 (MQ=255)
                                                                                                                                  gggcatgCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCaa                              >  6:34953/7‑134 (MQ=255)
                                                                                                                                  gggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACg                >  8:246996/1‑148 (MQ=255)
                                                                                                                                  gggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGACCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGc               >  3:204612/1‑149 (MQ=255)
                                                                                                                                          ccATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGCCGAACTTGCCGCCGATCACCCCGGTGTGGCCCGACGCCAGGCGTGGCGCCCCGCCGAGGCGt       >  4:166515/1‑149 (MQ=255)
                                                                                                                                          ccATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGt       >  5:85984/1‑149 (MQ=255)
                                                                                                                                            aTGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCacgac                                  <  4:250066/120‑1 (MQ=255)
                                                                                                                                            aTGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCacgac                                  >  3:250066/1‑120 (MQ=255)
                                                                                                                                               ggCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCtt  <  3:235471/149‑1 (MQ=255)
                                                                                                                                               ggCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCtt  >  2:168936/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGAGGCGTTGCTT  >  NC_002947/2508643‑2508934

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