Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913829,7251.G20.0% 67.5 / 20.8 30coding (248/999 nt)ybhGHlyD_D23 family protein YbhG
Reads supporting (aligned to +/- strand):  ref base . (8/16);  new base G (4/2);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCTGATTTCAGCGTTGCCTGCGCCTGGTCGCGCGAGGAGCGGGCATTTTCCAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCAGGCCGTTATCCTGGCGGCTTTGATACCACCAGTAGCCTCCGGCAACCACGGCGGCAAGTACCACTACCGCCAATCCGATCACGACAGGTTTT  >  NC_000913/829477‑829964
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                
gtgCTGATTTCAGCGTTGCCTGCGCCTGGTCGCGCGAGGAGCGGGCATTTTCCAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCCTAGTCATAGGCGGCTTTCGCCTGTTTCACCGCTGCGGCGGGCTGGGCGATTTTTTTTTTGCGGTACCCGGGCAGCATC‑AGGGCCTACTGCGCCTGTGCCACC‑GAAAAACCCGCCTTCGCCTGCATCA‑GGG‑cc                                                                                                                                                                                                                                                >  1:228756/1‑250 (MQ=255)
           aGCGTTGCCTGCGCCTGGTCGCGCGAGGAGCGGGCATTTTCCAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTTACCGCTGCGGCGGCCCTTGCGCTTTCTTTATTGCGATACCCCGCAAGCATCCAGG‑AATACTTCGCCTGTGTCACC‑GCAAAACCCGCTTTTGCCCGCATCA‑GGG‑CAATTTCaata                                                                                                                                                                                                                                       >  1:252931/1‑246 (MQ=255)
             cGTTGCCTGCGCCTGGTCGCGCGAGGAGCGGGCATTTTCCAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGGGCATAGTCCTAGGCGGGCTGCGCCTGTTTTACCGCTGCGGGGGGCTGAGGGGTTTTTTCATTGCGATACCCCGCAAGCCTC‑AGGGCATACTGCGCCTGGTCCACC‑GAAAAACCCCCTTTTGCCTGCATCA‑GGG‑GAATCTCAtccagcc                                                                                                                                                                                                                                   >  1:225438/1‑245 (MQ=255)
                 gCCTGCGCCTGGTCGCGCGAGGAGCGGGCATTTTCCAGGTCATTTGCCGAAAGAGTGCGGCTTTTCCACAACCCGTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCTCTTTCGCCTGCATCA‑TGG‑CAATCTCATACGGCTtgt                                                                                                                                                                                                                                <  2:159361/250‑1 (MQ=255)
                           ggTCGCGCGAGGAGCGGGCATTTTCCAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCg                                                                                                                                                                                                                     <  1:87135/251‑1 (MQ=255)
                           ggTCGCGCGAGGAGCGGGCATTTTCCAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCg                                                                                                                                                                                                                     <  1:83680/251‑1 (MQ=255)
                           ggTCGCGCGAGGAGCGGGCATTTTCCAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGGTATAGAAGTTCTGCGCATAGTCATAGGGGGCTTGCGCCTGTTTTACCGCTGCGGCGGCCCGAGGGATTTTTTCATTGCGATACCCGGCAAGCCTC‑CGGGCATACTGCGCCCGTGCCCCCCGAAAAACCC‑CGTTCGCCTGCCTAA‑GGG‑CAATTTCATACGGCGTGTGATCCACTTct                                                                                                                                                                                                                     >  1:194044/1‑250 (MQ=255)
                                     ggagCGGGCATTTTCCAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCt                                                                                                                                                                                                          <  1:25806/251‑1 (MQ=255)
                                          gggCAGTTTCCAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGT‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGccc                                                                                                                                                                                                     <  2:73886/251‑1 (MQ=255)
                                                    cAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAg                                                                                                                                                                                           <  1:160752/251‑1 (MQ=255)
                                                        tCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTTTGGGCATAGTCATAAGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTTTTCATTGCGATACCCGGCAAGCATC‑CGGTCATAATGCGCCTGTGCCCCC‑GACACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCCTGTGATCCAGTTCCCCCA‑GCACCTGGCCCGCTTTGGTAGCaac                                                                                                                                                                                       >  1:8128/1‑249 (MQ=255)
                                                             tGCCGAAATAGCGCGGCTGTTCCACATCCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCAcca                                                                                                                                                                                   <  2:38960/250‑2 (MQ=255)
                                                                 gAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGTTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGGGCAATCTCATACGGCGTGTGATCCAGTTCGCCCA‑GCTCCTGGCCCGCTTTGATAGCATCACCTTCGt                                                                                                                                                                               <  1:169889/251‑1 (MQ=255)
                                                                          cGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTTCGCATAGTCATAATCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGGTTTTTTTATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTTCCACC‑GAAACACCCGCTTTTGCCTGCATCA‑GGGGCAATCTCCTACGGCTTGGGATCCAGTTCGCCCC‑GCACCTTGCCCGCTTTGGTAGCAACACCTTCGTCCCCCGccc                                                                                                                                                                      >  1:194831/1‑250 (MQ=255)
                                                                           ggtttttCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTGCATTGCGATACCCGGCAGGCATC‑AGGTCATACTGCCCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCATTTCGCCCATGC‑CTGGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGc                                                                                                                                                                    <  1:169325/248‑1 (MQ=255)
                                                                                   cacaACCCTGGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGTAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACCCCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCACACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCCCCTTCGTCCACCGCCAGCGATTCAAc                                                                                                                                                            <  2:133236/251‑1 (MQ=255)
                                                                                       aCCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCAGAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGGGCAATCTCATACGGCGTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCAAcgcg                                                                                                                                                         <  1:153073/251‑1 (MQ=255)
                                                                                                            aGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGGGCAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTTAACGCCCCCCCCAACACGGAAACTa                                                                                                                                      >  1:153566/1‑249 (MQ=255)
                                                                                                                               ggcggcTTGCGCCTGTTTCACCGCTGCGGCGGCCTGCGCGCTTTCTTCATTGCGACACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGGGCAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATcc                                                                                                                 >  1:14898/1‑251 (MQ=255)
                                                                                                                                   gcTTGCCCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATTCCCGGCAAGCATC‑AGGTCATACTGCCCCTGTGCCCCC‑GAAACCCCCCCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCC‑GCCCCTGGCCCGCTTTGATAGCATCACCTTCCTCCACCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGtt                                                                                                            <  1:180934/251‑1 (MQ=255)
                                                                                                                                             gTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATGCTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGGGCAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCCCCTTCGTCCACCGCCCGCGCTTCAACGCGCCCCCCAACACGGAAACTCAGATTTACCGTACGAAAATCCTCGTTGCCATAc                                                                                                     >  2:70130/1‑249 (MQ=255)
                                                                                                                                              tttCCCCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATCCCCGGCAAGCATC‑AGGTCATACTGCGCCTGGGCCACC‑CAAACCCCCGCTTTCGCCTGCATCA‑GGG‑CCATCTCATACCGCTTGTGTTCCAGTTCGCCCA‑GCCCCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCCACCCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGc                                                                                                  <  1:100101/250‑1 (MQ=255)
                                                                                                                                                                    gCGATGTCTTCATTCCGATCCCCGGCAAGCATC‑AGTTCATACTGCGCCTGTCCCCCC‑GAACCACCCGCTTTCGCCTGCATCA‑GGG‑AAATCTCATCCGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCAGGCCGTTATCCTGGCGGCt                                                                           <  2:13545/251‑1 (MQ=255)
                                                                                                                                                                      gATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCAGGCCGTTATCCTGGCGGCttt                                                                         >  1:186494/1‑251 (MQ=255)
                                                                                                                                                                      gATTTCTTCATTGCGAGCCCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTTTGATCCAGTTCGCCCA‑GCCCCTGGCCCGCTTTGATAGCATCACCTTCGTCCCCCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCAGGCCGTTATCCTGGCGGCttt                                                                         <  2:206956/251‑1 (MQ=255)
                                                                                                                                                                              cATTGCGATCCCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCAGGCCGTTATCCTGGCGGCTTTGATaccac                                                                 <  1:204608/251‑1 (MQ=255)
                                                                                                                                                                                    gATACCCGGCAAGCATT‑AGGTCATAATGCGCCTGTTCTACC‑AAACCACCCGCTTTCGCCTGCATCATGGG‑AAATTTCATACGGCTTGTGATCCAGTTCTCCCC‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTTCACCGACAGCGATTCACCGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCACGCGGTTATCCTGGCGGCTTTGATatcccacgt                                                             >  2:162027/1‑242 (MQ=255)
                                                                                                                                                                                    gATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCAGGCCGTTATCCTGGCGGCTTTGATACCACCAGTAg                                                           <  1:193107/251‑1 (MQ=255)
                                                                                                                                                                                                         tCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCAGGCCGTTATCCTGGCGGCTTTGATACCACCAGTAGCCTCCGGCAACCACggcggc                                       >  1:157421/1‑251 (MQ=255)
                                                                                                                                                                                                              actgtgcctgttcacccgactcacCCGCTTTCGCTTGCATGA‑GGG‑AAATCTCATACGTCTTGTGACCCAGTGCGCCCA‑GCACCGGGCCCGCGTTAAGAGAATCACATCCGACCACCGCCAGCGATTCAACGCGCCCACCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCAGGCCGTTATCCTGGCGGCTTTGATACCACCAGTAGCCTTCGGCATCCACGGCGGCAAg                                    <  2:194410/229‑1 (MQ=255)
                                                                                                                                                                                                                                               ccTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCAGGCCGTTATCCTGGCGGCTTTGATACCACCAGTAGCCTCCGGCAACCACGGCGGCAAGTACCACTACCGCCAATCCGATCACGACAGGtttt  >  1:194379/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                
GTGCTGATTTCAGCGTTGCCTGCGCCTGGTCGCGCGAGGAGCGGGCATTTTCCAGGTCATTTGCCGAAATAGTGCGGCTTTTCCACAACCCTTGCTGGCGGTTATAGAAGTTCTGCGCATAGTCATAGGCGGCTTGCGCCTGTTTCACCGCTGCGGCGGCCTGAGCGATTTCTTCATTGCGATACCCGGCAAGCATC‑AGGTCATACTGCGCCTGTGCCACC‑GAAACACCCGCTTTCGCCTGCATCA‑GGG‑CAATCTCATACGGCTTGTGATCCAGTTCGCCCA‑GCACCTGGCCCGCTTTGATAGCATCACCTTCGTCCACCGCCAGCGATTCAACGCGCCCCCCAACACGGAAACTAAGATTTACCGTACGAATATCCACGTTGCCATACAGCGTCAGGCCGTTATCCTGGCGGCTTTGATACCACCAGTAGCCTCCGGCAACCACGGCGGCAAGTACCACTACCGCCAATCCGATCACGACAGGTTTT  >  NC_000913/829477‑829964

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