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.94.7% 72.4 / ‑1.0 19intergenic (‑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 . (13/5);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01
*NC_0009132,173,3640C.94.7% 72.3 / ‑1.0 19intergenic (‑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 . (13/5);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

CAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTT  >  NC_000913/2173134‑2173549
                                                                                                                                                                                                                                     ||                                                                                                                                                                                         
cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGCAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACg                                                                                                                                                                     >  1:65113/1‑251 (MQ=255)
             cccgacgccctcccgatCGCGACCCACGGGGTGATGCTCAAAATGTTGACCCCCGAGATTAAGGTGCGTAACAGATTCCCCCGATCCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCACCGTATTCCCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAACCAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACTCCGTAGCAAGGcc                                                                                                                                                        <  2:415521/235‑1 (MQ=255)
                  aaTAGTTTGCTGCGCGACCCACAGGTGGATGCTCAGAAGGATGACACCCGAGATTAAGGTGCGGAACAGATTCCCAGGATGCACGGCGCCCGCCATCGCCACGAAGAAGCAGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAATGGGGAGGGGGATAAAAAGCAGGCTTGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                                                                                   <  2:348263/251‑1 (MQ=255)
                      gTTTGCGTCGCGAGCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGAGGCACGGCGACCGCCATCGCCCCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGGGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                                                                                                                                                <  2:448915/250‑1 (MQ=255)
                               gcgATCCACAGGGTGATGCCCATAATGAGGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCCCGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                      <  2:457220/251‑1 (MQ=255)
                                cgATCCCCAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                        >  1:415521/1‑232 (MQ=255)
                                      acaGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                                                                               >  2:356942/1‑251 (MQ=255)
                                                         gatgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                            >  2:148250/1‑251 (MQ=255)
                                                                accgagattatcgtcaggAACAGATTTCCCCGATGCACCACCACCGCCATCGCCAAGAAGAAGCCGCTGGTGGAAAGATCGACAAACCGCAGCACCTAATTACCCGGCACACAACCAGCCATTAACATGGTGAGGGGGATACAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGCCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCTACgcgcttg                                                                                                     >  2:193888/17‑245 (MQ=255)
                                                                     gaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCACACGGCAGAACCTGATTACACGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCACAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                  >  1:427002/1‑249 (MQ=255)
                                                                                         ttCCACGATGCATGGCGACCGCCATCGCCACGCAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTTTAAACGACTACGAGCTTGCTTAGCGATGGTCGTTAAAccatcc                                                                            >  2:33310/1‑251 (MQ=255)
                                                                                                                                                        gcagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGc             >  1:150254/1‑251 (MQ=255)
                                                                                                                                                           gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt          >  1:195771/1‑251 (MQ=255)
                                                                                                                                                                   ttACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt  >  2:138592/1‑251 (MQ=255)
                                                                                                                                                                    tACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACg                                                                                                                 >  1:38379/1‑139 (MQ=255)
                                                                                                                                                                                             aaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                    >  2:154538/1‑159 (MQ=255)
                                                                                                                                                                                             aaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                    <  1:154538/159‑1 (MQ=255)
                                                                                                                                                                                                ggTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCTTGCTTAGCGATGGGCGTTAAAcc                                                                                >  2:227604/1‑144 (MQ=255)
                                                                                                                                                                                                ggTGAGGGGGATAAAAATCAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAcc                                                                                <  1:227604/144‑1 (MQ=255)
                                                                                                                                                                                                                                     ||                                                                                                                                                                                         
CAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTT  >  NC_000913/2173134‑2173549

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