Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,167,011 C→T 7.4% M155I (ATG→ATA)  creA ← creatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,167,0110CT 7.4% 154.5 / 2.3 54M155I (ATG→ATA) creAcreatinase
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base T (2/2);  total (24/30)
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.95e-01

CGCACCATCGCCTGGGTGGCATGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCAACGCCTGCTGGATGGCGACGAAATAGT  >  NC_002947/4166865‑4167159
                                                                                                                                                  |                                                                                                                                                    
cgcACCATCGCCTGGGTGGCATGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCAt                                                                                                                                                    <  2:125870/149‑1 (MQ=255)
  cACCATCGCCTGGGTGGCATGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATAc                                                                                                                                                  <  2:22645/149‑1 (MQ=255)
    ccATCGCCTGGGTGGCATGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGTCCACCTGTGCACCGTGACGGATAATTGCGTGCTCCTCGGCCGACTTGATCATGCGTATgcgc                                                                                                                                                <  1:38642/149‑1 (MQ=255)
     cATCGCCTGGGTGGCATGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGTATGCGCa                                                                                                                                               >  4:168060/1‑149 (MQ=255)
         gCCTGGGTGGCATGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGTATGCGCATga                                                                                                                                            <  3:168060/148‑2 (MQ=255)
                     tGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGAt                                                                                                                               >  1:300067/1‑149 (MQ=255)
                         gTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTcc                                                                                                                           >  5:256488/1‑149 (MQ=255)
                            ccACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTccac                                                                                                                         <  5:191427/148‑1 (MQ=255)
                             cACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTccac                                                                                                                         >  7:103799/1‑147 (MQ=255)
                              aCCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAg                                                                                                                      >  4:16918/1‑149 (MQ=255)
                              aCCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAg                                                                                                                      >  2:51408/1‑149 (MQ=255)
                                cTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCAt                                                                                                                                              >  6:85049/1‑123 (MQ=255)
                                   gTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTc                                                                                                                   <  2:107523/147‑1 (MQ=255)
                                         gggCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGGCTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTg                                                                                                            <  8:265005/148‑1 (MQ=255)
                                         gggCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTgg                                                                                                           <  3:94766/149‑1 (MQ=255)
                                         gggCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTgg                                                                                                           <  7:17365/149‑1 (MQ=255)
                                         gggCAACTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTg                                                                                                            <  2:244905/148‑1 (MQ=255)
                                          ggCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTggg                                                                                                          <  2:300067/149‑1 (MQ=255)
                                           gCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTgggg                                                                                                         >  8:231312/1‑149 (MQ=255)
                                                  cTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCa                                                                                                    <  6:256488/147‑1 (MQ=255)
                                                  cTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAgg                                                                                                  >  4:44240/1‑149 (MQ=255)
                                                      cGCAGGTCTTCGTCCACTGCCGACCCGCCGATGTCGGCCACCTGTGCACCGTGTCGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACTCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAggcggc                                                                                              <  8:305236/149‑1 (MQ=255)
                                                       tcAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAggcggc                                                                                              >  6:20185/2‑148 (MQ=255)
                                                       gCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGTGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAggcggc                                                                                              >  8:145701/1‑148 (MQ=255)
                                                        cAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGTGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGa                                                                                            <  7:145701/149‑1 (MQ=255)
                                                         aGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTACACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAg                                                                                           >  1:193509/1‑149 (MQ=255)
                                                          gggCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGt                                                                                          <  5:274929/149‑1 (MQ=255)
                                                                cGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGcc                                                                                     <  1:51408/148‑1 (MQ=255)
                                                                  aCCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCTGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTg                                                                                  <  1:165442/149‑1 (MQ=255)
                                                                   ccACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTgg                                                                                 <  5:233616/149‑1 (MQ=255)
                                                                     aCTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGtt                                                                               <  7:231312/149‑1 (MQ=255)
                                                                       tGCCGCCCCGCCGATGTCGGCCACCTATGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCt                                                                             >  3:178745/1‑149 (MQ=255)
                                                                              ccgccgATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGtt                                                                      >  4:45878/1‑149 (MQ=255)
                                                                                 ccgATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAg                                                                   >  2:320926/1‑149 (MQ=255)
                                                                                  cgATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGATCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAgg                                                                  >  5:220459/1‑149 (MQ=255)
                                                                                    aTGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAggtg                                                                >  6:242826/1‑149 (MQ=255)
                                                                                      gTCGGCCACCTATGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAggtgg                                                               <  4:178745/148‑1 (MQ=255)
                                                                                       tCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAggtggt                                                              <  3:16918/148‑1 (MQ=255)
                                                                                       tCGGCCACCTGTACACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGTCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTc                                                             <  2:193509/149‑1 (MQ=255)
                                                                                            cACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACt                                                        <  5:242826/149‑1 (MQ=255)
                                                                                                 gtgCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACg                                                   >  7:34512/1‑149 (MQ=255)
                                                                                                      ccGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGAt                                              <  8:103799/149‑1 (MQ=255)
                                                                                                          gACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATgcgc                                          >  1:171207/1‑149 (MQ=255)
                                                                                                                   gTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGCGGTCGTACTCCACGCCGATGCGCGAGGCCAGc                                 >  3:243574/1‑149 (MQ=255)
                                                                                                                     gacTGCTCCTCTTCCGACTTGATCATGCGAATACGCATGAACGTCGCAGAGGTGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTTCAGGTTCAGGTGGGCGTACTCCACGCCGATGCGCGAGGCCAGCg                                <  2:212929/146‑1 (MQ=255)
                                                                                                                     gCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCgg                               <  5:20185/149‑1 (MQ=255)
                                                                                                                       cTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCa                             <  3:44240/149‑1 (MQ=255)
                                                                                                                          ctcctcGGCCGACTTGATCATGCGCATACGCATGCACg                                                                                                                                         <  2:152905/38‑1 (MQ=255)
                                                                                                                          ctcctcGGCCGACTTGATCATGCGCATACGCATGCACg                                                                                                                                         >  1:152905/1‑38 (MQ=255)
                                                                                                                          ctcctcGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCAAc                           <  2:171207/148‑1 (MQ=255)
                                                                                                                              tcGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCAACGcc                        <  7:196207/147‑1 (MQ=255)
                                                                                                                              tcGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCAACGcc                        >  8:196207/1‑147 (MQ=255)
                                                                                                                                         gtgcATGCGTATGCGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCAACGCCTGCTGGATGGcga           >  6:268228/4‑149 (MQ=255)
                                                                                                                                         gATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCAACGCCTGCTGGATGGcga           >  7:209011/1‑149 (MQ=255)
                                                                                                                                                cgcATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCAACGCCTGCTGGATGGCGACGaaa      >  3:152939/1‑147 (MQ=255)
                                                                                                                                                cgcATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCAACGCCTGCTGGATGGCGACGAAATa    <  8:34512/149‑1 (MQ=255)
                                                                                                                                                  cATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCAACGCCTGCTGGATGGCGACGAAATAGt  <  8:209011/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCACCATCGCCTGGGTGGCATGCAGTGCCACCTCGTACTCGGGCACCTGCTCACGCAGGGCTTCGACCACTGCCGCCCCGCCGATGTCGGCCACCTGTGCACCGTGACGGATAAGTGCCTGCTCCTCGGCCGACTTGATCATGCGCATACGCATGCACGGCGCAGCGATGTCCACCAGCTCGGCCTTGGGGTAGCAGGCGGCGAGTTTGCCGTGGTTCTGCAGGTTCAGGTGGTCGTACTCCACGCCGATGCGCGAGGCCAGCGGCAACGCCTGCTGGATGGCGACGAAATAGT  >  NC_002947/4166865‑4167159

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