Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑490/+918) gatD ← / ← gatB galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3630C.88.9% 25.7 / ‑3.5 9intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (5/3);  total (5/4)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3640C.88.9% 25.7 / ‑3.5 9intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (5/3);  total (5/4)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173113‑2173606
                                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                                                  
tAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGCCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGCTACAAATCAGGCTTGCCGATACCAccc                                                                                                                                                                                                                                                     >  2:268067/1‑251 (MQ=255)
                                           gTTTGCGTCGAGCTCCACAAGGTTTTGCTCATAAGGATGACACCCGAGATTAAGGTGCGGAAAAGATTTCCACGATGCACGGCGACCGCCATCGCAACGAAGAAGACGGTGTAGGCAAGATCGCCAAACGGCCGCACCTGAGTACCCGGCCCAAAAACAGCAATTAAAATGGGGCGATGGAAAAAAATAAGGCTTGACGAGTCCACC‑‑GCAGTATGTCCAAGCAACAAGTCCGGAACAAGTCCAATCaagc                                                                                                                                                                                                         >  2:221907/1‑247 (MQ=255)
                                                   ggcgttCCACAGGGTGATCCTCATAATGATGCTCCCCGAGATTAAGGTGGGGAACAGATTCCCACGATGCACGGCGCCCGCCCTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                                                                                                                                 <  1:206310/245‑1 (MQ=255)
                                                     cgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCACCACCGGATCAAGGCCAATCAGGAACGCCTGAccc                                                                                                                                                                                              >  2:202755/1‑250 (MQ=255)
                                                                gTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                                                                                                                                   >  1:100766/1‑251 (MQ=255)
                                                                                                  tGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGc                                                                                                                                                 >  1:45131/1‑251 (MQ=255)
                                                                                                          aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGCCCAGATAGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAcc                                                                                                                                         >  1:196145/1‑251 (MQ=255)
                                                                                                               cccACGATCCCCGGCGACCGCAATCGCCACGAAGAAGCCGATGGTGGCAAGACCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatcca                                                                                                                                    <  1:140987/250‑1 (MQ=255)
                                                                                                                                                                              cagcaCCTGATTCCCCGACAAACAAACAGAATTAAAAATGGTGAGTGGGAAAAAAATAGGCCTTGCCTATCCCCCCGCCGTATGTCCCGGCAGCAACGCCGGATCAAGGCCAATCAGTAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCa                                                                     <  1:198231/251‑1 (MQ=255)
                                                                                                                                                                                      aattccccGGCACAAAAAAACCAATAAAAATGGGGGGTGGAATAAAAACCGGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCt                                                             <  1:202755/246‑1 (MQ=255)
                                                                                                                                                                                                                                                     gataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATAATGGGTTTAATAACCCGTGGCCTTAGCAGCATCACTGCCGCCGTTTTTAACGCCAGCTGCAATACACATTTGACATCATAACCCGCGAGGATGACGATAATCAGCCCAATCACAAAACCCACGGTGAc  >  2:207808/8‑251 (MQ=255)
                                                                                                                                                                                                                                                            ||                                                                                                                                                                                                                                                  
TAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAC  >  NC_000913/2173113‑2173606

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