Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,806,674 G→A 8.0% G406G (GGC→GGT)  PP_4234 ← bifunctional aldehyde oxidase/xanthine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,806,6740GA 8.0% 143.2 / 2.3 50G406G (GGC→GGT) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (18/28);  new base A (2/2);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTAACCGCGCATGGAGCCGGCCTCGACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCTTGATGCCGCTCTGGAACTGCTCGTAGCGGTCGTTGGCCAGGCGGATCGGCA  >  NC_002947/4806536‑4806819
                                                                                                                                           |                                                                                                                                                 
gTAACCGCGCATGGAGCCGGCCTCGACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCg                                                                                                                                          >  4:352047/1‑148 (MQ=255)
gTAACCGCGCATGGAGCCGGCCTCGACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCg                                                                                                                                          >  3:360050/1‑148 (MQ=255)
     cgcgCATGGAGCCGGCCTCGACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCt                                                                                                                                     >  6:310256/1‑148 (MQ=255)
      gcgcATGGAGCCGGCCTCGACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCACCGCCGTCGCAGCTCa                                                                                                                                   >  8:203945/1‑149 (MQ=255)
      gcgcATGGAGCCGGCCTCGACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCACCGCCGTCGCAGCTCa                                                                                                                                   >  8:209812/1‑149 (MQ=255)
           tGGAGCCGGCCTCGACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCg                                                                                                                              <  4:67671/149‑1 (MQ=255)
              aGCCGGCCTCGCCGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCAGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGgcg                                                                                                                           <  4:278896/149‑1 (MQ=255)
               gCCGGCCTCGACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGgcgc                                                                                                                          <  8:147399/149‑1 (MQ=255)
                        gACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCt                                                                                                                 <  3:213658/149‑1 (MQ=255)
                        gACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCt                                                                                                                 >  5:203128/1‑149 (MQ=255)
                         aCGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCtt                                                                                                                <  1:401682/149‑1 (MQ=255)
                         aCGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCACCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCtt                                                                                                                <  7:209812/149‑1 (MQ=255)
                         aCGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCACCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCtt                                                                                                                <  7:203945/149‑1 (MQ=255)
                           gCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGCCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGCTGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGa                                                                                                              <  6:298098/149‑1 (MQ=255)
                            ccACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGaa                                                                                                             <  5:180149/149‑1 (MQ=255)
                              aCGGGAGTAGTAGGCGGTTACCGTCAGGTCGTTCTGT‑GGCATTTTGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAg                                                                                                                                       <  8:107095/121‑1 (MQ=255)
                              aCGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAg                                                                                                                                       >  7:107095/1‑121 (MQ=255)
                              aCGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGc                                                                                                           >  3:837/1‑149 (MQ=255)
                                   agtagtAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCtgttg                                                                                                      >  5:257675/1‑149 (MQ=255)
                                       gtagGCGGTTACCGACAGGTCGTTCTGTGGGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTcc                                                                                                   >  1:230460/1‑149 (MQ=255)
                                          ggCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCcttct                                                                                               <  3:352047/149‑1 (MQ=255)
                                          ggCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCcttc                                                                                                <  1:198559/148‑1 (MQ=255)
                                          ggCGGTTACCGACAGGTCGTTCTAT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCcttc                                                                                                <  3:311034/148‑1 (MQ=255)
                                             ggTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGt                                                                                            <  3:168825/149‑1 (MQ=255)
                                                  ccGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCAc                                                                                        <  4:360050/148‑1 (MQ=255)
                                                   cGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACTGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTTTTGTCCTTCTTGTCCACCg                                                                                      <  3:153423/149‑1 (MQ=255)
                                                   cGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCg                                                                                      >  3:379347/1‑149 (MQ=255)
                                                       aGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAg                                                                                  >  4:368620/1‑149 (MQ=255)
                                                         gTCGTTCTGT‑GGCATGTAGTAGGTCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTAGGTGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCt                                                                                >  7:59847/1‑149 (MQ=255)
                                                          tCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTg                                                                               >  5:401535/1‑149 (MQ=255)
                                                                  t‑gGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCgg                                                                        <  6:257675/148‑1 (MQ=255)
                                                                     gCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCcttc                                                                                                <  3:385149/122‑1 (MQ=255)
                                                                     gCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCcttc                                                                                                >  4:385149/1‑122 (MQ=255)
                                                                     gCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGGCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGat                                                                      <  5:310256/148‑1 (MQ=255)
                                                                          tcggacaTCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGCTGATGCGCCCGCCGCCGTCGCAGCTCATGGCGGCGCGGAAGATCTTGCAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGaa                                                                <  8:182866/143‑1 (MQ=255)
                                                                            gtagATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAgg                                                              <  2:230460/149‑1 (MQ=255)
                                                                               gATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCAt                                                           >  5:314065/1‑149 (MQ=255)
                                                                                     cTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCt                                                     <  6:401535/149‑1 (MQ=255)
                                                                                       ctgtagcagTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCTTg                                                   <  1:349184/145‑1 (MQ=255)
                                                                                        ggCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCtt                                                    <  4:837/147‑1 (MQ=255)
                                                                                        ggCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCtt                                                    <  6:314065/147‑1 (MQ=255)
                                                                                             cagTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCctt                                                                                                 <  4:238771/97‑1 (MQ=255)
                                                                                             cagTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCctt                                                                                                 >  3:238771/1‑97 (MQ=255)
                                                                                                     cccACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCTCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCTTGATGCCGCTCTGGaa                                     <  6:203128/149‑1 (MQ=255)
                                                                                                     cccACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCTTGATGCCGCTCTGGaa                                     <  8:273933/149‑1 (MQ=255)
                                                                                                     cccACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCTTGATGCCGCTCTGGaa                                     >  7:273933/1‑149 (MQ=255)
                                                                                                               ccACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGTGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCTTGATGCCGCTCTGGAACTGCTCGTAg                           <  8:59847/149‑1 (MQ=255)
                                                                                                                cACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCTTGATGCCGCTCTgg                                       <  8:166067/136‑1 (MQ=255)
                                                                                                                cACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCTTGATGCCGCTCTgg                                       >  7:166067/1‑136 (MQ=255)
                                                                                                                                 tGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCTTGATGCCGCTCTGGAACTGCTCGTAGCGGTCGTTGGCCAGGCg          >  5:224562/1‑148 (MQ=255)
                                                                                                                                        cccgccGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCGGTTGTCCTGCTTGTCCACCGCCAGCTGGTAGCGGATATCGAACGCATGGCGCTTGATGCCGCTCTGGAACTGCTCGTATCGGTCGTTTGCCAGGCGGATAGGCa  >  7:300128/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
GTAACCGCGCATGGAGCCGGCCTCGACGCCACGGGAGTAGTAGGCGGTTACCGACAGGTCGTTCTGT‑GGCATGTAGTAGATCGACTGGGCAGCAGTGGCGCCCACGGCGGCCACCGACGGGCTGTAGTTGATGCGCCCGCCGCCGTCGCAGCTCATGTCGGCGCGGAAGATCTTGAAGCTGTTGTCCTTCTTGTCCACCGCCAGCTGGTAGCGGATATCGAAGGCATGGCGCTTGATGCCGCTCTGGAACTGCTCGTAGCGGTCGTTGGCCAGGCGGATCGGCA  >  NC_002947/4806536‑4806819

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