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.87.0% 76.1 / 6.0 23intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (10/10);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01
*NC_0009132,173,3640C.87.0% 77.6 / 4.5 23intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (10/10);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAC  >  NC_000913/2173130‑2173596
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                        
ccgcgAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGGGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATCCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCCCACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATg                                                                                                                                                                                                                                  <  2:13204/241‑1 (MQ=255)
ccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATg                                                                                                                                                                                                                                  >  1:13204/1‑241 (MQ=255)
       ttGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAAGTGAGGAACCGATTTCCACGATGCACGACCACAGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCg                                                                                                                                                                                                                 >  1:334309/1‑251 (MQ=255)
                            tggCGTCGCGAGCCACAGGGGGATGCTCAAAAGGATGACCCCCGAAATTAAGTGGCGGAACAGGTTTCCACGATGCACGGCGACCTCCATCGCCAGGAAGAAGCCGAGGGTGGAAAGATCTCCAAACGGCAGCACCGGATTACCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAACCGGGCTTGCCGATACCACC‑‑GCCGGATGTCCCAGCAGCAACGCCGGAGCAAGGCCAATCAGGAAc                                                                                                                                                                                             <  2:334309/248‑1 (MQ=255)
                                gTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                        >  2:277275/1‑251 (MQ=255)
                                  cgcgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                                                                                                                      >  1:243436/1‑251 (MQ=255)
                                     gATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                                                                                                                   >  2:80472/1‑251 (MQ=255)
                                     gATCCACAGGGAGATGCTCATAATGCTGACACCCGAGCTTAAGCTCCAGCACAGATTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGAAGCCGCTGGCGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACCCAACCAG‑CAATTAACA‑TGGTGAGTGCGATAAAACTCAGGCTTGCCGCTACCACC‑‑GCCGTATGTCCCAGCAGCCACGCCGGATCAAGGCCAATCCGGACCTCCTGAccgc                                                                                                                                                                                   >  2:345926/1‑251 (MQ=255)
                                            aGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                                                                                                                            >  1:2012/1‑251 (MQ=255)
                                                                gatacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                          >  1:183784/4‑249 (MQ=255)
                                                                   ccccGGGATTAAGGGGGGGAACAGTTTCCCCCGAGGCACGGCGCCCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                                                                                                                     <  2:2012/250‑1 (MQ=255)
                                                                      cGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCCCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                                                                                                                                  <  1:194132/251‑1 (MQ=255)
                                                                                             ttCCACGATGCACGGCGACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGTTCGCAAAACTGCAGCACCTGATTACCCGGCACACAAACAG‑CAATAAAAA‑TGGGGAGTGGTATAAAAAACAGACTTGCCGATGCGACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCGAGGCCAATCATGAACTCCTGACCGCCGAACTACGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                                                                                                                            <  1:178306/250‑1 (MQ=255)
                                                                                              ccccggaggCACGGCGCCCGCCATCGCCAAGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCCCACAAACAG‑CAATTAAAA‑GGGTGAGGGGGATAAAAATCAGGCTTGCCGAAACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                           <  2:175115/242‑1 (MQ=255)
                                                                                                        ccGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                <  1:277275/250‑1 (MQ=255)
                                                                                                                                      gtggtggCAAGATCGCCAAACGGCAGCACCTGATAACCCGGCACACAAA‑AGTCAATTAAAA‑GGGTGAGTGTGATAAAATCCAGGCTTGCAGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGTCCAATCAGGATCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgc                                                                                  <  2:329587/250‑1 (MQ=255)
                                                                                                                                       ggggggCAAGATCGCCAAACGGCAGCACCGTATTCCCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgcc                                                                                 <  1:55504/247‑1 (MQ=255)
                                                                                                                                                             cagccctggATTACCCGGCACAAAAACAG‑AAATTAAAA‑TGGGGAGGGGAATAAAAATCAGGCTGGCCGAT‑CC‑CCCCGCCGTATTGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCa                                                           <  1:25394/243‑1 (MQ=255)
                                                                                                                                                                 cccTGATTCCCCGGCACACAAAAAA‑CAATTAAAA‑TGGTGAGGGGGAAAAAAATCAGGGTTCCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTTACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                       <  1:27990/250‑1 (MQ=255)
                                                                                                                                                                  ccTGATTCCCCGGCACACAAACAC‑CAATTAAAA‑GGGTGAATGGCATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                                 <  2:189064/160‑1 (MQ=255)
                                                                                                                                                                  ccTGATTACCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                                 >  1:189064/1‑160 (MQ=255)
                                                                                                                                                                                ccccaaaaacacAATTAAAAATGGTGAGTGGAATAAAAACCAGTCTTCCTG‑‑ACCCCCCCGCCGTATGCCCCAGCAGCAACGCCGGAGCAGGGCAAATCAGAAACTCTGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCTTTAAACCATCCAGGATGGGTTTAATCGCCCGTGGCATTAGCAGCATCACTTCCGCCGTTTTTTCCGCCAGCTGCAATACACCTTTGTCATCg                                          <  2:208605/240‑1 (MQ=255)
                                                                                                                                                                                             aTTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGAAACCACG‑‑GAAGTATGTCCCAGCAGCAACGCCGGTTCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGATAAACGATCCATGAGGGGTTTAATCTCCCGTGGCATTAGCCGCAACACTGCCGCCGTTTTTACCGCCAGCTGGAATACACCTTTGACATGGTAACCcgcgtgaat                            >  2:19030/1‑246 (MQ=255)
                                                                                                                                                                                                                          aagccgtgcagatccacccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAAGTCCGGACCGCCGAACTTCGCCTGTAAAGAACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTGTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAc  <  1:41692/235‑1 (MQ=255)
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                        
CCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAA‑TGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAC  >  NC_000913/2173130‑2173596

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