Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,805,118 T→A 12.0% K925M (AAG→ATG)  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,805,1180TA12.0% 131.4 / 8.9 50K925M (AAG→ATG) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (21/23);  new base A (3/3);  total (24/26)
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.16e-01

CGTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGC  >  NC_002947/4804974‑4805258
                                                                                                                                                |                                                                                                                                             
cGTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGc                                                                                                                                            <  5:172887/148‑1 (MQ=255)
     aGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATgg                                                                                                                                      <  6:345813/149‑1 (MQ=255)
       gTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCa                                                                                                                                    >  5:357924/1‑149 (MQ=255)
       gTAGTCGATCATCGCCAGGTCCGCGGGAACATCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCa                                                                                                                                    >  1:324251/1‑149 (MQ=255)
           tCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCAt                                                                                                                                   <  6:158993/146‑1 (MQ=255)
           tCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCAt                                                                                                                                   >  5:158993/1‑146 (MQ=255)
             gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATgcgc                                                                                                                               >  4:111000/1‑148 (MQ=255)
             gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGCCATCGGCAAGcgcg                                                                                                                              >  8:187933/1‑149 (MQ=255)
               tcatcGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCCGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGCCATCGGCAAGcgcgtg                                                                                                                            <  7:187933/149‑3 (MQ=255)
                 atcGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGCCATCGGCAAGcgcgtgc                                                                                                                           <  7:186435/148‑4 (MQ=255)
                     ccAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCAtt                                                                                                                      >  5:97088/1‑149 (MQ=255)
                       aGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGc                                                                                                                    <  1:174627/149‑1 (MQ=255)
                       aGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGc                                                                                                                    <  3:111000/149‑1 (MQ=255)
                         gTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGcc                                                                                                                   >  2:216929/1‑148 (MQ=255)
                          tCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGACCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTg                                                                                                                     <  3:173792/145‑1 (MQ=255)
                          tCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGACCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTg                                                                                                                     >  4:173792/1‑145 (MQ=255)
                           ccGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGAt                                                                                                                <  5:91751/149‑1 (MQ=255)
                           ccGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGAt                                                                                                                <  5:284466/149‑1 (MQ=255)
                                          ccGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGTCGATGGCACCGACGACCgg                                                                                                 >  1:275073/1‑149 (MQ=255)
                                          ccGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCgg                                                                                                 >  8:220498/1‑149 (MQ=255)
                                             gCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGat                                                                                              <  7:20423/149‑1 (MQ=255)
                                                    cGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATcac                                                                                        >  7:228845/1‑148 (MQ=255)
                                                      gTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGAcc                                                                                                   <  6:295334/135‑1 (MQ=255)
                                                      gTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGAcc                                                                                                   >  5:295334/1‑135 (MQ=255)
                                                         gACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACtt                                                                                  <  5:20494/148‑1 (MQ=255)
                                                         gACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTc                                                                                 >  4:156674/1‑149 (MQ=255)
                                                          aCCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTc                                                                                 <  6:357924/148‑1 (MQ=255)
                                                              ttGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGa                                                                            >  6:265492/1‑149 (MQ=255)
                                                                gAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGAt                                                                           <  1:216929/148‑1 (MQ=255)
                                                                             aGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCAACTTCGGCGATGCCCTTGGAgggg                                                              <  7:220498/149‑1 (MQ=255)
                                                                             aGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGggt                                                             >  3:127089/1‑149 (MQ=255)
                                                                              gTGGACGGTTGGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGggtc                                                            <  8:11539/149‑1 (MQ=255)
                                                                               tGGACGGTAAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGggtcg                                                           <  6:105366/149‑1 (MQ=255)
                                                                                ggACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGggtcgg                                                          >  6:223895/1‑149 (MQ=255)
                                                                                 gACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATGGCCGATGGCACCGACGACCGGGATTATCACC‑ACTTCGGCGATGCCCTTGGAGGggtcggt                                                         >  8:68872/1‑149 (MQ=255)
                                                                                      ttAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGCCATCGGCAAGcgcgtgc                                                                                                                           >  7:240650/1‑76 (MQ=37)
                                                                                      ttAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGCCATCGGCAAGcgcgtgc                                                                                                                           <  8:240650/79‑4 (MQ=37)
                                                                                               gAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGTCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAgcggc                                            <  2:275073/148‑1 (MQ=255)
                                                                                                  ggCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGgga                                        <  1:18808/149‑1 (MQ=255)
                                                                                                        ctTGATACGGGCTGGGGTGATGGGGAGATCGCGAACGCGATTGCGTATGGAATTCGCGACGGCTTTGCCGATGGCACCGACGACCGGGATGATGAAA‑ACAACGGCGATGCGCTGGGAGGGGTCGGTCGGTGCGAGCGGCGggaggggtg                                  >  6:222751/1‑145 (MQ=255)
                                                                                                            aTACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGc                               <  5:223895/148‑1 (MQ=255)
                                                                                                                      gggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGt                           >  2:16469/1‑142 (MQ=255)
                                                                                                                      gggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGt                           <  1:16469/142‑1 (MQ=255)
                                                                                                                      gggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCtgctgc                    >  4:210995/1‑149 (MQ=255)
                                                                                                                           aTGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCacac               <  6:97088/149‑1 (MQ=255)
                                                                                                                             ggggAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCg             >  3:38896/1‑149 (MQ=255)
                                                                                                                                 agaTCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCcac         <  8:228845/149‑1 (MQ=255)
                                                                                                                                        gcACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATcgcgc  <  5:265492/149‑1 (MQ=255)
                                                                                                                                        cgtcgcgCATGCCATCGGCAAGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATcgcgc  >  2:347211/4‑149 (MQ=255)
                                                                                                                                          aCGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCtgctgc                    >  6:138037/1‑129 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CGTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGC  >  NC_002947/4804974‑4805258

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: ‑