Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,001,304 (G)7→8 coding (960/2259 nt) xdhA → xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,001,2971.G97.8% 162.6 / ‑3.2 45coding (953/2259 nt)xdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (20/24);  total (20/25)
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.65e-01

CGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTT  >  NC_000913/3001160‑3001427
                                                                                                                                          |                                                                                                                                   
cGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑g                                                                                                                                   >  2:749886/1‑139 (MQ=255)
cGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑g                                                                                                                                   <  2:449778/139‑1 (MQ=255)
cGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTACATCGCTTCTGC‑‑g                                                                                                                                   >  1:198305/1‑139 (MQ=255)
  cAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGTAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑ggg                                                                                                                                 >  1:42289/1‑139 (MQ=255)
  cAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑ggg                                                                                                                                 >  1:817208/1‑139 (MQ=255)
    aCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑ggggg                                                                                                                               <  2:327679/139‑1 (MQ=255)
           cccGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑ggggggggaata                                                                                                                        >  1:1192862/1‑134 (MQ=255)
           cccGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑ggggggggaata                                                                                                                        <  2:1192862/139‑6 (MQ=255)
           cccGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑ggggggggaata                                                                                                                        <  2:66174/139‑6 (MQ=255)
                cACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTTAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTc                                                                                                                    >  2:844293/1‑139 (MQ=255)
                   gCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATGGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                 <  2:882264/139‑1 (MQ=255)
                    cTTTTGCCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCTCTTCTGCGGGGGGGGGAATAAAGTCGCt                                                                                                                 <  1:749886/139‑1 (MQ=255)
                              tGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATcc                                                                                                      >  1:1055683/1‑139 (MQ=255)
                               gACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCt                                                                                                     >  1:44949/1‑139 (MQ=255)
                               gACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCt                                                                                                     >  1:577865/1‑139 (MQ=255)
                                 cGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCg                                                                                                   <  2:577865/139‑1 (MQ=255)
                                     cAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTtgtg                                                                                              <  2:220077/139‑1 (MQ=255)
                                        aTGGGCGTGAACCGCGACGGAACAGTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                            <  1:676152/139‑1 (MQ=255)
                                               tGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATTTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCtt                                                                                     <  1:622933/139‑1 (MQ=255)
                                                           gAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGa                                                                         <  2:979946/139‑1 (MQ=255)
                                                              cATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCt                                                                      <  1:337644/139‑1 (MQ=255)
                                                                 tGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCt                                                                   <  1:497449/139‑1 (MQ=255)
                                                                        ttATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTa                                                                                                               <  2:300741/88‑1 (MQ=255)
                                                                        ttATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTa                                                                                                               >  1:300741/1‑88 (MQ=255)
                                                                                 ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTc                                                   >  1:419350/1‑139 (MQ=255)
                                                                                  gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                 >  2:1122502/1‑76 (MQ=255)
                                                                                  gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                 <  1:1122502/76‑1 (MQ=255)
                                                                                  gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCaaa                                                                           >  1:173364/1‑114 (MQ=255)
                                                                                  gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCaaa                                                                           <  2:173364/114‑1 (MQ=255)
                                                                                      ttCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTg                                              <  1:506306/139‑1 (MQ=255)
                                                                                        cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGt                                            >  1:299486/1‑139 (MQ=255)
                                                                                        cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑gg                                                                                                                                  <  1:791191/52‑1 (MQ=255)
                                                                                        cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑gg                                                                                                                                  >  2:791191/1‑52 (MQ=255)
                                                                                             tAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGAt                                       >  1:671589/1‑139 (MQ=255)
                                                                                              aaCACCGGCGCTTACGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      <  1:844293/139‑1 (MQ=255)
                                                                                                     gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGCCTGGTGCGATGCGTGGt                                >  1:77925/1‑139 (MQ=255)
                                                                                                     gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTAc                                                                                                              >  2:399078/1‑60 (MQ=255)
                                                                                                     gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTAc                                                                                                              <  1:399078/60‑1 (MQ=255)
                                                                                                     gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGa                                                                         >  1:723582/1‑97 (MQ=255)
                                                                                                     gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGa                                                                         <  2:723582/97‑1 (MQ=255)
                                                                                                                 ctcACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa                                                                                    >  1:744400/1‑74 (MQ=255)
                                                                                                                 ctcACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa                                                                                    <  2:744400/74‑1 (MQ=255)
                                                                                                                   cACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAg                                                                                                                      <  1:781644/38‑1 (MQ=255)
                                                                                                                   cACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAg                                                                                                                      >  2:781644/1‑38 (MQ=255)
                                                                                                                           cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCa                                                                             <  2:613991/71‑1 (MQ=255)
                                                                                                                           cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCa                                                                             >  1:613991/1‑71 (MQ=255)
                                                                                                                           cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTg                                                                     <  1:90312/79‑1 (MQ=255)
                                                                                                                           cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTg                                                                     >  2:90312/1‑79 (MQ=255)
                                                                                                                           cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      >  1:450072/1‑110 (MQ=255)
                                                                                                                           cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      <  2:450072/110‑1 (MQ=255)
                                                                                                                              cATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGc      >  1:200245/1‑139 (MQ=255)
                                                                                                                                tCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCca                    <  1:369508/122‑1 (MQ=255)
                                                                                                                                tCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCca                    >  2:369508/1‑122 (MQ=255)
                                                                                                                                 cGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGt   <  2:419350/139‑1 (MQ=255)
                                                                                                                                  gCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTTCGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt  <  1:659573/139‑1 (MQ=255)
                                                                                                                                  gCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt  <  2:1173741/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                   
CGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTT  >  NC_000913/3001160‑3001427

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