Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,017,648 A→G 24.7% intergenic (+103/‑62) xdhD → / → xanQ putative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding/xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,017,6480AG24.7% 25.9 / 3.1 16intergenic (+103/‑62)xdhD/xanQputative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding/xanthine permease
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base G (2/2);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAA  >  CP009273/3017439‑3017878
                                                                                                                                                                                                                 |                                                                                                                                                                                                                                      
gcgcAGCTCCTGCGTATGCTACCGCAATTCACGAAGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGTAGAAAATACTAACTGCGCTGGAAAAAATATAAAGAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGt                                                                                                                                                                       <  1:96836/275‑1 (MQ=255)
                                                                      acCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCACATTTACt                                                                                                                                                                                            >  1:543/1‑184 (MQ=255)
                                                                        cGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAGATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAg                                                                                                                                                                                     >  2:97900/1‑189 (MQ=255)
                                                                            gaAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAGATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGt                                                                                                                                                                                 <  1:97901/189‑1 (MQ=255)
                                                                                               gAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGa                                                                                                                                                                               >  2:203636/1‑172 (MQ=255)
                                                                                                   aaaTATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAggg                                                                                                                                                                           <  1:203645/172‑1 (MQ=255)
                                                                                                                                                tGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAatttatt                                                                                                                                                                                                                       >  2:296344/1‑83 (MQ=255)
                                                                                                                                                    cGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTAc                                                                                                                                                                                                                   <  1:296355/83‑1 (MQ=255)
                                                                                                                                                                    aCTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAGATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGAt                                                                                                                        >  1:196845/1‑158 (MQ=255)
                                                                                                                                                                    aCTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGa                                                                                                                                                                               >  1:67610/1‑103 (MQ=255)
                                                                                                                                                                    aCTGCGAAAATGAGATAAGCCGTTCTCAAGGACTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTaa                                                                                                                                                                                                                >  1:424520/1‑70 (MQ=255)
                                                                                                                                                                       gCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAaacaa  >  1:198159/1‑273 (MQ=255)
                                                                                                                                                                        cGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAGATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGcc                                                                                                                    <  2:196838/158‑1 (MQ=255)
                                                                                                                                                                        cGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAggg                                                                                                                                                                           <  2:67610/103‑1 (MQ=255)
                                                                                                                                                                                                    gTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAggg                                                                                                                                                                           >  2:115268/1‑75 (MQ=255)
                                                                                                                                                                                                        tGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGt                                                                                                                                                                       <  1:115269/75‑1 (MQ=255)
                                                                                                                                                                                                                 |                                                                                                                                                                                                                                      
GCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTCCTGCGCCCCAAATTCGTGGGAATTTATTTTACTAATTCAGATGATCAAATTTACTTTAAAAGGAGTGAAGGGATGTCTGATATAAACCATGCAGGTTCTGACCTTATATTTGAACTGGAGGATCGCCCTCCCTTTCATCAGGCTCTCGTTGGTGCCATTACCCATCTGTTGGCAATTTTCGTTCCGATGGTAACCCCCGCGTTAATCGTGGGTGCGGCCTTACAGCTTTCCGCTGAAACAA  >  CP009273/3017439‑3017878

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