Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,003,834 G→C A116P (GCA→CCA)  xdhC → putative xanthine dehydrogenase iron‑sulfur‑binding subunit XdhC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,003,8340GC96.8% 103.0 / ‑3.4 31A116P (GCA→CCA) xdhCputative xanthine dehydrogenase iron‑sulfur‑binding subunit XdhC
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (14/16);  total (14/17)
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.26e-01

ATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGT  >  NC_000913/3003757‑3003908
                                                                             |                                                                          
aTGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCa                                                                          <  1:548258/80‑1 (MQ=255)
      agcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGGTTATGGCTACCACGCCAATGCTgg                                                                   >  2:548283/1‑81 (MQ=255)
      agcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTg                                                                    <  1:548259/80‑1 (MQ=255)
      agcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTgg                                                                   >  2:548299/1‑81 (MQ=255)
         agGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGa                                                                <  2:548261/81‑1 (MQ=255)
         agGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGa                                                                <  2:548260/81‑1 (MQ=255)
                tGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCAcg                                                         >  2:548317/1‑81 (MQ=255)
                   gAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCAcgc                                                        >  2:548310/1‑79 (MQ=255)
                    aaaTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCAcgc                                                        >  1:548311/1‑78 (MQ=255)
                     aaTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCgaga                                                    <  2:548262/81‑1 (MQ=255)
                      aTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGaa                                                   <  2:548263/81‑1 (MQ=255)
                       tCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAg                                                  <  1:548264/81‑1 (MQ=255)
                         cGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGcc                                                <  2:548265/81‑1 (MQ=255)
                            cgcAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCAtt                                             <  2:548266/81‑1 (MQ=255)
                                   cagtgcGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCAt                                       <  2:548267/80‑1 (MQ=255)
                                     gtgcGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTAc                                    <  1:548268/81‑1 (MQ=255)
                                       gcGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACg                                   >  1:548314/1‑80 (MQ=255)
                                        cGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACgg                                  >  2:548326/1‑80 (MQ=255)
                                         gggTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGa                                 >  1:548304/1‑80 (MQ=255)
                                          ggTTTTGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGaaa                               <  1:548269/81‑1 (MQ=255)
                                            ttttGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAAtt                             <  1:548270/81‑1 (MQ=255)
                                            ttttGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAAtt                             >  2:548322/1‑81 (MQ=255)
                                              ttGTACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATtcg                           >  1:548298/1‑81 (MQ=255)
                                                 tACGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATtcgtc                         <  1:548271/80‑1 (MQ=255)
                                                  aCGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTcgc                       >  1:548329/1‑81 (MQ=255)
                                                  aCGCCTGGCCTGATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTcgc                       >  1:548293/1‑81 (MQ=255)
                                                          cctgATTATGGCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGc               <  1:548272/81‑1 (MQ=255)
                                                                   ggCTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCt      >  1:548327/1‑81 (MQ=255)
                                                                     cTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCttt    <  2:548273/81‑1 (MQ=255)
                                                                     cTACCACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCttt    >  2:548302/1‑81 (MQ=255)
                                                                        ccACGCCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGt  <  1:548274/80‑1 (MQ=255)
                                                                             |                                                                          
ATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATCTTTGT  >  NC_000913/3003757‑3003908

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