Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3640C.77.4% 69.1 / 23.2 31intergenic (‑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 . (20/4);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 1.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGC  >  NC_000913/2173132‑2173610
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
gccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAgcagca                                                                                                                                                                                                                                     >  1:122162/1‑250 (MQ=255)
       gggTGTGCAAACCAATAGTTTGAGTCGCGATCCACAGGGTGATGCTCATAATGATGACACACGAGATTAAGGTGCGGAACA‑GATTTCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                                                              >  2:297568/1‑250 (MQ=255)
        ggTGTGAAACCCAATAGTTGGCGTCGGGACCCAAAGGGTGATGCTAAAAATGATGCCCCCCGAGATAAAGGGGCGGAACA‑GATTCCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑AAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATAAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                                            <  1:22071/251‑1 (MQ=255)
               aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                                     >  1:300393/1‑251 (MQ=255)
                  ccaatgcttggcgcagcgACCCAGAGGGTGAGGCCCAAAGTGATGACCCCCGAGATTAAGGTGCGGAAGA‑‑AATTCCAAGAAAGCACGCGGCCCCCCGCCGCCCCGAATAAGCCGTTGGTGG‑CAAGATCGCCAAACGGCTGTACCTTATTACTCGGCACAGAAGCAGCAAATAAAGGGGTGAGTGGGATAAAAAACAGCGATGCCAATTCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                                                                                                                                  <  1:168595/236‑1 (MQ=255)
                  ccAATAGTTTGCGTCGCGATTTACTGGGGGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATGTCCACGAT‑GCACGGCGACCGCCATCGCCAGGAAGCAGCAGATGGTGG‑CACGATCGCCAAACGGCAGGACCTGATTACCCGGCACACAAACAGCAAGTTAAATGGTGAGTGGGATCAAAATCCGGCTAGCCGATACCACC‑‑GCCGTATGTCCCAGCAGACACCCCGGATCAAAGCCacc                                                                                                                                                                                                                  >  2:25048/1‑249 (MQ=255)
                   cAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                 >  1:279080/1‑251 (MQ=255)
                   cAATAGTTTGCGGCGCGATCCCCAGGGGGATGCTCATAATGAGGACCCCCGAGATTAAGGTGCGGAACA‑GATTTCCCCGAT‑GCACGGCGACCGCCATCGCCACGTAGAGGCCGATGGTGG‑AAAGATGGCCAAACGGCAGCTCCTGATTACCCGGCACACAAACAGGAATTAAAATGGTCAGTGGGAAAAAAATCAGGCTTGCCGAT‑CC‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                 <  1:330950/251‑1 (MQ=255)
                      tAGTTTGCGTCGCGATCCACAGGGAGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                                              >  1:105282/1‑251 (MQ=255)
                                                        aaggatgaCACCCG‑GAATTAGGGGGGGAACAAGATTCCCCCGAT‑GCAGGGCGACCGCCACGGCCACGAAGAAGCCGATGGGGG‑CAAGATCGCCAAACGGCAGCCCCTGATTACCCGGAACAAAAAAAGCAATTAAAAGGGTGAGTGGGAAAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                                                                                                                                            <  2:122162/248‑1 (MQ=255)
                                                            atgaCACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                        >  2:57971/1‑251 (MQ=255)
                                                                 ccccGAGATTAAGGTGCGAAACA‑GATTCCCCCGAT‑GCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                                                                                                                                   <  1:166503/250‑1 (MQ=255)
                                                                      agaTTAAGGTGCGGAACA‑GATTTCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                                               >  1:117812/1‑250 (MQ=255)
                                                                            gagggggggaaacgtttcccACGAG‑GCGCGGCCACCGCCACCCCCACGAAGAAGCGGATGGTGG‑AAAGATCGCAAACCGGCAGCCCCTGATTCCCCGGCCCAAAAACAAAAATTAAAGTGGTGAGGGGAAAAAAAACCAGGCTTGCCGACAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACGTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                                         <  2:72269/233‑1 (MQ=255)
                                                                                    taCA‑GCTTTCCACGAT‑GCACGGCGACCGCCATCGCCAAGAAGAAGCCGATGGTGG‑CAAGATCGCCCAACGGCAGCACATGATTACCCGGCACACAAACAGCCATTACAATGGTGAGTGGGATAAAAATCCGGCTTGCCGATACAACC‑‑GCCGTATGTCCCAGCAGCAACGCAGGTTCAAGACCAATCAGGAAATACTGACAGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTCGCGATCGGCgctc                                                                                                                                                 >  2:49608/2‑247 (MQ=255)
                                                                                          aTTTCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                                                                                                                           >  1:259654/1‑251 (MQ=255)
                                                                                                          ggCGCCCGCCATCGCCACGAAGAAGCCGA‑GGTGGGAAAGATCGCAAAACGGCAGCCCCTGATTCCCCGGCACACAAACAGCAATAAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGA‑ACC‑CACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCGCCGAACTGCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                                                                                                                            <  1:297568/251‑1 (MQ=255)
                                                                                                                                aagCCGATGGTGG‑AAAGATCGCCAAACGGCAGCACCTAATTACCCGGCACAAAAACAGAAATAAAAATGGTGAGTGGGAAAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCAt                                                                                                      <  1:57971/251‑1 (MQ=255)
                                                                                                                                                aGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtttt                                                                                       >  1:119058/1‑251 (MQ=255)
                                                                                                                                                 gATCGCCAACCGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGAAAAAAATCAGGCGTGCCG‑TACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                      <  1:142510/251‑1 (MQ=255)
                                                                                                                                                    cGCCAAACGGCAGCCCCTGATTACCGGGCACACAAACAGCAATTAAAAGGGTGAGTGGAAAAAAAATCAGGCTTGCCAAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                                   <  1:321414/251‑1 (MQ=255)
                                                                                                                                                           cGGCAACACCTGGTTCACCGGCACAAAAACATCAATTAAACTGGCGAGCGGGATAATCATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTAGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCt                                                                            >  2:28184/1‑251 (MQ=255)
                                                                                                                                                             tcagaaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGCGGCATTAGCAGCATCACTGCAGCCGTTTTTACCGCCCGCTGc                                                                          >  2:38706/6‑251 (MQ=255)
                                                                                                                                                                          ccccGTCACACAAACAGCAATAAAAAGGGTGAGGGGAAAAAAAACCAGCCTGGCCGA‑ACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGAAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTGAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCATCTGCAATACACCTTTg                                                              <  1:215070/249‑1 (MQ=255)
                                                                                                                                                                                   acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAc                                                    >  2:108610/1‑251 (MQ=255)
                                                                                                                                                                                   acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAc                                                    >  1:47048/1‑251 (MQ=255)
                                                                                                                                                                                    caAACAGCAATTAAAATGGTGAGAGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcc                                                   >  2:102746/1‑251 (MQ=255)
                                                                                                                                                                                     aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAccc                                                  >  1:209321/1‑251 (MQ=255)
                                                                                                                                                                                                        gAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc                                >  1:56537/1‑250 (MQ=255)
                                                                                                                                                                                                                 aaaaaTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcc                      >  2:53933/1‑250 (MQ=255)
                                                                                                                                                                                                                           cTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCTTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGCTAATCAGCCCCATCACAAAACCCAc            >  1:62872/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ccaccgccgTATGTCCCAGCCGCAACGCCGGATCAAGGCCACTCAGGCCCTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTACACCATCCATGCTGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGc  >  2:176119/4‑251 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
GCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACA‑GATTTCCACGAT‑GCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGG‑CAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGC  >  NC_000913/2173132‑2173610

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