Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,001,304 (G)7→9 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.G100.0% 136.5 / NA 38A318G (GCG→GGG) xdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (21/17);  total (21/17)
*NC_0009133,001,2972.G92.3% 128.8 / 3.6 39coding (953/2259 nt)xdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (20/16);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

CCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTT  >  NC_000913/3001165‑3001427
                                                                                                                                     ||                                                                                                                                  
cccGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑gggggg                                                                                                                              <  1:1004356/139‑1 (MQ=255)
                 tttACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTa                                                                                                               >  2:855579/1‑139 (MQ=255)
                          gACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATcc                                                                                                      >  1:238725/1‑139 (MQ=255)
                              ggCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGt                                                                                                  >  1:342332/1‑139 (MQ=255)
                                       gCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTAc                                                                                         >  1:153945/1‑139 (MQ=255)
                                         gTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGc                                                                                       >  1:551299/1‑139 (MQ=255)
                                          tGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                            <  2:320913/133‑1 (MQ=255)
                                          tGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                            >  1:320913/1‑133 (MQ=255)
                                                 cgACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGtt                                                                               >  1:1057354/1‑139 (MQ=255)
                                                      gAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAg                                                                          <  2:272348/139‑1 (MQ=255)
                                                                   ttATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑gggggg                                                                                                                              <  2:551866/72‑1 (MQ=255)
                                                                   ttATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑gggggg                                                                                                                              >  1:551866/1‑72 (MQ=255)
                                                                     ataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAc                                                           <  2:1057354/139‑1 (MQ=255)
                                                                           tGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTCCTGTTCAAAGACCTCCTATCCCAACCTcccc                                                     >  1:931024/1‑139 (MQ=255)
                                                                             gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGt                                                                                                  <  1:293407/92‑1 (MQ=255)
                                                                             gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGt                                                                                                  >  2:293407/1‑92 (MQ=255)
                                                                               tGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCgg                                                 <  1:552143/139‑1 (MQ=255)
                                                                                     gTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAg                                                                          <  2:796629/108‑1 (MQ=255)
                                                                                     gTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAg                                                                          >  1:796629/1‑108 (MQ=255)
                                                                                          acacCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      <  2:153945/139‑1 (MQ=255)
                                                                                           cacCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTg                                           <  1:43699/133‑1 (MQ=255)
                                                                                           cacCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTg                                           >  2:43699/1‑133 (MQ=255)
                                                                                              cGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTg                                  >  1:407874/1‑139 (MQ=255)
                                                                                                gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGt                                >  2:605581/1‑139 (MQ=255)
                                                                                                gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGATAAAGTCGCt                                                                                                                 >  1:305083/1‑57 (MQ=255)
                                                                                                gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGGATAAAGTCGCt                                                                                                                 <  2:305083/57‑1 (MQ=255)
                                                                                                            ctcACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTtgtg                                                                                              >  1:112994/1‑65 (MQ=255)
                                                                                                            ctcACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCt                                                                      >  1:898862/1‑89 (MQ=255)
                                                                                                            ctcACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCt                                                                      <  2:898862/89‑1 (MQ=255)
                                                                                                            ctaaCGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTtgtg                                                                                              <  2:112994/62‑1 (MQ=255)
                                                                                                              cACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTAc                                                                                   <  2:704896/74‑1 (MQ=255)
                                                                                                              cACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCt                                                                      <  2:333247/87‑1 (MQ=255)
                                                                                                              cACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCt                                                                      >  1:333247/1‑87 (MQ=255)
                                                                                                              cACGGGCACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACACTTAc                                                                                   >  1:704896/1‑74 (MQ=255)
                                                                                                                cGGGCCCTCCATCGCTTCTGCGGGGGGGGGAATAAATTCCCTTACCTTTATCCTCGTTGTTCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAg                <  2:931024/139‑1 (MQ=255)
                                                                                                                    cACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTtgtg                                                                                              >  2:618511/1‑57 (MQ=255)
                                                                                                                    cACTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTtgtg                                                                                              <  1:618511/57‑1 (MQ=255)
                                                                                                                      cTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCt                                                                                                                 <  2:79627/36‑1 (MQ=255)
                                                                                                                      cTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCt                                                                                                                 >  1:79627/1‑36 (MQ=255)
                                                                                                                      cTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGt                                                                                                  <  2:423255/51‑1 (MQ=255)
                                                                                                                      cTCCATCGCTTCTGCGGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGt                                                                                                  >  1:423255/1‑51 (MQ=255)
                                                                                                                             gCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt  <  1:1039186/139‑1 (MQ=255)
                                                                                                                                     ||                                                                                                                                  
CCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTT  >  NC_000913/3001165‑3001427

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