Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,008,3730GT32.3% 28.1 / 13.3 31G33C (GGC→TGC) yedEputative selenium transporter YedE
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base T (10/0);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 4.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCAT  >  NC_000913/2008139‑2008623
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                               
aTGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTAATTCCTATTTTGGCATTACTGGCA‑CC‑tt                                                                                                                                                                                                                                               >  1:74673/1‑251 (MQ=255)
                          gagaAATTTGATCTAACTATTCTTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTTTCTCTTCTTTCTTTTTTGGCCTTTCTTGGA‑CC‑CTTTTGGCTGTC‑ACCGGTGTAttttacc                                                                                                                                                                                                                    >  1:114430/1‑247 (MQ=255)
                                   gATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTTTTTTCTTCACTTTCCATTTTTGCATTACTGGCC‑CC‑TTTTTGGCTGTCCACGGG‑GAATTTACC‑CGTTGGggcg                                                                                                                                                                                                          >  2:109515/1‑251 (MQ=255)
                                        aaCTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGcc                                                                                                                                                                                                       <  1:305302/249‑1 (MQ=255)
                                          cTATTATTTTTGTGAATTATTACCCCGAAAATCAATTTGCCTGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTTGCCTTTATGGCC‑CCCTTTTGGG‑TTTC‑ACGGGTGCATTTAcc‑ctcttggggggccacgc                                                                                                                                                                                                   >  2:173207/1‑235 (MQ=255)
                                                       gAATTATTACCCCGAAACTAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGATCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTTCCTCCTCTTTTTGCCTTACTTGCC‑CC‑TTTTTGGGTGGC‑ACGGGGGAATTTTTC‑CGTTGGGGGGGcccggcccctccactgttt                                                                                                                                                                                      >  2:479057/1‑234 (MQ=255)
                                                       gAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑c                                                                                                                                                                                     <  2:131848/251‑1 (MQ=255)
                                                                   cGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAAACATCTCAGTTCAGCGGCTACATTCATGCCATGGCCGCACTTCAAACACGCCCGGTTGATTACATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGCCTTACTCTCTACATACCATTTTTGCAATTCTGGCA‑CC‑TTTTTGGCAGCT‑ACGGGTGGATTTACC‑CGTTGGGGGGGCCAGCCCCTGCACCTGTT‑CGGCGtacctgct                                                                                                                                                                         >  1:503469/1‑244 (MQ=255)
                                                                          aaatTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTTTCTCTACTTACTATTTTGGGATTACTGGCA‑CC‑TTTTTGGCTGTC‑CCGGGGGAATTTACC‑CGTTTGGGCGGCCCGCTCCTGCCACTGTT‑CGGCGTCCCTGCTGCagagg                                                                                                                                                                  >  1:232491/1‑250 (MQ=255)
                                                                                    ggTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTTTTTTTGCCTTTCTGGCA‑CC‑TTTTTGGGTGTC‑CCGGGTGGATTTTCC‑CGTTCGGGCGGCCAGCTCCTGCCACTGTT‑GGGCGTCCATGCTGGAGAGTggggggtctt                                                                                                                                                        >  2:449686/1‑245 (MQ=255)
                                                                                     gtgtATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTCCTTACTATTTTTGGATTACTGGCA‑CC‑TTTTGGGGTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACCGTT‑CGGGGGCCATGCTGAAGAGGGGGGTTACttt                                                                                                                                                       >  1:181032/1‑251 (MQ=255)
                                                                                       gtATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCTTCCATGCTGAAGAGTGGGGTTACTTTa                                                                                                                                                      <  1:203962/250‑1 (MQ=255)
                                                                                              tCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTAt                                                                                                                                               <  1:110163/250‑1 (MQ=255)
                                                                                                  gACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGGCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATc                                                                                                                                          <  1:192350/251‑1 (MQ=255)
                                                                                                                            cAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTACATTCTGGGCCCCCATCACTGCGGTCATCGCGGCGGGGATTTTCTCCACTTTCTATTTTTGGATTACTTGCA‑CC‑TTTTTGGCTGTC‑ACGGGGGGATTTTCC‑CGTTTGGGCGGCACGCTCCTACAACTGTT‑CGGCGGCCATGCTTCAGAGTGGGGGTACTTTTAAATTATCCCTCTGGGAGGATTGCCATTAAACCGCacc                                                                                                                >  2:452294/1‑249 (MQ=255)
                                                                                                                            cAGAGGCTACATTCATGTAATGGCCGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCAATCACTGCGGTCATCGCGGCGGGGATTTTCTCTTATTACCATTTTTGCATTACTTGCC‑CC‑TTTTTGGCTTTT‑ACGGGGGAATTTTCC‑CGTTGGGGGGGCCAGCTCCCGGAACTGTTTGGGGGTCC‑TGCTGCAGAGTGGGGTTTATTTTAAATTAACCATCTGGAAgggtacgcatatacccgcatc                                                                                                                >  2:354705/1‑232 (MQ=255)
                                                                                                                                         catgtcatgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGgatgatga                                                                                                   <  2:232491/251‑1 (MQ=255)
                                                                                                                                              catgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTAATTTCTCTACTTACTATTTTTGCATTTCTGGCA‑CC‑TTTTTGGGTTTC‑CCGGGGGGATTTTCC‑CGTTGGGGGGGCCAGCTCCTGCAACTGTT‑CGGCGGCCCTGCTGAAGAGTGGGGTTACTTTATAATTATCCATCTGGAAAGATCGCCCTTAACCACCCTCGACGGGATGATGATcccc                                                                                              >  2:120963/1‑249 (MQ=255)
                                                                                                                                                          tCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTgg                                                                                   <  1:449686/250‑1 (MQ=255)
                                                                                                                                                            aaaCACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑CCGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGGT‑CGGCGGCCATGCTGAAGAGTGGGGGTACTTTTAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGCCGGGGTGATGATCCTCGGGATGTTTggcgg                                                                                >  1:10320/1‑251 (MQ=255)
                                                                                                                                                                      ggTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGCTTCTCCCTACTTCCTATTTTGGCATTACTGGCC‑CC‑CTTTGGGCCGTC‑AAGGGGGAATTTAAC‑CGTTGGGGCGGCCAGCCCCTGCCACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGGTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCCACGGGATGGTGATCCTCGGTATGTTTGGCGGCTGCTTTGcc                                                                      >  1:267958/1‑251 (MQ=255)
                                                                                                                                                                        ttGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCATTAATTAATATTTTGGCAGTACTGGCA‑CC‑CTTTTGGCTGTC‑ACCAGGGAATTTACCACGTTGGG‑CGGGCAGCTCCTGCCACTGTT‑TGGGGTCCAAGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGGTCGCCACTAAACCCCCTCCACGGGTAGATGATCCTCGGTATGTTTGGCGGCTTCTTTGcagc                                                                    >  1:43199/1‑248 (MQ=255)
                                                                                                                                                                                      ggCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGGGAATTTACC‑CGTTGGGGCGGCCCGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGGTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGccc                                                      >  1:387138/1‑250 (MQ=255)
                                                                                                                                                                                          cccATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTTCTATTTTTGCATTACTGGCAACC‑TTTTGG‑CTGTC‑ACGGGTGAATTTAAC‑CGTTGGGGCGGCCAGCTCCTGCCACTGTT‑CGGCGTCCATGCTGAAGAGTGGAGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAAACCGCATCGAAGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCCGTGGGCcaaca                                                   >  1:28102/1‑250 (MQ=255)
                                                                                                                                                                                                        cATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGgct                                                                              <  1:138950/209‑1 (MQ=255)
                                                                                                                                                                                                        cATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCAACTGTAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTGGTCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGa                                     <  1:354705/250‑1 (MQ=255)
                                                                                                                                                                                                                     aTTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCg                       >  2:140512/1‑251 (MQ=255)
                                                                                                                                                                                                                       tctctcTACTTACGTTTTTGGCGTTACTGGCA‑CC‑TTTGTGCCTGTA‑ACGGGTGAATTTTCC‑CGTTGGGGCGGCCAGCTCCGGCAACTGTT‑CGTCTTCCATGCGGAGGAGTGGTGTTACTATAAAATTATCCATCTGGAAGGATCGCCATTAGCCCGCATCGACGGGATGATGATCCTCTGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGGGGGCTAACAATGTGAAGCTGCGAATGGCGGGCAGCCGt                      <  2:452430/250‑1 (MQ=255)
                                                                                                                                                                                                                          ctctACTTACTATTTTGGCATTACTGGCA‑CA‑TTTTGGGCTGTC‑ACGGGAGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCCACTGTT‑CGGCGTCCCTGCTGAAGAGTGGGGTTACTTTAACATTCTCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGACTGCCGCGCAGCCGTATcc                  >  2:348004/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tactATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCAtt             <  1:120963/249‑1 (MQ=255)
                                                                                                                                                                                                                                  actATTTTGGCCTTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGACGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCACCAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCTCATTATg          >  2:358431/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ggCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCcat  <  1:56666/251‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                               
ATGTCGCCAGCCTTTCAAGGTTCAGAGAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCA‑CC‑TTTTGGGCTGTC‑ACGGGTGAATTTACC‑CGTTGGGGCGGCCAGCTCCTGCAACTGTT‑CGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCAT  >  NC_000913/2008139‑2008623

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