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,4580CT98.3% 202.8 / NA 60P16L (CCC→CTC) cofthiamine pyrimidine pyrophosphate hydrolase, HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (33/26);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCCCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCG  >  NC_000913/467323‑467591
                                                                                                                                       |                                                                                                                                     
aaCAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGa                                                                                                                                    >  2:28468/1‑139 (MQ=255)
aaCAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGa                                                                                                                                    >  1:853470/1‑139 (MQ=255)
 aCAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGAc                                                                                                                                   <  1:1096824/139‑1 (MQ=255)
               cTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTga                                                                                                                     >  1:315490/1‑139 (MQ=255)
               cTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTg                                                                                                                      >  1:1094003/1‑138 (MQ=255)
               cTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTg                                                                                                                      <  2:1094003/138‑1 (MQ=255)
               cTATTCACAATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCAGCATTTAGGTga                                                                                                                     >  2:720587/1‑139 (MQ=255)
                     aCTATTACTTTATTATAAAGAAGGTGATACTCAATCACTATTAACCGTGTCAAAGGGTGAAAAAAGAAATGCCTCTTCTGGCAGAATTTGATATGGATGGCATTTTATTGATGCCCGACAATCTTTTAGGTGAAAAAcc                                                                                                              <  1:475582/139‑1 (MQ=255)
                      cTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAAcc                                                                                                              <  1:28468/139‑1 (MQ=255)
                               ccGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACt                                                                                                     >  2:511825/1‑139 (MQ=255)
                                 gTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACttt                                                                                                   >  1:985211/1‑139 (MQ=255)
                                          aGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACt                                                                                          >  1:946750/1‑139 (MQ=255)
                                            gTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACca                                                                                                                                 >  2:15132/1‑98 (MQ=255)
                                            gTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACca                                                                                                                                 <  1:15132/98‑1 (MQ=255)
                                            gTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGc                                                                                        >  1:475130/1‑139 (MQ=255)
                                            gTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGc                                                                                        >  2:614382/1‑139 (MQ=255)
                                            gTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGc                                                                                        >  2:462940/1‑139 (MQ=255)
                                               aTACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTg                                                                                     <  1:18185/139‑1 (MQ=255)
                                               aTACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTg                                                                                     >  1:85840/1‑139 (MQ=255)
                                                tACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGa                                                                                    >  2:181596/1‑139 (MQ=255)
                                                 aCTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGaa                                                                                   >  1:841125/1‑139 (MQ=255)
                                                 aCTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGaa                                                                                   >  1:68592/1‑139 (MQ=255)
                                                 aCTCAATCACCATTAACCGGGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGaa                                                                                   <  1:720587/139‑1 (MQ=255)
                                                   tCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAAcg                                                                                 <  2:985211/139‑1 (MQ=255)
                                                      aTCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGgc                                                                                                <  1:362479/121‑1 (MQ=255)
                                                      aTCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGgc                                                                                                >  2:362479/1‑121 (MQ=255)
                                                      aTCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGa                                                                              >  1:759882/1‑139 (MQ=255)
                                                        cACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACa                                                                            <  1:285570/139‑1 (MQ=255)
                                                               aaCCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCt                                                                     <  1:418850/139‑1 (MQ=255)
                                                               aaCCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCt                                                                     >  2:418850/1‑139 (MQ=255)
                                                                    tgtCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGaaaa                                                                                                                <  1:810050/91‑1 (MQ=255)
                                                                    tgtCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGACATGGATGGCCCTTTATTGATGCTCGACCATCATTTAGGTGAGaaaa                                                                                                                >  2:810050/1‑91 (MQ=255)
                                                                     gtCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGAATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACttt                                                               >  2:705731/1‑139 (MQ=255)
                                                                              tGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAAcc                                                                                                              >  2:154766/1‑83 (MQ=255)
                                                                              tGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAAcc                                                                                                              <  1:154766/83‑1 (MQ=255)
                                                                              tGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACggg                                                      >  2:747446/1‑139 (MQ=255)
                                                                                 agaaagaaATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCg                                                   >  2:759726/1‑139 (MQ=255)
                                                                                        aaTGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTTAACGCGACATTACCCTCACTTTTGCCACgggg                                                     <  1:659140/130‑1 (MQ=255)
                                                                                        aaTGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGATGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACgggg                                                     >  2:659140/1‑130 (MQ=255)
                                                                                              tCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGgag                                      <  2:759882/139‑1 (MQ=255)
                                                                                                 tCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGAc                                                                                                                                   <  1:1025690/43‑1 (MQ=255)
                                                                                                 tCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGAc                                                                                                                                   >  2:1025690/1‑43 (MQ=255)
                                                                                                 tCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATg                                   <  1:747446/139‑1 (MQ=255)
                                                                                                 tCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATg                                   <  1:462940/139‑1 (MQ=255)
                                                                                                          aTTTGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATAtt                          <  2:68592/139‑1 (MQ=255)
                                                                                                           tttGATATGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATtc                         >  1:814824/1‑139 (MQ=255)
                                                                                                               atatGGATGGCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCgg                     <  2:841125/139‑1 (MQ=255)
                                                                                                                 atggatggCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCgggg                   <  1:791530/139‑1 (MQ=255)
                                                                                                                    gatggCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCt                                                                     >  2:536327/1‑86 (MQ=255)
                                                                                                                    gatggCACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCt                                                                     <  1:536327/86‑1 (MQ=255)
                                                                                                                         cACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCctct                                                                                                          >  1:665373/1‑44 (MQ=255)
                                                                                                                         cACTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCctct                                                                                                          <  2:665373/44‑1 (MQ=255)
                                                                                                                          aCTTTATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGc          >  1:1038709/1‑139 (MQ=255)
                                                                                                                              tATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGa      >  2:143365/1‑139 (MQ=255)
                                                                                                                              tATTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGa      <  1:511825/139‑1 (MQ=255)
                                                                                                                               aTTGATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGAt     >  1:416644/1‑139 (MQ=255)
                                                                                                                                  gATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCg  <  2:1038709/139‑1 (MQ=255)
                                                                                                                                  gATGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCg  <  2:946750/139‑1 (MQ=255)
                                                                                                                                   aTGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACa                                                                            >  1:914595/1‑64 (MQ=255)
                                                                                                                                   aTGCTCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACa                                                                            <  2:914595/64‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
AACAATATTATTTAACTATTCACTATTACTTCCGTATATATCAGGTGATACTCAATCACCATTAACCGTGTCACAGAGTGGAGAAAGAAATGGCTCGTCTGGCAGCATTTGATATGGATGGCACTTTATTGATGCCCGACCATCATTTAGGTGAGAAAACCCTCTCTACTTTGGCGCGACTGCGTGAACGCGACATTACCCTCACTTTTGCCACGGGGCGTCATGCGCTGGAGATGCAGCATATTCTCGGGGCGCTATCGCTGGATGCG  >  NC_000913/467323‑467591

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