Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 467,458 C→T P16L (CCC→CTC)  cof → thiamine pyrimidine pyrophosphate hydrolase, HMP‑PP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913467,4580CT94.8% 179.9 / 1.5 58P16L (CCC→CTC) cofthiamine pyrimidine pyrophosphate hydrolase, HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (33/22);  total (35/23)
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.87e-01

AATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCCCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCGT  >  NC_000913/467326‑467592
                                                                                                                                    |                                                                                                                                      
aaTATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACca                                                                                                                                  <  1:284986/139‑1 (MQ=255)
       tttAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTa                                                                                                                           >  1:628700/1‑139 (MQ=255)
       tttAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTa                                                                                                                           <  1:279272/139‑1 (MQ=255)
       tttAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTa                                                                                                                           >  2:585888/1‑139 (MQ=255)
                      tcaCTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATACGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCctc                                                                                                            >  1:481691/3‑139 (MQ=255)
                         cTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCctctct                                                                                                         >  2:280104/1‑139 (MQ=255)
                            ccGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACt                                                                                                      >  1:123977/1‑139 (MQ=255)
                            ccGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTCATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACt                                                                                                      >  1:485297/1‑139 (MQ=255)
                              gTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGaaa                                                                                                                  <  1:385469/125‑1 (MQ=255)
                              gTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGaaa                                                                                                                  >  2:385469/1‑125 (MQ=255)
                                     tCAGGTGATACTCAATCACCATTAACCGAGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGa                                                                                             <  2:481691/139‑1 (MQ=255)
                                      cAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACcatca                                                                                                                               <  1:398503/104‑1 (MQ=255)
                                      cAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACcatca                                                                                                                               >  2:398503/1‑104 (MQ=255)
                                         gTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGc                                                                                         >  2:693526/1‑139 (MQ=255)
                                           gATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAAcc                                                                                                               >  1:49153/1‑115 (MQ=255)
                                           gATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAAcc                                                                                                               <  2:49153/115‑1 (MQ=255)
                                             tACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGCTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGa                                                                                     >  2:502862/1‑139 (MQ=255)
                                               cTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAAc                                                                                   >  1:62893/1‑139 (MQ=255)
                                               cTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAAc                                                                                   >  1:138307/1‑139 (MQ=255)
                                                tCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAAcg                                                                                  <  2:485297/139‑1 (MQ=255)
                                                tCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAAcg                                                                                  <  2:123977/139‑1 (MQ=255)
                                                tCAATCACCATTAACCGTGTCACAGAGTCGAGACAGAAATGGCTCGTCTGGCAGCATTTCAAATGGATGGCCCTTTATTGCCGCTCGACCATCATTTAAGTGCGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAAcg                                                                                  <  1:481826/139‑1 (MQ=255)
                                                    tCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGAc                                                                              >  2:431486/1‑139 (MQ=255)
                                                                 tgtCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACtt                                                                 >  1:275282/1‑139 (MQ=255)
                                                                  gtCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACttt                                                                >  1:343212/1‑139 (MQ=255)
                                                                  gtCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACttt                                                                >  1:641325/1‑139 (MQ=255)
                                                                  gtCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACttt                                                                >  2:629239/1‑139 (MQ=255)
                                                                  gtCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACttt                                                                >  2:217678/1‑139 (MQ=255)
                                                                  gtCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCCCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACAt                                                                            <  1:656571/127‑1 (MQ=255)
                                                                  gtCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCCCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACAt                                                                            >  2:656571/1‑127 (MQ=255)
                                                                           tGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACggg                                                       >  2:87712/1‑139 (MQ=255)
                                                                           tGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACggg                                                       <  2:628700/139‑1 (MQ=255)
                                                                           tGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACggg                                                       >  1:516420/1‑139 (MQ=255)
                                                                           tGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACggg                                                       <  1:585888/139‑1 (MQ=255)
                                                                             gagaAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGc                                                     >  2:165065/1‑139 (MQ=255)
                                                                                   gaaaTGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCCCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATg                                               >  2:182056/1‑139 (MQ=255)
                                                                                     aaTGGCTCGTCTGGCAGCATTTGATATTGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATgcg                                             >  1:581712/1‑139 (MQ=255)
                                                                                       tGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCt                                           >  2:295068/1‑139 (MQ=255)
                                                                                          cTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCctct                                                                                                           <  1:126335/72‑1 (MQ=255)
                                                                                          cTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCctct                                                                                                           >  2:126335/1‑72 (MQ=255)
                                                                                              tCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATg                                    >  2:461088/1‑139 (MQ=255)
                                                                                              tCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATg                                    <  1:280104/139‑1 (MQ=255)
                                                                                              tCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATg                                    <  2:138307/139‑1 (MQ=255)
                                                                                                 ggCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATgcag                                 <  1:629239/139‑1 (MQ=255)
                                                                                                       aTTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATAtt                           <  2:275282/139‑1 (MQ=255)
                                                                                                        tttGATAGGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATtc                          >  1:463491/1‑139 (MQ=255)
                                                                                                                ggatggCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGgcg                  <  1:461088/139‑1 (MQ=255)
                                                                                                                ggatggCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGgcg                  <  1:115880/139‑1 (MQ=255)
                                                                                                                  atggCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCt                >  1:788250/1‑139 (MQ=255)
                                                                                                                    ggCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTAt              <  1:87712/139‑1 (MQ=255)
                                                                                                                           tATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTAt              <  2:446906/132‑1 (MQ=255)
                                                                                                                           tATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTAt              >  1:446906/1‑132 (MQ=255)
                                                                                                                           tATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGa       >  2:472038/1‑139 (MQ=255)
                                                                                                                           tATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGa       >  2:182828/1‑139 (MQ=255)
                                                                                                                             ttGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATg     <  1:502862/139‑1 (MQ=255)
                                                                                                                                aTGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGAGATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATg                                    >  1:524172/1‑105 (MQ=255)
                                                                                                                                aTGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATg                                    <  2:524172/105‑1 (MQ=255)
                                                                                                                                aTGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCGt  >  2:321642/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                      
AATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCCCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCGT  >  NC_000913/467326‑467592

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