Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,805,005 Δ1 bp 8.3% coding (154/618 nt) PP_3363 ← thymidylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,805,0050C. 8.3% 147.9 / 7.3 48coding (154/618 nt)PP_3363thymidylate kinase
Reads supporting (aligned to +/- strand):  ref base C (25/19);  new base . (2/2);  total (27/21)
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.32e-01

CGGTCAATCAGGATGATGCTGTCAGCGGGCTGCCTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGG  >  NC_002947/3804866‑3805153
                                                                                                                                           |                                                                                                                                                     
cGGTCAATCAGGATGATGCTGTCAGCGGGCTGCCTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACCCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGtctc                                                                                                                                              <  2:265058/149‑1 (MQ=255)
       tCAGGATGATGCTGTCAGCGGGCTGCCTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGGTTCTGGAAAGCTGGTTGACCAGGG‑CATGGTCTCGCCCCGGt                                                                                                                                      >  5:321300/1‑149 (MQ=255)
       tCAGGATGATGCTGTCAGCGGGCTGCCTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCgg                                                                                                                                       >  7:294271/1‑149 (MQ=255)
       tCAGGATGATGCTGTCAGCGGGCTGCCTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCgg                                                                                                                                       <  2:26792/149‑1 (MQ=255)
       tCAGGATGATGCTGTCAGCGGGCTGCCTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGG‑CATGGTCTCGCCCCGGt                                                                                                                                      >  5:327293/1‑149 (MQ=255)
          ggATGATGCTGTCAGCGGGCTGCCTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAGGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGt                                                                                                                                      >  5:239405/1‑147 (MQ=255)
                tgCTTTCAGCGGGCTGCCTGGGCAGTACATGTCGGTAAATCAAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATc                                                                                                                               <  8:285712/148‑1 (MQ=255)
                         cGGGCTGCCTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAGGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATtctttctt                                                                                                                     <  6:239405/149‑1 (MQ=255)
                         cGGGCTGCCTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGATAGCTGGTTGACGAGGGCCATGGTCTCGCCGCGGTGAGGATCAttctttct                                                                                                                      <  2:206251/148‑1 (MQ=255)
                                ccTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCa                                                                                                                >  3:2676/1‑147 (MQ=255)
                                      cAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGa                                                                                                        >  5:306019/1‑149 (MQ=255)
                                            aTGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCt                                                                                                  >  2:222431/1‑149 (MQ=255)
                                            aTGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCt                                                                                                  >  1:93247/1‑149 (MQ=255)
                                               tCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTAGCAAAGGCGGATCACCTTGt                                                                                               >  5:285213/1‑149 (MQ=255)
                                                 gggAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCTGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTc                                                                                             >  5:327561/1‑149 (MQ=255)
                                                  ggAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTcc                                                                                            >  3:59137/1‑149 (MQ=255)
                                                     aaTCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCGCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGt                                                                                         >  6:24828/1‑149 (MQ=255)
                                                     aaTCCAGGCTCGGCTTTAAGCGAAGCGTTCGCAAACCTAGAACTCAACGTCCTGGCAAGGGCTGCTGGAAAGCTGGGTGCCCAGGGCCATGGTCTCTCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGCCGGATCACCTTGTTCCCg                                                                                          <  8:331306/148‑1 (MQ=255)
                                                         cAGGCTCGGCTTTCAGCAAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGATGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAg                                                                                     >  7:101795/1‑149 (MQ=255)
                                                                ggCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAAc                                                                              >  4:106198/1‑149 (MQ=255)
                                                                   tttCAGCGAAGCTTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTccc                                                                                           <  3:290759/133‑1 (MQ=255)
                                                                   tttCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTccc                                                                                           >  4:290759/1‑133 (MQ=255)
                                                                   tttCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCACGTCCAGCCGTAACGCt                                                                           >  1:63055/1‑149 (MQ=255)
                                                                      cAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGG‑CATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGa                                                                        <  6:327293/148‑1 (MQ=255)
                                                                      cAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGG‑CATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGa                                                                        <  6:321300/148‑1 (MQ=255)
                                                                          gAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGa                                                                    >  1:74767/1‑149 (MQ=255)
                                                                          gAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGa                                                                    >  1:80055/1‑149 (MQ=255)
                                                                           aaGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGAc                                                                   <  4:2676/149‑1 (MQ=255)
                                                                                   caggaCCTTGAGTTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGtt                                                                                              >  8:52841/5‑114 (MQ=255)
                                                                                   caggaCCTTGAGTTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGtt                                                                                              <  7:52841/110‑1 (MQ=255)
                                                                                        ccTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCagcagc                                                      <  2:93247/149‑1 (MQ=255)
                                                                                                  aaGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGc                                           >  5:165628/1‑149 (MQ=255)
                                                                                                  aaGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGc                                           >  3:350798/1‑149 (MQ=255)
                                                                                                  aaGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGc                                           <  3:106198/149‑1 (MQ=255)
                                                                                                     gTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCACCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACATTGTTCACGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCgg                                        >  8:54751/1‑149 (MQ=255)
                                                                                                         tGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTAcc                                    >  7:74477/1‑149 (MQ=255)
                                                                                                                ggCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGc                             >  6:329557/1‑149 (MQ=255)
                                                                                                                        gAAAGCTGTTTGACCAGGGCGATGGTCTCGCCCCGGTGAGTATCATTCTTTCTTTCGCAATGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCGAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGaa                     <  4:342064/149‑1 (MQ=255)
                                                                                                                        gAAAGCTGGTTGACCAGTGCCATGGTGTCGCCCCGGGGAGGATCATTCTTTCTTTCGCAAAGGCTGATCACCTTGTTCCCGTCCAGCCGTTACGCTTGAGTGACGGCCTCCATCAGC‑GTGGTCTTGCCGGTACCCTTGTGCCCATc                        <  6:285213/146‑1 (MQ=255)
                                                                                                                         aaaGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGAAg                    <  5:217756/149‑1 (MQ=255)
                                                                                                                          aaGCTGGTTGACCAGGTCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGAAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGAAg                    <  2:63055/148‑1 (MQ=255)
                                                                                                                                gTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATaa             <  6:177779/149‑1 (MQ=255)
                                                                                                                                 ttGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGCCGGCCTCCAGCAGCAAAGGTCTTGCCGCTACCCTTGGGCCCATCGAAGGAAATa              >  1:311867/1‑148 (MQ=255)
                                                                                                                                     ccAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGc                                           >  7:42643/1‑114 (MQ=255)
                                                                                                                                     ccAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGc                                           <  8:42643/114‑1 (MQ=255)
                                                                                                                                      cAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAgcg       <  2:74767/149‑1 (MQ=255)
                                                                                                                                      cAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAgcg       <  2:80055/149‑1 (MQ=255)
                                                                                                                                        gggCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAgcggg     >  6:245322/1‑149 (MQ=255)
                                                                                                                                           ccATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAgcgggcgg  <  6:165628/149‑1 (MQ=255)
                                                                                                                                           ccATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAgcgggcgg  <  5:329557/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
CGGTCAATCAGGATGATGCTGTCAGCGGGCTGCCTGGGCAGTACATGTCGGGAAATCCAGGCTCGGCTTTCAGCGAAGCGTTCGCAAACCTTGAACTCAAGGTCCTGGCAAGGGCTTCTGGAAAGCTGGTTGACCAGGGCCATGGTCTCGCCCCGGTGAGGATCATTCTTTCTTTCGCAAAGGCGGATCACCTTGTTCCCGTCCAGCCGTAACGCTTGAGTGACGGCCTCCAGCAGC‑GTGGTCTTGCCGGTACCCTTGGGCCCATCGAAGGAAATAAACAGCGGGCGG  >  NC_002947/3804866‑3805153

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