Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 909,353 G→A 74.0% A566V (GCG→GTG)  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,3530GA74.0% 86.3 / 34.1 56A566V (GCG→GTG) poxBpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base A (19/22);  total (27/29)
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.94e-01

GTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCAC  >  NC_000913/909205‑909483
                                                                                                                                                    |                                                                                                                                  
gTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCg                                                                                                                                    <  3:221410/149‑1 (MQ=255)
              cAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACt                                                                                                                      >  8:153248/1‑149 (MQ=255)
               aaCCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCAGCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACtt                                                                                                                     >  1:330640/1‑149 (MQ=255)
                aCCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACtt                                                                                                                     >  6:182070/1‑148 (MQ=255)
                aCCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTtc                                                                                                                    <  6:298615/149‑1 (MQ=255)
                               atcGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGc                                                                                                      <  2:194768/148‑1 (MQ=255)
                                  gCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCtgat                                                                                                  <  5:182070/149‑1 (MQ=255)
                                    ttCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCGCTgatgat                                                                                               >  5:336176/1‑149 (MQ=255)
                                       gAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCGCTGATGATTgc                                                                                            <  6:336176/149‑1 (MQ=255)
                                              gcgcGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCATTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTAATCACCGCGTCCGCTGATGATTCCGCgcag                                                                                      >  6:217236/1‑149 (MQ=255)
                                                 cgGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCat                                                                                   <  8:364369/149‑1 (MQ=255)
                                                    cTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCAtata                                                                                >  1:289126/1‑149 (MQ=255)
                                                        gTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAg                                                                             <  1:12095/148‑1 (MQ=255)
                                                            ggTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGa                                                                        <  2:289126/149‑1 (MQ=255)
                                                            ggTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGAAGATTGCGCGCAGCATATACAGGCTg                                                                         >  7:356859/1‑148 (MQ=255)
                                                               aaCGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCtgatg                                                                                                 <  3:273775/121‑1 (MQ=255)
                                                               aaCGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCtgatg                                                                                                 >  4:273775/1‑121 (MQ=255)
                                                                    tATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCTCCCTTTTTACCTTAGCCATTTTGTTTTCAACAGTTCGATCACTTCATCACCGCGTCCGCTTATGATTGCGCGCAGCATATACAGGCTGAAACCTTTg                                                                <  8:138936/149‑1 (MQ=255)
                                                                    tATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTg                                                                <  7:153248/149‑1 (MQ=255)
                                                                    tATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTg                                                                <  5:131408/149‑1 (MQ=255)
                                                                        aCTGCGTAAATCAATCATGGCATGTCCTTATTATTACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATAACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCt                                                            >  6:198267/1‑149 (MQ=255)
                                                                         cTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCatat                                                                                 <  1:219049/127‑1 (MQ=255)
                                                                         cTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCatat                                                                                 >  2:219049/1‑127 (MQ=255)
                                                                               aaaTCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATc                                                                                                                         <  7:24067/81‑1 (MQ=255)
                                                                               aaaTCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATc                                                                                                                         >  8:24067/1‑81 (MQ=255)
                                                                               aaaTCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCAAGTTTGTTTTCTCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAg                                                     >  7:295285/1‑149 (MQ=255)
                                                                                 atcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACt                                                                                                                      <  1:313468/82‑1 (MQ=255)
                                                                                 atcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACt                                                                                                                      >  2:313468/1‑82 (MQ=255)
                                                                                 atcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGtt                                                   >  3:81980/1‑149 (MQ=255)
                                                                                   caatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGttt                                                  >  7:338907/1‑148 (MQ=255)
                                                                                      tcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATc                                              >  4:242938/1‑149 (MQ=255)
                                                                                         tGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGttt                                                  <  8:148699/142‑1 (MQ=255)
                                                                                         tGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGttt                                                  >  7:148699/1‑142 (MQ=255)
                                                                                         tGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGAtt                                                                                              <  4:179752/98‑1 (MQ=255)
                                                                                         tGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGAtt                                                                                              >  3:179752/1‑98 (MQ=255)
                                                                                             aTGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTg                                       >  7:313210/1‑149 (MQ=255)
                                                                                             aTGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTg                                       <  3:242938/149‑1 (MQ=255)
                                                                                                tCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCGGATGATTGCGCGCAGCATATACAGGCTTAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGaa                                    <  5:185892/149‑1 (MQ=255)
                                                                                                  cTTATTATGACGGGAAATGCCACCCTTTTTACATTAGCCAGTTTGTTTTCACCAGTTCGAACACTTAAACACCGCGACCGCTGATGATTCCGCGCAGCATAAACACGCTGACACCTTTGGCCTGTTCGAGTTTGATCTGCGTTGGAATg                                  >  4:310074/1‑149 (MQ=255)
                                                                                                   ttattaTGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATgg                                 >  4:343557/1‑149 (MQ=255)
                                                                                                     attatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCgca                                                                                       >  4:294781/1‑93 (MQ=255)
                                                                                                     attatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCgca                                                                                       <  3:294781/93‑1 (MQ=255)
                                                                                                     attatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGc                                <  8:356859/148‑1 (MQ=255)
                                                                                                           aCGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAACCCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAActc                          <  4:7723/148‑1 (MQ=255)
                                                                                                            cGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTg                                       <  2:336336/134‑1 (MQ=255)
                                                                                                                aaaTGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCttt                    <  8:313210/149‑1 (MQ=255)
                                                                                                                 aaTGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTg                   <  4:81980/149‑1 (MQ=255)
                                                                                                                     ccACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAAt                                   >  7:49715/1‑129 (MQ=255)
                                                                                                                     ccACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAAt                                   >  7:55680/1‑129 (MQ=255)
                                                                                                                     ccACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAAt                                   <  8:49715/129‑1 (MQ=255)
                                                                                                                     ccACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAAt                                   <  8:55680/129‑1 (MQ=255)
                                                                                                                      cACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGac              >  5:279790/1‑149 (MQ=255)
                                                                                                                      cACCCTTATTACCTTAGACAGTTTCTTTTCACCAGTTCGATCACTTCATCACCGCGTCAGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGac              >  6:358898/1‑149 (MQ=255)
                                                                                                                         ccTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACa                                                                              <  2:386823/82‑1 (MQ=255)
                                                                                                                                  cTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATccac  <  3:343557/149‑1 (MQ=255)
                                                                                                                                     aGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACtcttct                      <  2:86379/126‑1 (MQ=255)
                                                                                                                                     aGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACtcttct                      >  1:86379/1‑126 (MQ=255)
                                                                                                                                                    |                                                                                                                                  
GTAAACGTCGTCCCCAACCGGGGCGGCCATCATCGCTTCGAGCATGGCGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCAC  >  NC_000913/909205‑909483

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