Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,805,130 T→A 13.2% H921L (CAT→CTT)  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,1300TA13.2% 113.7 / 7.2 44H921L (CAT→CTT) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base A (3/3);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGC  >  NC_002947/4804985‑4805258
                                                                                                                                                 |                                                                                                                                 
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)
                                                                                                                                                 |                                                                                                                                 
TCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACC‑ACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGC  >  NC_002947/4804985‑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: ‑