Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,665,6980TC33.3% 44.0 / 10.3 31D38G (GAT→GGT) gadXDNA‑binding transcriptional dual regulator GadX
Reads supporting (aligned to +/- strand):  ref base T (15/5);  new base C (0/10);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 2.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGATTAATGCTCTTTTACGTAGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATAT‑CTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATAATGTGTGATGTTCTACGGGCAGG  >  NC_000913/3665455‑3665945
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       
aaGATTAATGCTCTTTTACGTAGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATAT‑CTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCaaa                                                                                                                                                                                                                                                  <  1:170144/251‑1 (MQ=255)
       aTGCTCTTTTACGTAGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATAT‑CTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAgg                                                                                                                                                                                                                                           >  1:119870/1‑251 (MQ=255)
         gCTCTTTTACGTAGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATAT‑CTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTc                                                                                                                                                                                                                                         >  1:304998/1‑251 (MQ=255)
                             tGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATAT‑CTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGat                                                                                                                                                                                                                     >  2:431641/1‑251 (MQ=255)
                                     gAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATAT‑CTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc                                                                                                                                                                                                              >  2:371035/1‑250 (MQ=255)
                                                      aatagttcggcAACTCCCTAAAAAAAAGTAATAT‑CTCTACTACATTAAGAAACGAGAAATAAAAAAACTTTTTT‑ATGT‑ACATGAAGATTTAAGCCCCCC‑CCCTTGAG‑CT‑TCGGAAAAAATAATAAAG‑G‑CGTTCCCCTT‑GACACAAATAAAAAATTTCC‑AACCCCTTTATCTCCTCAATTTTGG‑C‑TGGC‑ACCCGCAAAACCCAGGTCTCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat                                                                                                                                                                                            <  1:236583/242‑1 (MQ=255)
                                                           gtactcAACTTCCGAAAAAACAGGA‑TATTCTCAAATACAGTGAAGACCGAGTAAGGAA‑AAATTTTTTT‑ATTTTATATGATGTTTAAATGCCTCC‑TCTTTGAG‑AA‑GCGGGAAAAAAAATGAAA‑G‑CGGGTCCTTT‑GACAC‑AAAAAAAAATTTTCCAACACTCTTATCAACT‑AAATTTGGGA‑TTTC‑CCCCGAAAAAACCGGGACACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                                                                                                                       <  2:221477/246‑1 (MQ=255)
                                                               ccaactccttaaaaaaCAGGAAAAT‑CTCTAATAAAGTAATAAAAAAGTAATGAAGATCTTTTTTT‑GTGT‑CCAGGAA‑ACTTAATCCCCCCCTCCTTGAG‑CA‑‑CGGGGAAAAAAATGAAAAG‑CGGGCCCTTT‑GCAACAAATAAAAAATTTAC‑AACACACTTACCTACACAAATTGGT‑C‑TTGC‑ACCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTAtt                                                                                                                                                                                   <  2:291293/240‑2 (MQ=255)
                                                                   acttccgaaaaaacgggaagttgtctACTACAGTGATAAACGAATAAGAAAAAAATTTTTTT‑AGGT‑GAAGGAAGATTTAAT‑CCCCC‑CCCGTGGGACA‑TGGGAAGAAA‑AAAGAAA‑GGTGTGCCCCTT‑GACACAAATAAAAATTTTTA‑AACACAATTATCTCCTCAAATTTGC‑T‑TTGC‑ATCCGAAAAAACCAGTTAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCt                                                                                                                                                                                <  1:237254/228‑1 (MQ=255)
                                                                            aaaCAGGAATAT‑CTCTACTACAGTGATAAACGAGTAAGGAAGAAATTTTTTT‑ATGT‑GAATGAAGATTTAATGCCTCC‑TCCTGGAG‑CA‑TCGGAAAAAATAATGAAA‑G‑CGTGTCCCTT‑GACACAAAAAAAAATTTTTC‑AACACCTTTATCTCCTCGATTTTGG‑C‑TTCC‑ATCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCa                                                                                                                                                                      <  1:323404/251‑1 (MQ=255)
                                                                                    attt‑ctcTACTAAAGTGAGAAACGAGAAATAAAAAAACTTTTTT‑ATGT‑GCAGAAAGATTTAATGCCTCC‑‑CCTTGGA‑CA‑ATGGGAAAAAAAATGAAA‑GGCGT‑TCCCTTTGACACAAATAAAAAATTTCA‑AACACCCTTATCTCCTAAATTTTGT‑C‑TTCC‑ACCAGAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACa                                                                                                                                                              <  2:38352/248‑1 (MQ=255)
                                                                                                   gTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTg                                                                                                                                                >  1:330048/1‑251 (MQ=255)
                                                                                                     aatgaACGAGTAAAGAA‑AAACTTTTTTTATGT‑GCATGAAGATTTAATCCCTCC‑TCCTGGAG‑CA‑TCGGGAAAAATAATGAAG‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑ACCACCTTTATCTCCTCGAATTTGG‑C‑TTGC‑ACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGAt                                                                                                                                              <  1:412398/250‑1 (MQ=255)
                                                                                                         aaCGAGTAATGAAGAAACTTTTTT‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGc                                                                                                                                           >  2:10798/1‑250 (MQ=255)
                                                                                                            gAGTAATGAAGAAACTTTTTT‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCAt                                                                                                                                         >  1:56961/1‑249 (MQ=255)
                                                                                                               tAATGAAGAAACTTTTTT‑ATGT‑GAATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTg                                                                                                                                    >  2:386626/1‑251 (MQ=255)
                                                                                                                 agaaaaaatttttttt‑ATGT‑GAATGAAGATTTATTGCCTCC‑TC‑TTGGA‑CAGTCGGGAGAAATAATGAAC‑G‑CTTGCCCCTT‑CACACGAATACAAATTTTTC‑AACACACTTATCT‑CCCATATTTGG‑CCTATC‑AACCGAAAAACCCAGGTCGCCGCCATTAAAATAGCGATATTCGCCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGAc                                                                                                                                  <  1:185588/241‑1 (MQ=255)
                                                                                                                    aataaaatttttt‑ATGT‑GCAGAAAGATTTAAGGCCTCC‑TGCGTGGA‑CA‑TGGGGAAAAAAAATGAAA‑G‑CGTGTCCGTG‑GACAGGAATACAAAATTTTC‑AGCCCGCTTATGTAGTCGA‑TTTGG‑G‑TTGCCGTCCGAAAAAGCGGGGGCGCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATGCGCAATGAGACAATTCCCATGTAGTGATTGCATAGTTGACTTa                                                                                                                               <  2:390684/244‑1 (MQ=255)
                                                                                                                                t‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatatat                                                                                                                                                                                        >  2:136255/1‑182 (MQ=255)
                                                                                                                                t‑ATGT‑GCAAGAAGATTAAATCCCCCC‑CCCTTGAG‑CA‑TCGGGAAAAATAATAAAA‑G‑CGTGTCCCTT‑GACACAAATAAAAAATTTTC‑AACACACTTATCTCCTCGAATTGGG‑T‑TTGC‑ACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAatatat                                                                                                                                                                                        <  1:136255/182‑1 (MQ=255)
                                                                                                                                           aaaaaTTTAATGCCCCC‑TCCTTGAG‑CA‑TCGGAAAAAATAATGAAA‑C‑CGTGTCCCTG‑AACACAAATAAAAAATTTCC‑AACCCACTTTCCTACTAGAATTTGG‑C‑TTCC‑ATCCACAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACt                                                                                                          <  1:386626/247‑1 (MQ=255)
                                                                                                                                                    aTGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTttatta                                                                                                  >  1:406844/1‑250 (MQ=255)
                                                                                                                                                                     agcatgggaaaaaaaaatGAAA‑G‑CGTGTCCCTT‑GACACAAATAAAAAATTTCC‑AACCCACTTTTCTCCTCGATTTTGG‑T‑TTGC‑ACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACgccgcc                                                                                   <  1:431641/236‑1 (MQ=255)
                                                                                                                                                                                       gAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTg                                                                                                                                    >  1:406455/1‑184 (MQ=255)
                                                                                                                                                                                                        gaacaaaaaaaaaaaaTTTTC‑AACCCACTTATCTTTAAAAATTTGG‑G‑TTCC‑ACCCGAAAACCCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATaaataaat                                                     <  2:261934/237‑1 (MQ=255)
                                                                                                                                                                                                           cacaaaTAAAAAATTTTA‑AACACACTTATCTACCGGATTTTGG‑T‑TTGC‑ACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTaaa                                                <  1:135852/247‑1 (MQ=255)
                                                                                                                                                                                                               aaTACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGa                                             >  1:96637/1‑250 (MQ=255)
                                                                                                                                                                                                                  aCAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTCCTGTTTATTAATGTAGCACGCCGCCCACTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCg                                         >  2:429353/1‑251 (MQ=255)
                                                                                                                                                                                                                            taaaCACACTTACCTACCAGAATTTGG‑C‑TTGC‑ACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACa                                 <  2:119870/248‑1 (MQ=255)
                                                                                                                                                                                                                                              cGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATAATGTGTGATGt              >  1:457266/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            gC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATAATGTGTGATGTTCTAcaggcagg  >  2:177042/1‑251 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       
AAGATTAATGCTCTTTTACGTAGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGGAATAT‑CTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGT‑GCATGAAGATTTAATGCCTCC‑TCCTTGAG‑CA‑TCGGGAGAAATAATGAAA‑G‑CGTGTCCCTT‑GACACGAATACAAAATTTTC‑AACACACTTATCTACTCGAATTTGG‑C‑TTGC‑ATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTGTTTATTAATGTAGCACGCCGCCCTCTTAAGTCAAATTGATTGATAAATAAATTTAAACGAAGCGGGAGAACAGGATGATAATGTGTGATGTTCTACGGGCAGG  >  NC_000913/3665455‑3665945

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