Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,472,879 A→C 10.2% Y226D (TAT→GAT) ‡ nadD ← nicotinate‑nucleotide adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,472,8790AC10.2% 130.8 / 8.0 49Y226D (TAT→GAT) ‡nadDnicotinate‑nucleotide adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (26/18);  new base C (3/2);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

CCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGA  >  NC_002947/5472734‑5473026
                                                                                                                                                  |                                                                                                                                                   
ccTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAtata                                                                                                                                                   <  1:60465/149‑1 (MQ=255)
 cTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGt                                                                                                                                                 >  5:41374/1‑149 (MQ=255)
    gATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGTGCGCACGATATAGTTCg                                                                                                                                              >  7:300976/1‑149 (MQ=255)
    gATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGtt                                                                                                                                                <  7:291529/147‑1 (MQ=255)
      tGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTg                                                                                                                                            <  1:38030/149‑1 (MQ=255)
        tCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTTCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTgc                                                                                                                                           <  7:168608/148‑1 (MQ=255)
            gggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAATCTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTc                                                                                                                                      >  7:22724/1‑149 (MQ=255)
             ggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCt                                                                                                                                       <  5:17639/147‑1 (MQ=255)
              gCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAACTATATCGTGCGcccaac                                                                                                                                    <  3:255151/149‑5 (MQ=255)
               ccTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGat                                                                                                                                   <  2:283857/149‑1 (MQ=255)
               ccTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGat                                                                                                                                   <  4:10852/149‑1 (MQ=255)
                cTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGata                                                                                                                                  >  5:203168/1‑149 (MQ=255)
                 ttGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTTGTGCGCCTCGata                                                                                                                                  >  6:341772/1‑148 (MQ=255)
                   gACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAg                                                                                                                               >  6:177821/1‑149 (MQ=255)
                     cGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGc                                                                                                                             <  6:203168/149‑1 (MQ=255)
                      gTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGc                                                                                                                             >  1:43239/1‑148 (MQ=255)
                            cAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACg                                                                                                                      >  6:283813/1‑149 (MQ=255)
                                    ccTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCgg                                                                                                              >  4:134596/1‑149 (MQ=255)
                                                cGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAAcc                                                                                                   <  6:41374/148‑1 (MQ=255)
                                                cGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCt                                                                                                  >  5:68149/1‑149 (MQ=255)
                                                 gACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTc                                                                                                 <  6:68149/149‑1 (MQ=255)
                                                    cAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCAcc                                                                                              <  2:43239/149‑1 (MQ=255)
                                                      gTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGa                                                                                            >  1:59263/1‑149 (MQ=255)
                                                          ttcGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTg                                                                                        >  3:191430/1‑149 (MQ=255)
                                                          ttcGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAACTATATCGTgcgc                                                                                                                                         <  1:156006/100‑1 (MQ=255)
                                                          ttcGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGAACTATATCGTgcgc                                                                                                                                         >  2:156006/1‑100 (MQ=255)
                                                            cGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGc                                                                                       <  4:191430/148‑1 (MQ=255)
                                                             gCCGTAAATTTTCTGCTTGGTCATATAAAAATCGATCAACTCGTTGCATG‑AGGCGCCAGAGGCACCGTCAGATGGGCGCACGATATACTTCGTGCGCATCGATATAGGCCAGTACGGAGTCCGGCACCAGGAACCTCAACGAttttcag                                                                                     >  7:265559/1‑145 (MQ=255)
                                                               cGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCt                                                                                   <  5:283813/149‑1 (MQ=255)
                                                                 tAAATTTTTTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGTTG‑AGGAGCCAGAGGCACAGTCATTTGGGCGGACGATATAGTTAGTGCGCCTCGATATATGCGAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTg                                                                                  <  8:265559/148‑1 (MQ=255)
                                                                   aaTTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGc                                                                                <  5:177821/148‑1 (MQ=255)
                                                                     ttttCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCcag                                                                             >  5:169389/1‑149 (MQ=255)
                                                                                 tCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGa                                                                  >  6:148828/1‑148 (MQ=255)
                                                                                   atatAAATCTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCt                                                               <  8:22724/149‑1 (MQ=255)
                                                                                     ataAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCtgt                                                             >  1:90388/1‑149 (MQ=255)
                                                                                         atcTCGTTCAACTAGTTGGATT‑AGTCGCCAGAGGCACCGGCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAg                                                         <  8:300976/147‑1 (MQ=255)
                                                                                              gTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGa                                                     >  5:659/1‑148 (MQ=255)
                                                                                                  aaCTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCg                                                >  7:154160/1‑149 (MQ=255)
                                                                                                        ttGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAgcg                                          >  2:65499/1‑149 (MQ=255)
                                                                                                                 ggCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGc                                  >  5:5469/1‑149 (MQ=255)
                                                                                                                     ccAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGa                              >  6:257750/1‑149 (MQ=255)
                                                                                                                     ccAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGa                              <  2:59263/149‑1 (MQ=255)
                                                                                                                              ccGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGacgaac                         >  2:47666/1‑145 (MQ=255)
                                                                                                                              ccGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGacgaac                         <  1:47666/145‑1 (MQ=255)
                                                                                                                                gTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAata                   >  6:321051/1‑149 (MQ=255)
                                                                                                                                   aGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATAtt                 >  4:160917/1‑148 (MQ=255)
                                                                                                                                    gTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTc               >  8:237563/1‑149 (MQ=255)
                                                                                                                                       aggCGCACGAACTATATCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccg            >  4:19108/2‑149 (MQ=255)
                                                                                                                                          cgcACGAACTATATCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCccgccg         >  5:268598/1‑149 (MQ=255)
                                                                                                                                           gcACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAAAGAAATATTTCCCgccgg        >  5:99815/1‑149 (MQ=255)
                                                                                                                                                 tataGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGa  <  5:148828/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATG‑AGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCGCTGGCCAGCAGCTGTCGGATCTGTGTAGCCGACACCGCAAGCGGCGTCTGCCAGACGAACGAAATATTTCCCGCCGGGCCGGA  >  NC_002947/5472734‑5473026

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