Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 909,753 Δ1 bp 7.3% coding (1297/1719 nt) poxB ← pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913909,7530C. 7.3% 183.5 / 8.3 55coding (1297/1719 nt)poxBpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (31/20);  new base . (2/2);  total (33/22)
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.70e-01

CGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAGATAACGTGCCGCCCACACCGT  >  NC_000913/909615‑909900
                                                                                                                                          |                                                                                                                                                   
cGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACAt                                                                                                                                            <  1:241340/148‑1 (MQ=255)
cGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACAt                                                                                                                                            <  3:312284/148‑1 (MQ=255)
     gtcaAATAGCCACCAGCTTTCATCTCCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGa                                                                                                                                       <  7:200493/148‑1 (MQ=255)
             gCCACCAGCTTTCATCTCCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGa                                                                                                                               <  4:69774/148‑1 (MQ=255)
                accaGCTTTCATCTCCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGt                                                                                                                            >  8:161793/1‑148 (MQ=255)
                 ccaGCTTTCATCTCCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGtt                                                                                                                           >  2:109456/1‑148 (MQ=255)
                  caGCTTTCATCTCCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTc                                                                                                                          >  8:254915/1‑148 (MQ=255)
                            ctcCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTc                                                                                                                >  2:236497/1‑148 (MQ=255)
                             tcCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCg                                                                                                               >  3:265507/1‑148 (MQ=255)
                              ccATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGAAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGc                                                                                                              >  7:90812/1‑148 (MQ=255)
                               cATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGt                                                                                                                 <  7:89337/144‑1 (MQ=255)
                                  cGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTgc                                                                                                          >  8:224800/1‑148 (MQ=255)
                                  cGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTg                                                                                                           <  2:348694/147‑1 (MQ=255)
                                  cGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTg                                                                                                           >  1:348694/1‑147 (MQ=255)
                                     cacaAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACAt                                                                                                                                            >  1:49505/1‑111 (MQ=255)
                                     cacaAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACAt                                                                                                                                            <  2:49505/111‑1 (MQ=255)
                                       caAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCAcc                                                                                                     >  6:362595/1‑148 (MQ=255)
                                         aaGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATTCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCAcc                                                                                                     >  4:214170/1‑146 (MQ=255)
                                          aGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAg                                                                                                  <  7:161793/148‑1 (MQ=255)
                                            accAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTCAGAGGAAATCGCCAATCAACATGCTAAAACCGCAATCGCCGCACATGGCGACCACCAGACGTTCTGGCAATCTCGACAGCGCACCCAgct                                                                                                >  6:349099/2‑147 (MQ=255)
                                              cAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGcc                                                                                              >  3:311925/1‑148 (MQ=255)
                                                gCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCTCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTg                                                                                            <  8:223250/148‑1 (MQ=255)
                                                     cTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCa                                                                                       >  4:339613/1‑148 (MQ=255)
                                                              aaGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTa                                                                               <  4:311925/147‑1 (MQ=255)
                                                                gacgacAATTTTCACTGGCAGTTTCATATGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGcc                                                                            <  8:90812/148‑1 (MQ=255)
                                                                    acAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTACGGCATGGCGTTAGCCATc                                                                         >  4:351154/1‑147 (MQ=255)
                                                                      aaTTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGaa                                                                      <  7:224800/148‑1 (MQ=255)
                                                                                ggCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTATGTCGCCTGCGCACCAAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTa                                                             >  5:346174/1‑147 (MQ=255)
                                                                                      ttCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAAc                                                      >  2:201354/1‑148 (MQ=255)
                                                                                      ttCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCAACTGACGTTCTCTCTCTCTCGCCTGCGCACCCACCGCCTGCGGCATGGCGTTAGCCATCGAACCCTGGTTAAACGAAc                                                      >  3:143770/1‑148 (MQ=255)
                                                                                           cTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTAGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTaa                                                  >  1:354515/1‑147 (MQ=255)
                                                                                             gCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAg                                               >  5:359610/1‑148 (MQ=255)
                                                                                              cACTACTGAGAGGAAATCGCCCATCAACATTCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTaa                                                  <  3:214170/144‑1 (MQ=255)
                                                                                              cACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAAc                                                 <  7:254915/145‑1 (MQ=255)
                                                                                                   ctGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACg                                         <  1:109456/148‑1 (MQ=255)
                                                                                                      agagGAAATCGCCCATCAACATGCTAAAACCGCAATCGCCGCAAATGGCGACCACCTGACGTTATGGCTCTGTCGCCTGCGCACCCACCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCtt                                      >  5:254763/1‑148 (MQ=255)
                                                                                                                   cATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACTCTTGCCGTTCAttttt                         <  5:362595/148‑1 (MQ=255)
                                                                                                                   cATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCAttttt                         <  3:339613/148‑1 (MQ=255)
                                                                                                                    aTCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCAtttt                          >  7:282714/1‑146 (MQ=255)
                                                                                                                      cAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCAtttt                          >  3:275717/1‑144 (MQ=255)
                                                                                                                      cAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCAtttt                          >  2:343132/1‑144 (MQ=255)
                                                                                                                      cAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCAtttt                          >  3:226780/1‑144 (MQ=255)
                                                                                                                        aCATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCAATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTTCCGTTCATTTTTAGAt                     >  4:61557/1‑147 (MQ=255)
                                                                                                                         cATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAGATaa                   <  4:143770/148‑1 (MQ=255)
                                                                                                                         cATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAGATaa                   >  2:169233/1‑148 (MQ=255)
                                                                                                                         cATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAGATaa                   <  1:236497/148‑1 (MQ=255)
                                                                                                                          attcTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAGATAAc                  >  6:18457/4‑148 (MQ=255)
                                                                                                                          aTGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGTCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAGATAAc                  >  5:11915/1‑148 (MQ=255)
                                                                                                                          aTGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGCCGCTTGCCGTTCATTTTTAGATAAc                  >  7:78930/1‑148 (MQ=255)
                                                                                                                                  aacGCCAT‑GCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAg                       >  6:49424/3‑134 (MQ=255)
                                                                                                                                  aacGCCAT‑GCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAg                       >  6:44238/3‑134 (MQ=255)
                                                                                                                                  aacGCCAT‑GCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAg                       <  5:49424/132‑1 (MQ=255)
                                                                                                                                  aacGCCAT‑GCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAg                       <  5:44238/132‑1 (MQ=255)
                                                                                                                                      ccATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAg                                               <  2:33064/107‑1 (MQ=255)
                                                                                                                                      ccATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAg                                               >  1:33064/1‑107 (MQ=255)
                                                                                                                                          cgccgcACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAGATAACGTGCCACCCACACCGt  >  7:201036/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
CGTCAGTCAAATAGCCACCAGCTTTCATCTCCATCGCCACAAAGCCCAGCACGCTGTTGTTAAAGACGACAATTTTCACTGGCAGTTTCATCTGCACTACTGAGAGGAAATCGCCCATCAACATGCTAAAACCGCCATCGCCGCACATGGCGACCACCTGACGTTCTGGCTCTGTCGCCTGCGCACCCAGCGCCTGCGGCATGGCGTTAGCCATCGAACCGTGGTTAAACGAACCTAACAGGCGACGCTTGCCGTTCATTTTTAGATAACGTGCCGCCCACACCGT  >  NC_000913/909615‑909900

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑