Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,521,4060AT35.3% 51.3 / 12.1 34F205L (TTT→TTA) yjhFputative transporter YjhF
Reads supporting (aligned to +/- strand):  ref base A (3/19);  new base T (12/0);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 8.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑AAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGA  >  NC_000913/4521168‑4521644
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
aGATCGCGAAGAGCAGTGAAATCAGCATGGCTATCGTTGAATTACCGAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGTCTCTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACATCCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTGGCTCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGa                                                                                                                                                                                                                                        <  1:419602/251‑1 (MQ=255)
    cgcgAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCACTTTTCCCCTCTTTTTTTTAGGCGGCCgggaaa                                                                                                                                                                                                                                     >  2:26479/1‑246 (MQ=255)
                        gCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGGCGCCCCGGGGGAGGCTATATTCTCCCCTTTTTTTAAAGCGCCCCCAAAAATGCGCCGGCAGAAattattt                                                                                                                                                                                                                >  1:10035/1‑250 (MQ=255)
                            ggCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAAt                                                                                                                                                                                                             <  1:26479/250‑1 (MQ=255)
                                            ccTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGCGAGAGGGTATATTCCCCTCTTTTTTTAAGGCGGCCGGGATAATGCGCTGGCAGAATTTATTTAATAACGGGCCGC‑CAAGaa                                                                                                                                                                                            >  1:373063/1‑250 (MQ=255)
                                                   aCAGTATCACGTTGAAAAATAAGCCGGAACCCGCGGTACTTTGGT‑CA‑TTAGC‑CATCTGGATACTGTTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTGTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGTTGGCTGAATTTATTCAATACCGGTCCGC‑TCAGAATAAAGGtt                                                                                                                                                                                    <  2:285474/251‑1 (MQ=255)
                                                    cAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTTAAAGCGTCCGGGATAATGTGCTGGCAGAATTTTTTTAATACCCGGCCGC‑CAAGaaaaaaaggtg                                                                                                                                                                                   >  1:37496/1‑246 (MQ=255)
                                                        aTCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGGCGCGGCGAGAGAGGCTATATTCCCCTCTTTTTTTAAGGCGGCCGGTATAAAGCGCTGGGAGGATTTATTTAATAAAGGGCCGCCAAAGA‑TAAAAGTTgggt                                                                                                                                                                                >  1:10016/1‑248 (MQ=255)
                                                           agtttGTAAATTAAGCCTGACTCCGCGTC‑CTTTTGGGCAGTG‑GT‑CATCTGGATGATGGGGACAACCGCCATCAGTATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGTCTGGCATTTCTCTTTCTCTTAATCGGCGGGGCGCGCCGATAGATGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATTCTCTGGCATAATTTATTCAATACCGGGCCGC‑ACAGATTAAAGGTTGGTATGcc                                                                                                                                                                            <  2:372844/248‑1 (MQ=255)
                                                                  aaaaTAAGCCAGACTCCGCGGCGCTTTTTG‑CA‑TG‑GTGCATCGGGATGAGGGTGACAACCGCCATCAGGAGCACCGGCAGCAGTGCCGTGAGAAAACTGATCCCAAAGCCGGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGcagtcaat                                                                                                                                                                      <  2:305167/250‑8 (MQ=255)
                                                                     aTAAGCCGGAATCCGCGGCGCTTTTGG‑CA‑TGAGT‑CATCGGGATGATGGGGACAACCGCCATCAGGATCCCCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATa                                                                                                                                                                  <  2:336717/251‑1 (MQ=255)
                                                                                  cgcgGCACTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTa                                                                                                                                                    <  1:386802/251‑1 (MQ=255)
                                                                                         cTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGGCGCGCCGGGAGAGGCTATATTCCCCTCTTTTTTTAAGGCGGCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑CAAGAATAAAAGGTGGAATGCCGACAATGATAACGTATAT‑CAGTACCAATcc                                                                                                                                             >  1:303517/1‑251 (MQ=255)
                                                                                             tGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTAc                                                                                                                                          <  1:416707/250‑1 (MQ=255)
                                                                                                                    atgGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGGGAGAGGAGATATTCCCCTCTTTTTTTAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTTAATACCGGGGCGC‑AAAGAATAAAGAGTGGGATGCCGACAATGATAACGCATAT‑CAGTAACAATCCTACATCAGCGCCAAATTCCCTGGc                                                                                                                     >  1:95432/1‑251 (MQ=255)
                                                                                                                          aCAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAAc                                                                                                               <  1:28648/251‑1 (MQ=255)
                                                                                                                              ccGCCACCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGGCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTTAAAGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAAGGTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCCCCATATTCCCTGGCGATAAAAGTc                                                                                                           >  1:115853/1‑251 (MQ=255)
                                                                                                                                   aTCAGGATCACCGGAAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGAATTTCGCTTTCGCTGAATCGGCGGGGCGCGCCGGGAGAAGGAATATTCCCCTCTTTTTTTAAAGCGTTCCGAATAAACCGCTGGCAGAATTTATTCAATACCGGGGCGC‑CAAGAATAAAGGGTGGGATGCCGACAATGATAACGGATAT‑CAGTACCACACCTACAACAGCGCCCTATTTCCTGGGGATaaaaatctgggc                                                                                                      >  1:270309/1‑241 (MQ=255)
                                                                                                                                                  agcaTTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGGGAGAGGCTATATTCCCCCCTTTTTTTAAGGCGTCCCGAATAATGCGCTGGCCGAAATTTTTCAAAAACCGGCCGC‑AAAGAATAAAGGTTGGGAAGGCGACAATGATAACGTATAT‑CAGTAACAAACCTAAATCAACGGCAAATTTCCTGGCGATAACAGTC‑GGGCCGGGGTGCggcggcc                                                                                       >  1:31266/1‑249 (MQ=255)
                                                                                                                                                    caTTGCTGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTTAAAGCGTCCGGGATAATGCGCTGGCAGAATTTTTTCAATAAACGGCCGC‑ACAGAATAAAAGGTGGAATGCCGTCAAAGATACCGTA‑ATACAGTACCAATACTCAATCAGCGCCATATTTCCTGGCGGTAAACGTCTGGGCCGGGg‑gcggcggccaa                                                                                     >  1:119402/1‑247 (MQ=255)
                                                                                                                                                             gagaAAACCGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATAGTCCCCTCGTTTTTAAAGGCGTCCGGAATAATGCGCTGGCATAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATTCCGTATAT‑CAGTGCCTATCCTTCATCAGCGCCATATTCCTTGGCGATAACAGTC‑GGGCCGGGTTGCGGCGGC‑AAAAAACCgtga                                                                           <  2:259539/250‑2 (MQ=255)
                                                                                                                                                                     tGATCCCCAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTCTATTCCCCTCTTTTTTAAAAGCGTCCGGGATAATGCGATGGCAGAATTTATTCAATACCGGGCCCC‑CAAAAATAAAGGTTGGGATGCCGACAATGTAAACGTATAT‑CAGTACCAATCCTACATCAGCGCCATAATTCCTGGCGATAAAAATC‑AGGCCGGGGGGCGGCGGC‑AAAAAAACCTGTGCgccacac                                                                  >  1:362683/1‑245 (MQ=255)
                                                                                                                                                                      gATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCGCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑AAAAAACCGTGTGCGACCAg                                                                   <  2:115853/249‑1 (MQ=255)
                                                                                                                                                                         cccAAAGCCTGGCATTTCGCTTTCGCTGAATCTGCTGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑AAAAAACCGTGTGCGACCAGCAAAc                                                              <  2:37496/251‑1 (MQ=255)
                                                                                                                                                                          ccAAAGTCTGGCATTTCGCTTTCGCTGAATCGGCGGGGCGCCCCGAGAGAGGCTATATTCCCCTCTTTTTTTAAGGGGCCCGGAATAATGCGCTGGCAGAATTTTTTTAAAACCGGGCCCC‑AAAGAATAAAAGTTGGGATGGCGACAATGATAACCAATAT‑CAATAACAAAACCACATCATCGCCATATTCCCTCGCGGTAAAAGTC‑GGGCCGGGGGGCGGCGGCCAAAAAAcctggggggcccccaaacc                                                              >  1:116972/1‑234 (MQ=255)
                                                                                                                                                                                           gctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATc                                                                                                                                              <  1:235870/155‑1 (MQ=255)
                                                                                                                                                                                           gctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑AAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGcacac                                            <  1:194913/251‑1 (MQ=255)
                                                                                                                                                                                                             ggTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑ACAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑a                                                                                      >  2:437124/1‑191 (MQ=255)
                                                                                                                                                                                                             ggTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑a                                                                                      <  1:437124/191‑1 (MQ=255)
                                                                                                                                                                                                              gTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATAGTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑AAAAATCCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGatat                         <  1:389956/251‑1 (MQ=255)
                                                                                                                                                                                                                                  atTCCCCTCTTTTTTAAAGGCGGCCGGAATAATGCGCTGGCAGAATTTATTCAATCCCGGGCCGC‑AAAGAAGAAAGGTTGGAAAGCCGACAATGATACCGTATAT‑CATTACCATGCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑AAAAAACCGTGAGCGACCAGCAATCCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATa      <  2:270309/250‑1 (MQ=255)
                                                                                                                                                                                                                                     ccccTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑AAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGa  <  1:202232/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tttAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑AAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATAttttt                     >  2:1180/1‑223 (MQ=255)
                                                                                                                                                                                                                                              tttAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑AAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATAttttt                     <  1:1180/223‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
AGATCGCGAATAGCAGTGAAATCAGCATGGCGATCGTTGAATTACCTAAGAACAGTATCACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGG‑CA‑TGAGT‑CATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGC‑AAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATAT‑CAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTC‑GGGCCGGGGTGCGGCGGC‑AAAAAACCGTGTGCGACCAGCAAACCAGAAAGCATCGGCACACACATAAACATCGGTGATATTTTTGCTTCACGGGCAATAGCGA  >  NC_000913/4521168‑4521644

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