Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,043,722 T→A 11.4% intergenic (+18/+16) gmk → / ← PP_5297 deoxyguanylate kinase/guanylate kinase/amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,043,7220TA11.4% 124.4 / 7.8 44intergenic (+18/+16)gmk/PP_5297deoxyguanylate kinase/guanylate kinase/amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (23/16);  new base A (2/3);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

TACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACACCAGC  >  NC_002947/6043576‑6043870
                                                                                                                                                  |                                                                                                                                                    
tACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGc                                                                                                                                                   >  2:940945/1‑150 (MQ=255)
 aCGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCa                                                                                                                                                  <  2:212828/150‑1 (MQ=255)
        catcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCat                                                                                                                                           <  1:387927/150‑2 (MQ=255)
         atcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          <  1:763331/150‑1 (MQ=255)
         atcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCatc                                                                                                                                          <  1:267555/150‑3 (MQ=255)
         atcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGCTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          <  1:239582/150‑1 (MQ=255)
          tcaATGATGTTTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATGGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATg                                                                                                                                         <  2:98015/150‑1 (MQ=255)
             atgatgATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      <  1:458944/150‑1 (MQ=255)
              tgatgaTTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAATTAATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGttt                                                                                                                                     >  1:674366/1‑150 (MQ=255)
                  gatTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTc                                                                                                                                 <  2:363279/150‑1 (MQ=255)
                   atTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTcc                                                                                                                                >  2:49271/1‑150 (MQ=255)
                     ttCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTg                                                                                                                              <  1:84195/150‑1 (MQ=255)
                        gACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAg                                                                                                                           >  2:164401/1‑150 (MQ=255)
                                     cAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGcc                                                                                                              >  2:476111/1‑150 (MQ=255)
                                          ttGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCctgctg                                                                                                         >  1:706455/1‑150 (MQ=255)
                                                 gaGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                                  <  1:2561/149‑1 (MQ=255)
                                                    tgtTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTa                                                                                               >  1:556856/1‑150 (MQ=255)
                                                          gCTCGAATCGGCCGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCg                                                                                         >  2:197389/1‑150 (MQ=255)
                                                               aaTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGaa                                                                                    <  1:164401/150‑1 (MQ=255)
                                                                            tgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGcgca                                                                       >  2:871932/1‑150 (MQ=255)
                                                                                gaagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGcgcacgca                                                                   >  1:443991/1‑150 (MQ=255)
                                                                                 aagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCCATTTAGGTTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCAt                                                                  <  1:476111/150‑1 (MQ=255)
                                                                                       cagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGa                                                            >  1:1014976/1‑150 (MQ=255)
                                                                                          cagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAg                                                         >  1:169051/1‑150 (MQ=255)
                                                                                                        aTTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGa                                           <  1:940945/150‑1 (MQ=255)
                                                                                                         ttGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGCGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGAt                                          <  2:771114/150‑1 (MQ=255)
                                                                                                             tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCa                                      >  2:667576/1‑150 (MQ=255)
                                                                                                               aaGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGa                                                                                        <  1:827909/98‑1 (MQ=255)
                                                                                                               aaGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGa                                                                                        >  2:827909/1‑98 (MQ=255)
                                                                                                                 gCAAATGATCGGCTGAAGATCGAGGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCa                                  >  1:830937/1‑150 (MQ=255)
                                                                                                                     aTGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCCGGCCGATGCCAAGCACCga                              >  2:478495/1‑150 (MQ=255)
                                                                                                                           ggCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAgatgat                        >  1:113085/1‑150 (MQ=255)
                                                                                                                              tGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATcgc                     <  1:49271/150‑1 (MQ=255)
                                                                                                                                  gATCGATGGGGCTGCATTGCAGCCCCCATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCa                 >  2:215807/1‑150 (MQ=255)
                                                                                                                                  gATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCa                 >  2:967510/1‑150 (MQ=255)
                                                                                                                                      attGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCaa                                     >  2:487051/3‑126 (MQ=255)
                                                                                                                                      attGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCaa                                     <  1:487051/124‑1 (MQ=255)
                                                                                                                                       ttGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCAc            >  2:117845/2‑150 (MQ=255)
                                                                                                                                       aTGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATAAc            >  2:737746/1‑150 (MQ=255)
                                                                                                                                         ggggCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATg                                         <  1:964363/119‑1 (MQ=255)
                                                                                                                                         ggggCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATg                                         >  2:964363/1‑119 (MQ=255)
                                                                                                                                            gCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAaca       <  2:443991/150‑1 (MQ=255)
                                                                                                                                              tGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCCGCAATCACGAACAcc     >  2:134950/1‑150 (MQ=255)
                                                                                                                                                 aTTGCAGCCCCAATGGTTTAGTCCTGCAGCGGAAAGGTGGATTGGTGCAGTTGGGTAGCAACCATTAGGAACCCCCGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGACGATCGCGGCAATCACGAACAccagc  >  1:300257/1‑150 (MQ=255)
                                                                                                                                                  ttGCAGCCCCAATGGTTTGGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                       <  2:204112/64‑1 (MQ=255)
                                                                                                                                                  ttGCAGCCCCAATGGTTTGGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                       >  1:204112/1‑64 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACACCAGC  >  NC_002947/6043576‑6043870

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