Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,712,952 T→A 8.0% intergenic (+35/‑105) adk → / → tsaB adenylate kinase/endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,712,9520TA 8.0% 143.5 / 2.5 50intergenic (+35/‑105)adk/tsaBadenylate kinase/endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base A (2/2);  total (25/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 = 9.69e-01

TACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGT  >  NC_002947/1712810‑1713098
                                                                                                                                              |                                                                                                                                                  
tACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGc                                                                                                                                              >  1:193016/1‑149 (MQ=255)
  ccAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCggg                                                                                                                                               <  8:83746/146‑1 (MQ=255)
     gAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACTTCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGtt                                                                                                                                          <  6:70120/148‑1 (MQ=255)
     gAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGtt                                                                                                                                          >  5:70120/1‑148 (MQ=255)
       aGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTc                                                                                                                                       <  5:211864/149‑1 (MQ=255)
       aGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTc                                                                                                                                       <  2:148620/149‑1 (MQ=255)
           gTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGTAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGTTGTCga                                                                                                                                     <  1:80268/147‑2 (MQ=255)
           gTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGGAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCAAGCGGGCCGTTGTCga                                                                                                                                     >  2:80268/1‑146 (MQ=255)
              ggCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGt                                                                                                                                 <  4:205023/148‑1 (MQ=255)
                         gcggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAgc                                                                                                                     <  5:317/149‑1 (MQ=255)
                           ggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGt                                                                                                                                 >  1:242733/1‑135 (MQ=255)
                           ggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGcc                                                                                                                            <  5:110894/140‑1 (MQ=255)
                           ggcAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGcc                                                                                                                            >  6:110894/1‑140 (MQ=255)
                                  cGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGtt                                                                                                                                          >  2:239186/1‑119 (MQ=255)
                                  cGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGtt                                                                                                                                          <  1:239186/119‑1 (MQ=255)
                                              acaTCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTGGTTTGTAGGCCGCCCCGCGCCCCCGCAGGACCCAAcgcc                                                                                                >  2:236895/1‑149 (MQ=255)
                                                       gtgtCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCgctgctgct                                                                                       <  7:178927/149‑1 (MQ=255)
                                                               tCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGc                                                                                                                             >  7:160686/1‑103 (MQ=255)
                                                               tCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGc                                                                                                                             <  8:160686/103‑1 (MQ=255)
                                                                   ttGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGc                                                                           >  3:140570/1‑149 (MQ=255)
                                                                         ccATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCtttt                                                                     >  3:182393/1‑149 (MQ=255)
                                                                           aTCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTTAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCttttt                                                                    >  3:153534/1‑148 (MQ=255)
                                                                           aTCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTc                                                                   <  2:193016/149‑1 (MQ=255)
                                                                              aCCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCAt                                                                 <  7:92337/148‑1 (MQ=255)
                                                                                cGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCg                                                              >  2:207238/1‑149 (MQ=255)
                                                                                cGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCg                                                              >  8:192492/1‑149 (MQ=255)
                                                                                     aGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCt                                                         <  7:192492/149‑1 (MQ=255)
                                                                                          cTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCCCTGCAGCAAAATACAGCCGCTCTTTTTCATTAGCACCTGGAt                                                     >  4:139468/1‑148 (MQ=255)
                                                                                               aGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGt                                               <  1:236895/149‑1 (MQ=255)
                                                                                                gCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGtt                                              <  4:182393/149‑1 (MQ=255)
                                                                                                  ccTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTc                                             >  5:15678/1‑148 (MQ=255)
                                                                                                    tGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGAt                                          >  5:236972/1‑149 (MQ=255)
                                                                                                        cTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGacc                                      <  6:236972/149‑1 (MQ=255)
                                                                                                          gATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGGCGTTGTAGTTTGTAGGCCGCACAGCGCCCCCTAAGGAAACAACGCCGGTGATGCTCTATAATGCCGGTCTTTTTCATTAGCACCTTGATACCCGTTGGATGgccac                                    >  7:168654/1‑149 (MQ=255)
                                                                                                              cGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCct                                 <  1:43920/148‑1 (MQ=255)
                                                                                                                       aCGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTg                       <  4:140570/149‑1 (MQ=255)
                                                                                                                          cACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGAt                    <  3:139468/149‑1 (MQ=255)
                                                                                                                             cTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATa                   >  5:74423/1‑147 (MQ=255)
                                                                                                                               cGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTg                       >  1:176580/1‑141 (MQ=255)
                                                                                                                               cGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTg                       <  2:176580/141‑1 (MQ=255)
                                                                                                                                gACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTTAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGc               <  4:153534/148‑1 (MQ=255)
                                                                                                                                 aCAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCa             >  5:150663/1‑149 (MQ=255)
                                                                                                                                  cAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCAc            <  4:70583/149‑1 (MQ=255)
                                                                                                                                  cAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCAc            >  3:70583/1‑149 (MQ=255)
                                                                                                                                   aaCGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCAcc           >  7:210695/1‑149 (MQ=255)
                                                                                                                                      ggCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTgg                            >  1:134187/1‑129 (MQ=255)
                                                                                                                                      ggCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTgg                            <  2:134187/129‑1 (MQ=255)
                                                                                                                                      ggCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCg          >  8:196506/1‑147 (MQ=255)
                                                                                                                                        cccGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGc      >  6:213741/1‑149 (MQ=255)
                                                                                                                                            cTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGt  <  6:150663/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TACCAGAAGCTGTCGGCTGCCAATGGCGGCAAGCCGAAGTACAGCCACATCGAAGGTGTCGGCTCGGTTGAAGCCATCACCGCCAAGGTCCTGGCAGCCCTGAGCTGATCCGGCGCCTGACGCACCTCGACAACGGCCCGCTTGCGGGCCGTTGTCGTTTGTAGGCCGCACAGCGCCCCCTCAGGACACAACGCCGCTGCTGCTCTATACTGCCGCTCTTTTTCATTCGCACCTGGATACCCGTTCGATGACCACCCTGCTGGCCCTGGATACCGCCACCGAAGCCTGT  >  NC_002947/1712810‑1713098

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