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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,001,2971.G100.0% 150.5 / NA 43A318G (GCG→GGG) xdhAputative xanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (22/21);  total (22/21)
*NC_0009133,001,2972.G95.3% 143.5 / 0.2 43coding (953/2259 nt)xdhAputative xanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (21/20);  total (22/21)
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.45e-01

ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTT  >  NC_000913/3001170‑3001427
                                                                                                                                ||                                                                                                                                   
aCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑‑gggggggggaat                                                                                                                        >  2:857367/1‑135 (MQ=255)
 cccGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑‑gggggggggaata                                                                                                                       <  2:868837/140‑7 (MQ=255)
  ccGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑‑gggggggggaataa                                                                                                                      <  1:403043/140‑8 (MQ=255)
     ccACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGt                                                                                                                     >  1:323885/1‑140 (MQ=255)
      cACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTc                                                                                                                    >  2:292385/1‑140 (MQ=255)
      cACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTc                                                                                                                    >  1:936937/1‑140 (MQ=255)
        cGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGc                                                                                                                  <  1:365593/140‑1 (MQ=255)
        cGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGc                                                                                                                  <  2:323885/140‑1 (MQ=255)
               aCCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCtt                                                                                                           >  2:685489/1‑140 (MQ=255)
                     gACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCt                                                                                                     <  1:685489/140‑1 (MQ=255)
                         ggCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGtt                                                                                                 <  1:292385/140‑1 (MQ=255)
                              aTGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                            <  1:941267/140‑1 (MQ=255)
                              aTGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                            <  1:213041/140‑1 (MQ=255)
                                   cGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGc                                                                                       >  1:621880/1‑140 (MQ=255)
                                    gTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCt                                                                                      >  1:782187/1‑140 (MQ=255)
                                    gTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTATCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCt                                                                                      <  1:777704/140‑1 (MQ=255)
                                           gcgACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGtt                                                                               >  1:827885/1‑140 (MQ=255)
                                            cgACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGtt                                                                               >  1:496741/1‑140 (MQ=255)
                                            cgACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa                                                                                    >  1:623194/1‑134 (MQ=255)
                                            cgACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa                                                                                    <  2:623194/134‑1 (MQ=255)
                                             gACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCa                                                                             <  2:782187/140‑1 (MQ=255)
                                                ggAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAg                                                                          <  2:415753/140‑1 (MQ=255)
                                                         aaaGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCtat                                                                 >  1:616600/1‑140 (MQ=255)
                                                            ggTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCG‑‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCa                                                             >  2:867212/1‑140 (MQ=255)
                                                                 taGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCG‑‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTc                                                        <  1:867212/140‑1 (MQ=255)
                                                                      tGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCt                                                    >  2:360565/1‑140 (MQ=255)
                                                                       ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAcc                                                          >  2:460789/1‑133 (MQ=255)
                                                                       ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAcc                                                          <  1:460789/133‑1 (MQ=255)
                                                                        gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCg                                                  >  1:318112/1‑140 (MQ=255)
                                                                           gTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCt                                               >  2:609634/1‑140 (MQ=255)
                                                                             tCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTgg                                             >  2:836481/1‑140 (MQ=255)
                                                                              cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGt                                            >  1:74803/1‑140 (MQ=255)
                                                                              cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGt                                            >  2:742744/1‑140 (MQ=255)
                                                                                 tCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCg                                         <  2:647297/140‑1 (MQ=255)
                                                                                 tCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCg                                         <  2:827885/140‑1 (MQ=255)
                                                                                    aaCACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      <  2:936937/140‑1 (MQ=255)
                                                                                           gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGtt                               >  1:541618/1‑140 (MQ=255)
                                                                                                       ctcACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCcac                   >  1:863519/1‑140 (MQ=255)
                                                                                                       ctcACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCcac                   <  1:360565/140‑1 (MQ=255)
                                                                                                        tcACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCcaca                  <  2:541618/140‑1 (MQ=255)
                                                                                                         cACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACaa                 >  2:19664/1‑140 (MQ=255)
                                                                                                                 cTCCATCGCTTCTGCGGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTAt          <  2:496741/140‑1 (MQ=255)
                                                                                                                      tCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCg    <  2:318112/140‑1 (MQ=255)
                                                                                                                        gCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt  >  2:984794/1‑140 (MQ=255)
                                                                                                                        gCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt  <  1:836481/140‑1 (MQ=255)
                                                                                                                        gCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt  <  1:742744/140‑1 (MQ=255)
                                                                                                                                ||                                                                                                                                   
ACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTT  >  NC_000913/3001170‑3001427

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