Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 543,263 2 bp→GT coding (783‑784/786 nt) allE ← (S)‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,2630AG100.0% 88.0 / NA 29*262Q (TAA→CAA) allE(S)‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/13);  total (16/13)
*NC_000913543,2640AT93.1% 81.1 / NA 29I261I (ATT→ATAallE(S)‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (16/11);  minor base G (0/2);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGT  >  NC_000913/543033‑543486
                                                                                                                                                                                                                                       ||                                                                                                                                                                                                                              
aCCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTa                                                                                                                                                                                                              <  1:148069/250‑1 (MQ=255)
aCCCTTACCTTTAAAGCGGATGCGATTTATATCACATACGAGTCCATTACTCATGCCGGAAGAGGCATGAGTACCATATCCTTCCTGAAAATCGAGCAAATTCTATATATTGCAGAGCTCATGTAGGCCTGATACGCGAAGCGCATCAGGAAATGTTACAAAAAAAGCCACGGTATAAAACGTGGCAAAATCCAACATAGCTAAAAATAATCAAGAGAGTGGAATGACTTGTATCTCTACGTCGCgggat                                                                                                                                                                                                              >  2:148069/1‑247 (MQ=255)
                      cGATTTATATCGCATAAGAGTGCAGTACTCATTCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGACCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCAACGGAATAAACAGTGGCAAACTCCAACATAGCTAAAAATAATCAGGCGAGTGGAATGAATTGTATCACTACGTCGCGGTTACAATCTTTCGAGTAAAAATAGCt                                                                                                                                                                                       >  2:301410/1‑251 (MQ=255)
                         tttATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCTCATCAGGCAATGTTACAAAAAAAGCCACGGAATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAATGGTATGAATTGTATCTCTACGACGCGGTTACAATCTTTAGAGTAAATATATCTGa                                                                                                                                                                                     >  2:284213/1‑250 (MQ=255)
                                     aaGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAACGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACg                                                                                                                                                                        >  2:305390/1‑251 (MQ=255)
                                                aCTCAGGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCAt                                                                                                                                                             <  1:259035/251‑1 (MQ=255)
                                                 cTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCAt                                                                                                                                                             >  1:88686/1‑250 (MQ=255)
                                                    cagCCGGACGCGGCAGACGTACCAGAACCTGCCTGAAAATCGCTCAAAATCTATATATTGCAGAGATCATGTAGGCCTTATAAGCGAAGCTCATCAGGCAATGTTACAAAAAAATCCACGGTATAAACCGTGGCAAAATCCATCATAGCTAAAAATAATCAGTCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCCTCTCCACGACCTACACCaacccc                                                                                                                                                         <  2:228212/249‑6 (MQ=255)
                                                        cGGATGCGGCCTGAGTACCATAACCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGAGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCAAACATAGCTAAAAATAATCAGGCAAGTGGTATGACTTGTATCTCTAAGTCGCGGTTACCATCTTTCGAGTAAATATAGCTGAACGCTTCCCCCCGCACAACACCATAACCAGcc                                                                                                                                                     >  2:88048/1‑251 (MQ=255)
                                                                                                gcAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCtttt                                                                                                              <  1:267914/251‑1 (MQ=255)
                                                                                                             atTGCAGAGCTAATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAACGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGAGGAATGACTTGTATCTCTCCGTCGCGGTTACAATCTTTCGAGTAAAAATAGCTGAACGCTTCACCACGCCCTACACCATAACAAGCCTGTAAAGAATACGCGCCCATAAAGATGTAATCGCCTTTTTTCAACGGGATcc                                                                                                 >  2:141563/1‑251 (MQ=255)
                                                                                                             atTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAACAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATcc                                                                                                 >  1:227219/1‑251 (MQ=255)
                                                                                                                                   gATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCCCCTTTTTTCCCCGGGATCACGTTATTGACGAGGTTATAAAc                                                                           >  2:7488/1‑251 (MQ=255)
                                                                                                                                            agcgCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGTCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATACCCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGAcc                                                                  <  2:88686/251‑1 (MQ=255)
                                                                                                                                                cATCAGGCAATGTTACAAAAAAAGACACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGCGGTTATAAACCCCCTGACCGGa                                                               >  2:177492/1‑250 (MQ=255)
                                                                                                                                                      gCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAata                                                        <  1:7488/251‑1 (MQ=255)
                                                                                                                                                              aCAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCAcc                                                <  1:305390/251‑1 (MQ=255)
                                                                                                                                                                      aGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCt                                        <  1:177492/251‑1 (MQ=255)
                                                                                                                                                                                 aaaCCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTCTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTg                                       >  1:43945/1‑241 (MQ=255)
                                                                                                                                                                                 aaaCCGGGGCAAAAACCAACATAGCTAAAAATAATCAGGCGAGTGGAATGACTTGTTTCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATTTAATCGCCTTTCTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTg                                       <  2:43945/241‑1 (MQ=255)
                                                                                                                                                                                         gCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGCACGCTTCACCACGACCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGatat                      >  1:61276/1‑250 (MQ=255)
                                                                                                                                                                                             aaTCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAg                                                     <  2:111531/215‑1 (MQ=255)
                                                                                                                                                                                             aaTCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAg                                                     >  1:111531/1‑215 (MQ=255)
                                                                                                                                                                                              aTCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTAGAGTAAATATAGCTGAACGCTTCACCCCGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCCTAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCgg                >  2:294951/1‑250 (MQ=255)
                                                                                                                                                                                                  aaCATAGCTAAAAATAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGtggct            <  1:141563/251‑1 (MQ=255)
                                                                                                                                                                                                           aaaaaTAATCAGGCGAGTGGTATGACTTGGATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTGCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCt     <  2:227219/249‑1 (MQ=255)
                                                                                                                                                                                                            aaaaTAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGt  >  1:237661/1‑251 (MQ=255)
                                                                                                                                                                                                             aaaTAATCAGGCGAGTGGTATGTTGTGGATCTCTACGTCGCGGTTACAATCTTGCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGAGATAACCGtggctg           <  2:19090/241‑1 (MQ=255)
                                                                                                                                                                                                             aaaTAATCAGGCGAGTGGTATGACTTGTATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGtggctg           >  1:19090/1‑241 (MQ=255)
                                                                                                                                                                                                                                       ||                                                                                                                                                                                                                              
ACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATAGCTAAAAATAATCAGGCGAGTGGTATGACTTAAATCTCTACGTCGCGGTTACAATCTTTCGAGTAAATATAGCTGAACGCTTCACCACGCCCTACACCATAACCAGCCTGTAAAGAATAAGCGCCCATAAAGATGTAATCGCCTTTTTTCACCGGGATCCAGTTATTGTCGAGGTTATAAACCCCCTGACCGGAAAGAATATAGGCACCGTGTTCCTGAACGTGTGTTTCGATATAACCGTGGCTGGCACCTGGT  >  NC_000913/543033‑543486

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