Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2750TC59.3% ‑3.2 / 16.5 27intergenic (‑1005/‑707)yohF/yohPputative oxidoreductase YohF/putative membrane protein YohP
Reads supporting (aligned to +/- strand):  ref base T (8/3);  new base C (3/13);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 1.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑TTTGCTCGCGTTCCGCCGCCCGTGCTTTAA‑CCGTA‑CCCAGGCGGGCAGTAACCTCCGCC‑CGATTCTTTCCCC‑AGATATCGAGATGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTG  >  NC_000913/2228027‑2228504
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        
gcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑GTTGCTCGCGGTCCCCCGCCCGGGCTTTAA‑ACGTATCCCAG‑CGGGCGGTAACCTcc                                                                                                                                                                                                                                        >  1:229995/1‑251 (MQ=255)
                                  acCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGCCCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑TTTGCTCGCGTTCCGCCGCCCGTGCTTTAA‑CCGTA‑CCCAGGCGGGCAGTAACCTCCGCC‑CGATTCTTTCCCC‑AGATATCGAGATGCCAGc                                                                                                                                                                                                    <  2:101314/251‑1 (MQ=255)
                                            gaagaaGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑TTTGCTCGCGTTCCGCCGCCCGTGCGTTTA‑CCGTA‑CCCCGGCGGGGAGGAACCTCCGCC‑CGATTCTTTCCCC‑CGAAATCGAGAGGCCCGCccgccggtt                                                                                                                                                                                           >  1:168530/1‑247 (MQ=255)
                                                            tCTGATATAGCTGGAGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGG‑GTTTCCCGCGTTTCGCCCCCCGCGATTTAAACGGTA‑CC‑AGGCCGGCAGGCACCTCCGCC‑CGATTATTTTCCC‑CGATAACCAGATGCCAccccccccttacgcccacataaacgttg                                                                                                                                                                          >  1:235428/1‑227 (MQ=255)
                                                                atatAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑GTTGCTCGCGTTCCGCCGCCCCCGCTTTTA‑TCGTA‑CCCCGGCGGGCAGTAACCCCCGCC‑CGATTTTTTCCCC‑AGATATCGGGATGCCAGCCCCCCGTTCA‑ACCAAAAGTAccccttgtgg                                                                                                                                                                     >  1:207462/1‑243 (MQ=255)
                                                                                  cgtaaTAATGGTGTTCTGTTCTTTTTCTATTTGCTGCAAGCCGGTGTTTACCGCCGCCGGGGTTTGCCCCTCCCAGTACAGGGGGCCTACGCTGCCAGCCAGCAATTGGCGGG‑TTGGCGCGCGTTCCGCCGCCCTTGCTTTAA‑CCGTC‑CCCAGGCGGGAAGTACCCCCCCCC‑CGATTCTTCCCCC‑AGATATCGAGAGGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTAc                                                                                                                                                   <  1:89519/247‑1 (MQ=255)
                                                                                           ggTGTTCTGCTCTTTTTCTATTTGCTGCAAGCCCGAGTTTAACGCCGCCTGGGGTTGACACTCCCAGTACAGGCGGGCTAAGCTTCCAGCCAGCAATTGGCGGG‑TTTACTCGCGTTCCGCCGCCCGTGCTTTAA‑ACGCA‑CCCAGGCGGGCAGTAAACTCCGCC‑CAATTCTTTCCCC‑AGATTACGGGCTGCCAGCCCACCGTTAAAGCTAAAAGAACC‑TTTGTGTACCCCGGGGCggccagaacttaggccg                                                                                                                                          >  1:284882/1‑236 (MQ=255)
                                                                                               ttCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑TTTGCTCGCGTTCCGCCGCCCGTGCTTTAA‑CCGTA‑CCCAGGCGGGCAGTAACCTCCGCC‑CGATTCTTTccc                                                                                                                                                                                                                        <  2:258331/171‑1 (MQ=255)
                                                                                                     tctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑TTTGCTCGCGTTCCGCCCCCCGTGCTTTAA‑CCCTA‑CCCAGGCGGGCCGTAACCTCCCCC‑CGATTCTTTCCCC‑CGATATCGAGATGCCAACCCGCCGTTAA‑GCCAAAAATACCCTTGGTGTACCACCGGCCGGTCGTACCCGCGGCCGGGTCggtc                                                                                                                                  >  1:90534/1‑246 (MQ=255)
                                                                                                         tttCTATTGCCTGCAAGCCCGTGTTTAACCCCGCCGGGGTTTGCCCCTCCCATTACAGGCGGGCTCCGCTGCCAGCCAGCAATTGGCGGG‑TTTGCTCGCTTTCCGCCGCCCTTGCTTTAC‑CCGTC‑CCCAGGCGGGCAGTACCCCCCGCC‑CGATTCTTCCCCC‑AGATATCGAGATGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGc                                                                                                                             <  1:9492/250‑1 (MQ=255)
                                                                                                          ttCTATTTGCTGCAAGACCGTGTTTAACGACGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑ATTTCTCGCGGTCCGCCCCCCGTGCTTTTA‑CCCGG‑CCCAAGCGGGGAGTAACCCCCCCC‑CGATTTTTTTCCC‑CGATATCGAGAGGCCAGCCCGCCGATAA‑AACAAAAAAACCGTTGGTGTACCCCCGGGCGGGCGTACCCGCGGCCGGATCGGTCagagga                                                                                                                            >  1:112363/1‑248 (MQ=255)
                                                                                                                                    aCGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑TTTTCTCGCGTTTCGCCCCCCCTGCTTTAA‑CCCTA‑CCCAGGCGGGGCGGAACCCCCGCC‑CGATTTTTTTCCC‑AGATATTGAGAAGCCCGCCCCCCCTTAA‑ACCCAAAATACCCTTGGGGTACCCCCGGGCGGGCGTACCTCCGGCCGGGTCGTTTAGAGCAAAAGGCCCCCCTAAGACTTCTGCCggc                                                                                                 >  1:69676/1‑249 (MQ=255)
                                                                                                                                        cgccTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑TTTGCTCGCGTTCCCCCGCCCGTGCTTTAA‑CCGTA‑CCCCGGCGGGCCGTAAACTCCCCC‑CGATTCTTTCCCC‑CGATATCGAGATGCCAGCCCGCCGCTAA‑GCCCAAAATAACCTTGGTGTACCACGGGCCGGGCGTACCCGCGGCCGGGTCGGTCAGAGCAAAAAGCCCCCTTAAGCCTTCTGCCGGCAttt                                                                                             >  1:301407/1‑251 (MQ=255)
                                                                                                                                                  tcccccccccATTACCGGCGGGCTACGCTGCCAGCCAGCAAGGGGCGGT‑TTTGCTCTCGTCCCGCCGCCCGTGTTTTAA‑CCGTC‑CCCAGGCGGGCAGTACCCCCCGCC‑CGATTCTTCCCCC‑AGATATCGAGATGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCt                                                                                    <  1:27634/243‑1 (MQ=255)
                                                                                                                                                         cccgttcagggcgggcgACTCTGCCACCCAGCATTTGGGGGG‑TTGGCTCGCGTTCCGCCGCCCTTCCTTAAC‑CCGAC‑CCCGGGGGGGCAGTACCCCCCGCC‑CGATTTTTCCCCC‑AAATATCGAGAGGCCGGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATcc                                                                            <  2:128309/242‑1 (MQ=255)
                                                                                                                                                                 gggcggccgccccggccagccaG‑AATTGGGGGGGTTTGCTCGCTTCCCGCCCCCCGGGCTGTAA‑CCGTC‑CCCGGGCGGGCGGTACCCCCCTCC‑CGTTTTTTTCCCC‑AGACCTCGGGATGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGTGCCGGTCGTACCTTCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTGCTGCCGACATTTTTTGCCGCTCCATATCCGCCGaaaa                                                                    <  1:143199/237‑1 (MQ=255)
                                                                                                                                                                 gggCGGGCTACGCTCCCGCCCAGCAATGGGGGGT‑TTTCCTGGCTTCCCGCCCCCCGTTTTTTAC‑CCGCC‑CCCGGGGGGGCATTACCCCCCCCC‑CGTTTTTTTCCCC‑AGATATCGAGATGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGaaaa                                                                    <  2:10948/250‑1 (MQ=255)
                                                                                                                                                                       gccaCGCTCCCAGCCAGCAATTGGGGGT‑TTTGCTCGCGTGCCGCCGCCCGTGCTTTAC‑CCGGA‑CCCGGGCGGGCAGAACCCCCCGCC‑CGATTCTTCCCCA‑AGATATGGAGATGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGAt                                                              <  1:153995/248‑1 (MQ=255)
                                                                                                                                                                                  ggccagaAATTGGCGGG‑TTGGCTCGTTTCCCCCCGCCCGTGCTTTAC‑CCGCA‑CCCAGGCGGGCAGAACCCCCCGCC‑GGTTTCTTCCCCA‑AAATATCGAGATGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCAt                                                   <  1:309659/250‑1 (MQ=255)
                                                                                                                                                                                     ccgcaattgcggtg‑TTTGCGCGCGTGCCGCCGCCCGTGCTTTAC‑CCGTC‑CACGGGCGTGCAGTACCCCCCGCC‑CGATTCTTCCCCC‑AAATATCGAGTTGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCTCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCt                                                <  1:292815/238‑1 (MQ=255)
                                                                                                                                                                                             ggCGGG‑TTTGCTGGCTTTCCCCCGCCGGTGTTGTAA‑CCGTA‑CCCAGGGGGGGAG‑AACCCCCCCCACTATGCTTCCCCC‑AAATAGCGGGATCCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCGTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCa                                        <  1:301596/251‑1 (MQ=255)
                                                                                                                                                                                                       gccggggtgccgcccCCCGTGCTTAAC‑CCGTA‑CCCGGGGGGGAAGTAACCTCCGCC‑CAATTCTTCCCCC‑AGATATCGAGATGCCATCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTc                               <  2:254423/242‑1 (MQ=255)
                                                                                                                                                                                                       gcccgctttCCGCCGCCCGTGTTTTAA‑CCGTC‑CCCAGGCGGTCAGAACCCCCCCCC‑CGATTCTTCCCCC‑AGATATCGAGATCCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCtt                                <  1:253457/243‑1 (MQ=255)
                                                                                                                                                                                                                   ccgccCGTGCTTTAA‑CCGTA‑CCCAGGCGGGCAGTAACCTCCGCC‑CGATTCTTTCCCC‑AGATATCGAGATGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTc                    >  1:35980/1‑250 (MQ=255)
                                                                                                                                                                                                                     ccccccggcttttacccgcccccAGGGGGGCAGTACCCCCCCCC‑CTTTTTTTCCCCC‑AGATACCGAGAGGCCCGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCa                   <  2:168530/231‑1 (MQ=255)
                                                                                                                                                                                                                             tttAA‑CCGTA‑CCCAGGCGGGCAGTAACCTCCGCC‑CGATTCTTTCCCC‑AGATATCGAGATGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGAGGCAACTGCCTTCGCCTGGGCTTCAGCTAGATGa          >  2:245278/1‑250 (MQ=255)
                                                                                                                                                                                                                                     gca‑CCCGGGCGGGCATTACCCCCCCCC‑CG‑TTCTTTCCCCCAAATATCGAGAGCCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTg  <  2:243841/249‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                        
GCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑TTTGCTCGCGTTCCGCCGCCCGTGCTTTAA‑CCGTA‑CCCAGGCGGGCAGTAACCTCCGCC‑CGATTCTTTCCCC‑AGATATCGAGATGCCAGCCCGCCGTTAA‑GCCAAAAGTACCGTTGGTGTACCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTAGCTGAATGCGCTG  >  NC_000913/2228027‑2228504

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