Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3640C.75.0% 51.1 / 24.1 28intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/7);  new base . (13/8);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 6.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCT  >  NC_000913/2173122‑2173611
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
aGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGAAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTc                                                                                                                                                                                                                                               <  1:204218/251‑1 (MQ=255)
     tAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                                          >  1:11191/1‑251 (MQ=255)
         cgccAGATGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGGGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGATTGCCGATACCACC‑‑GCAGTATGTCCCagcagc                                                                                                                                                                                                                                       >  2:336797/1‑250 (MQ=255)
                                                  acaGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                                                                                                                                             >  1:101247/1‑251 (MQ=255)
                                                        ggatcttCATAATGAGGACCCCCAAGATTAGGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCACCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAACCGGCAGCACCTG‑ATTACCCGGCACCCAAACAGCAATTAAAATGGGGAGGGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                                                                                                                                       <  1:140745/245‑1 (MQ=255)
                                                           ggCTCATAATGATGCCCCCCGAGATTAAGGTGCGAAACAGTTTCCCACGATGCCCGGCGCCCGCCATCGCCCCGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTG‑ATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                                                                                                                                    <  1:292728/250‑1 (MQ=255)
                                                              tCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGCAGAAGCAGATGGTGGCAAGAACGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCACTTAAAATGGGGAGCGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAACAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTaa                                                                                                                                                                                  >  2:89603/1‑250 (MQ=255)
                                                                 aatggaagaCACCCGATATTAGGGTGAGGACCAGTTTCCCACGATGCACGGCGCCCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCACCAC‑TGTATTACCCGGCACACAAACAGCAATAAAAATGGTGATTGGGATAAAACTCAGGCTGGCCGA‑ACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                                                                                                                                              <  1:376730/244‑1 (MQ=255)
                                                                  aaTGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAATTCGCCTGTAAACGAc                                                                                                                                                                             >  2:179669/1‑251 (MQ=255)
                                                                   atgatgACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                            >  1:197286/1‑251 (MQ=255)
                                                                                         tGCGGAAAAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCCATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCGGCAGCAACGCCGGATAAGGGCCAATCAGGAACTCCTGACCGACGACCTTCGCCTGTAAACGCCTACGCGCTTGCTTAGCGCTgggg                                                                                                                                                      >  2:33853/1‑250 (MQ=255)
                                                                                              ataaaatgtcacgtagccACGGCCACCGCCACCGCCACGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCACCCCCTG‑ATTCCCCGCCCCACAAACAGCAATAAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTaa                                                                                                                                                 <  1:74274/234‑1 (MQ=255)
                                                                                                             tGCACGGCGACCGCCATCGCCACGAAGAAGCAGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGAAACCACC‑‑GCCGTATGTCCCAGCAGCAACGCAGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGAAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                                                                                  >  2:305010/1‑251 (MQ=255)
                                                                                                               cACGGCGCCCGCCAGCCCCCCGAAGAAGCCTTGGGGGGCAAGATCGCCAAAGGGCAGCCCCTT‑ATTACCCGGCACCAAAACAGCATTAAAAATGGTGAGTGGGATAAAAATCAGTCTGGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACTCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt                                                                                                                                <  1:381691/251‑1 (MQ=255)
                                                                                                                      cccgccaccgccaGGAAGAAGCAGATGGGGTCAAGATCGCCAATCGGCAGAACCTG‑ATTCCCCGTAACACAAACAGCAATGAAAAGGTTGGGTGGGATAAATATCAGGCTTGCTGAT‑CCACCC‑GCCGTATGTCCCATCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCTCCTAACTTCGCCTGTAAACGTCTACGAGCTTGCTTATCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc                                                                                                                         <  1:248843/250‑1 (MQ=255)
                                                                                                                                             gATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTg                                                                                                  >  1:291157/1‑251 (MQ=255)
                                                                                                                                                gggggCAAGACCCCCAAACGGCAGCCCCTG‑ATTCCCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAGACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccg                                                                                               <  1:143701/248‑1 (MQ=255)
                                                                                                                                                  gggaaaGATCGCCAACCGGCAGCACCGG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                             <  1:261006/247‑1 (MQ=255)
                                                                                                                                                  gggCAAGATCGCCAAACGGCAGCACCTG‑ATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                             <  1:336797/250‑1 (MQ=255)
                                                                                                                                                    gcaatacccccaaccggaaccacttG‑ATTACCCGTCACACAACCAGCAATTAAAAGGGTGAGAGGAAAAAAAATCAGGTTTGCCTAT‑CC‑CCCCGCCGTATGTCCCAGCAGCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAACCCATCCATGATGGGTTTAATCACCCGTGGCATTATCAGCATCACTGCCGCCGt                                                                                           <  1:41348/227‑1 (MQ=255)
                                                                                                                                                        gATCGCCAAACGGCAGCCCCTG‑ATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGAATAAAAACCAGGCTTGCCGAGCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                       <  1:305010/251‑1 (MQ=255)
                                                                                                                                                                       agcccgagtgtccccGGCACACAAACAGCAATTAAAATGGGGAGTGGAAAAAAAATCAGCCTGGCCGATAC‑‑CTCCGCCGTATGGCCCAGCAGCAACTCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCAGCCATGATGGGTTAAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCaa                                                                         <  1:354989/239‑1 (MQ=255)
                                                                                                                                                                                     cgccacacaAACAGTAATAAAAATGGTGAGTGGGATAAAAACCAGGCTGCCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAt                                                           <  1:266884/248‑1 (MQ=255)
                                                                                                                                                                                       ccacacaAAAAGCAATAAAAATGGTGAGTGGAAAAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                                         <  1:383297/250‑1 (MQ=255)
                                                                                                                                                                                                 aGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcg                                                >  1:81447/1‑250 (MQ=255)
                                                                                                                                                                                                   cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAgg                                             >  2:291792/1‑251 (MQ=255)
                                                                                                                                                                                                            tGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCTTTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                     >  2:296289/1‑250 (MQ=255)
                                                                                                                                                                                                                           aaaacccagcttgcccatcccccccGCCGTATGTCCCACCAGCACCGCCGGACCAAGGCCAATCAGGATCTCCTGACCTCCGATCTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGTCGTTAAACCATCCATGATGGGTTTAATCTCCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCTTAACCCGCGAGGATGCCTATAATCAGCCCCATca                       <  1:37810/230‑1 (MQ=255)
                                                                                                                                                                                                                                    ttGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCAGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGAAATACACCTTTGACATCGTACACCGCGAGGATGCCGATACTCAGCACCATCAaaaatcacacg            >  2:350428/1‑244 (MQ=255)
                                                                                                                                                                                                                                           gatacccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc       <  1:341151/246‑1 (MQ=255)
                                                                                                                                                                                                                                              accaccgccgTATGTCCCAGCAGCAACGCCGGACCAAGGCCAATCAGGAACTCCTGACAGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGCGTTTAATCACACGAGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCATTACACCTTTGACATCGTAACCCGACAAGATGCCGACAATCAGCCCCATCACAACAACCACGGTGACAgg    >  2:364735/5‑251 (MQ=255)
                                                                                                                                                                                                                                                caccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGcc  >  1:103884/3‑250 (MQ=255)
                                                                                                                                                                                                                                                  ccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATg                                                                                                                                                         >  1:325446/1‑98 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       
AGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTG‑ATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCT  >  NC_000913/2173122‑2173611

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