Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 467,839 T→C 74.8% V143A (GTC→GCC)  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,8390TC74.8% 113.6 / 33.3 67V143A (GTC→GCC) cofthiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base T (11/6);  new base C (24/26);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTT  >  NC_000913/467743‑467933
                                                                                                |                                                                                              
ttACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCAc                                                                                             <  2:1424784/100‑1 (MQ=255)
   ccGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCaa                                                                                          >  2:1558594/1‑100 (MQ=255)
   ccGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCaa                                                                                          >  2:510214/1‑100 (MQ=255)
      ggAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGAt                                                                                       >  1:1824788/1‑100 (MQ=255)
      ggAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGAt                                                                                       >  1:1824802/1‑100 (MQ=255)
       gAAAGCGAGCCCTGCGTTGTTGCAGGCTTTTGTCTATAGCGGGTTTCGTTACCAGATACTCGATGTCCAAAATATGCCACTGGGCAGCGCCACCAAGATc                                                                                      <  1:2268025/100‑1 (MQ=255)
       gAAAGAGATCCCTGCGTTGTTGTAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATc                                                                                      <  2:2417032/100‑1 (MQ=255)
         aaGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAGTGCCACTACGCAGCGCCACCCAGATCTg                                                                                    >  1:1781590/1‑100 (MQ=255)
           gagaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCt                                                                                  <  2:1553344/100‑1 (MQ=255)
             gaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTTAAAAAAATGCCTCTCGGCAGCGCCACCGAGATCTGCTTc                                                                                <  1:1482222/100‑1 (MQ=255)
              aTCCCTGCGTTTTGGCAGGCATTCGTTTTTAGCTTTTTTCGTTATCAGATAATGGATGTTTAGAATATGCCTCGCGGCAGCGCCACCAAGATCTGCTTCt                                                                               <  1:1390728/100‑1 (MQ=255)
                cccTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCtg                                                                              >  1:362868/1‑99 (MQ=255)
                 ccTGCGTTGTTGCAGGCATATGTCTATCGCGGTTTTCGTTATCAGATAATCAATGTAAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCtgtg                                                                            <  1:2225942/100‑1 (MQ=255)
                   tGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc                                                                          <  1:1000026/100‑1 (MQ=255)
                    gCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCg                                                                         >  1:1352615/1‑100 (MQ=255)
                        tgttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCa                                                                     <  2:307765/100‑1 (MQ=255)
                        tgttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCa                                                                     <  1:254512/100‑1 (MQ=255)
                          ttgCAGGCATTTGGCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCTCTCGGCAGCTCCACCAAGATCTGCTTCTGTGGCGATCacg                                                                   <  1:1181221/100‑1 (MQ=255)
                            gCAGGCGTTTGTATATAGCGGTTTTCGTTATCAGATAATCGATGTCAAACACCTGCCACTCGGCAACGCCACTGAGATCTCCtt                                                                                 >  1:2508845/1‑84 (MQ=255)
                            gCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCacgac                                                                 <  2:2088542/100‑1 (MQ=255)
                            gCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCacgac                                                                 <  2:680219/100‑1 (MQ=255)
                               ggCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCCCGACGAt                                                              >  1:2508070/1‑100 (MQ=255)
                                 cATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGCTGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCacgcc                                                                 >  1:2195233/1‑93 (MQ=255)
                                   tttGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTa                                                          >  2:320158/1‑100 (MQ=255)
                                         tataGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCt                                                    >  1:1541776/1‑100 (MQ=255)
                                           taGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCAGCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTg                                                  >  1:1653666/1‑100 (MQ=255)
                                            aGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGc                                                 >  2:1072614/1‑100 (MQ=255)
                                             gCGGTTTTCGTTATCAGATAATCGCTGTCAAAAAAATGCCGCTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCa                                                <  2:2415281/100‑1 (MQ=255)
                                             gCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCa                                                <  1:2565363/100‑1 (MQ=255)
                                              cGGTTTTCGTTATCATATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAg                                               <  2:2557502/100‑1 (MQ=255)
                                              cGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAg                                               <  2:1286154/100‑1 (MQ=255)
                                               ggTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGa                                              >  1:1037798/1‑100 (MQ=255)
                                                gTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGAt                                             >  1:2459942/1‑100 (MQ=255)
                                                gTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGAt                                             <  2:1541785/100‑1 (MQ=255)
                                                  tttCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATcc                                           <  2:1037801/100‑1 (MQ=255)
                                                     cGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGc                                        <  2:587667/100‑1 (MQ=255)
                                                     cGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGc                                        >  2:512997/1‑100 (MQ=255)
                                                     cGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGCCGATCACGACGATCTTACACGCTTGCAGATCCAGc                                        >  1:421556/1‑100 (MQ=255)
                                                      gTTCTCAGATAAAAGAAATAAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAg                                         <  1:2130401/98‑1 (MQ=255)
                                                        tATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCtat                                     >  2:1453248/1‑100 (MQ=255)
                                                        tATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCtat                                     <  1:809603/100‑1 (MQ=255)
                                                          tCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATAc                                   >  2:1332622/1‑100 (MQ=255)
                                                            aaaTAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCATGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGa                                 >  2:1750219/3‑100 (MQ=255)
                                                            aGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACg                                  >  2:1840948/1‑99 (MQ=255)
                                                              aTAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAg                               <  2:372954/100‑1 (MQ=255)
                                                                      tGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGc                       >  2:17352/1‑100 (MQ=255)
                                                                       gTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCg                      >  2:1909029/1‑100 (MQ=255)
                                                                         cAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCg                      >  1:2175085/1‑98 (MQ=255)
                                                                         cAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCTTTAGGc                       >  1:709712/1‑97 (MQ=255)
                                                                         cAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAg                    >  2:389838/1‑100 (MQ=255)
                                                                         cAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCAATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAg                         >  1:2354179/1‑95 (MQ=255)
                                                                          aaaaaaaTGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTCGACGCTTGCAGTTTCAGCTATCCGAAGCATTAGTCg                      >  1:1254151/1‑97 (MQ=255)
                                                                           aaaaaaTGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCg                  <  2:1722738/100‑1 (MQ=255)
                                                                           aaaaaaTGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCg                  >  2:2243708/1‑100 (MQ=255)
                                                                             aaaaTGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTg                <  1:470234/100‑1 (MQ=255)
                                                                              aaaTGTCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGc               <  2:582696/100‑1 (MQ=255)
                                                                               aaTGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGca              <  1:2261025/100‑1 (MQ=255)
                                                                                 tGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGcaca            <  1:2199798/100‑1 (MQ=255)
                                                                                   ccACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGTACAt           >  1:2216179/1‑99 (MQ=255)
                                                                                   ccACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACAtt          <  1:2541430/100‑1 (MQ=255)
                                                                                   ccACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACAtt          >  2:361088/1‑100 (MQ=255)
                                                                                   ccACTCGGCAGCGCCACCAAGATCTGCTTCTGTGACGATCACGACGAACTAACACGCTTGCAGCTCCAGCTCTGCGATGCATTAGCCGAGCGCGCACAt           >  1:543122/1‑99 (MQ=255)
                                                                                     aCTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTtg        <  2:1195898/100‑1 (MQ=255)
                                                                                     aCTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTtg        >  1:1970306/1‑100 (MQ=255)
                                                                                           cAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGtttt  <  2:347196/100‑1 (MQ=255)
                                                                                           cAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATGCAGCTATACGAAGGATTATGCGAGCGGGCACATGTGTGtttt  >  1:2524768/1‑100 (MQ=255)
                                                                                           cAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGtttt  <  2:1978842/100‑1 (MQ=255)
                                                                                                |                                                                                              
TTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTT  >  NC_000913/467743‑467933

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: