Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,024,339 A→G 12.9% intergenic (+153/‑12) xdhD → / → xanQ fused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit/xanthine:H(+) symporter XanQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,024,3390AG12.9% 71.5 / 3.7 31intergenic (+153/‑12)xdhD/xanQfused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit/xanthine:H(+) symporter XanQ
Reads supporting (aligned to +/- strand):  ref base A (9/18);  new base G (2/2);  total (11/20)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

GCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCC  >  NC_000913/3024209‑3024476
                                                                                                                                  |                                                                                                                                         
gcCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGAGAAGGGa                                                                                                                               >  1:573514/1‑143 (MQ=255)
 acTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGGGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGAt                                                                                                                              <  1:118380/142‑1 (MQ=255)
           tGTCTGAACGATGACTAAATGGTGACTGCGACGTTGCCATAAACCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGAtata                                                                                                                    <  2:594103/143‑1 (MQ=255)
             tttGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATaaa                                                                                                                  <  1:244024/143‑1 (MQ=255)
               tGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAAcc                                                                                                                <  2:269898/143‑1 (MQ=255)
                   cGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGc                                                                                                            >  1:671843/1‑143 (MQ=255)
                     aTGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAg                                                                                                          <  2:568848/143‑1 (MQ=255)
                      tGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAgg                                                                                                         >  2:258507/1‑143 (MQ=255)
                           aaaTAGTGACTGCGGAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTg                                                                                                    <  2:184225/143‑1 (MQ=255)
                           aaaTAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTg                                                                                                    <  2:508697/143‑1 (MQ=255)
                            aaTAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGa                                                                                                   <  2:238223/143‑1 (MQ=255)
                             aTAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGAc                                                                                                  <  2:285345/143‑1 (MQ=255)
                             aTAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGAc                                                                                                  <  1:781879/143‑1 (MQ=255)
                                  gACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTGTTGGGTCCCAAAAGCTTTGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTtat                                                                                             <  1:611423/143‑1 (MQ=255)
                                       cGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGAGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTg                                                                                        >  1:663689/1‑143 (MQ=255)
                                           aaTGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACt                                                                                    <  1:678176/143‑1 (MQ=255)
                                               agaTAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATATGAACTggag                                                                                >  2:205458/1‑143 (MQ=255)
                                                gaTAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTggagg                                                                               <  1:297723/143‑1 (MQ=255)
                                                  tAAGCCGTTCTCAAGGCGTCCTGCGCCCCACATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGAt                                                                             <  1:310273/143‑1 (MQ=255)
                                                     gCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGGTCAAATTTACGTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGc                                                                          <  2:642387/143‑1 (MQ=255)
                                                       cGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAGGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGccc                                                                        >  2:74832/1‑143 (MQ=255)
                                                       cGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAGGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGccc                                                                        >  2:486911/1‑143 (MQ=255)
                                                                         cgacccAACTTCGTGGGAATTTATTTTACTAATTCAGAGGATCAAATTTACTTTAAAAGGAGTGCAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGctct                                                      <  2:830466/140‑1 (MQ=255)
                                                                              cAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGTTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTg                                                 <  1:9149/143‑1 (MQ=255)
                                                                              cAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTg                                                 >  1:313030/1‑143 (MQ=255)
                                                                                                 tttACTAATTCAGATGATCAAATTTACTTTAAAGGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGtt                              <  1:74832/143‑1 (MQ=255)
                                                                                                 tttACTAATTCAGATGATCAAATTTACTTTAAAGGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGtt                              <  1:486911/143‑1 (MQ=255)
                                                                                                        aTTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCaaca                       >  2:286537/1‑141 (MQ=255)
                                                                                                                          aCTTTAAAAGGAGTGAAGGGCTGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAAcc     <  1:593787/143‑1 (MQ=255)
                                                                                                                            tttAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAAcccc   >  1:685440/1‑143 (MQ=255)
                                                                                                                             ttAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTCTTGGCACTATTCGTTCCGATGGTAGccccc  >  2:827387/1‑143 (MQ=255)
                                                                                                                                  |                                                                                                                                         
GCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCC  >  NC_000913/3024209‑3024476

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: