Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913377,4500TA65.0% 3.9 / 16.4 20intergenic (+139/+85)yaiL/frmBDUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base A (11/2);  total (13/7)
Fisher's exact test for biased strand distribution p-value = 2.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCT  >  NC_000913/377222‑377694
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ataGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACg                                                                                                                                                                                                                                >  1:23646/1‑251 (MQ=255)
                                gCAGGATGAAGACGACCCGTAGGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTTCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCg                                                                                                                                                                                                <  1:137955/251‑1 (MQ=255)
                                                                                                                                                           aTGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACggg                                                                                                                                                                                                                       >  1:167035/1‑105 (MQ=255)
                                                                                                                                                           aTGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACggg                                                                                                                                                                                                                       <  2:167035/105‑1 (MQ=255)
                                                                                                                                                                    gggtacaCCGGATAAGGCGTCCCCGCCCGCTCCGGCAGGCGGGCCCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTaa                                                            <  2:147800/245‑1 (MQ=255)
                                                                                                                                                                                  atttgatgtggtaaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTATTGATAtt                                                >  1:534026/24‑251 (MQ=18)
                                                                                                                                                                                     tgatgtggtaaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCAt                                             >  1:513249/21‑251 (MQ=18)
                                                                                                                                                                                           cccccccacccGGCCGTCGTGCGCCATGCCTGATGCGACGTTTTTCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGTGTTCACGCCGCATCCGGCATTCGTGCTATTATCAACGCATATTCAGTTTATTGTCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTTATAACGGATTAACCTCTTGATATTCATCTCCTgg                                      <  2:246382/242‑1 (MQ=255)
                                                                                                                                                                                             taaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGc                                     >  1:144981/13‑251 (MQ=21)
                                                                                                                                                                                               aacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGACGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGACAATATTCTCGCCAATAAAACTAGAGACAAAATAATAGCTGAGATCATAACCCTCGTAATAACGGCTTAACGACTTGGTATTCCTCTCCttggcg                                   >  1:81996/11‑246 (MQ=17)
                                                                                                                                                                                                 catagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGc                                     >  1:36749/9‑247 (MQ=25)
                                                                                                                                                                                                       taaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCCCGGGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTc                           >  1:352085/3‑251 (MQ=25)
                                                                                                                                                                                                       aaaggggtgcgcagaTGCCTTATGCGACGTTATGCGCGTTTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGGGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTc                           <  2:473379/237‑1 (MQ=25)
                                                                                                                                                                                                         aGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGGAGTGTAGACCGGCTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGAAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCCTAACCCTCGTGATAAAGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTcc                          >  2:187822/1‑249 (MQ=25)
                                                                                                                                                                                                          gTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTTTGATCCTAACCCTCGTGATAACGGATTAACGTCTTGATATTAATCTCCTGGCAGATCTTTTCa                          >  2:359808/1‑249 (MQ=25)
                                                                                                                                                                                                            cGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGAtt                     <  2:23646/251‑1 (MQ=255)
                                                                                                                                                                                                               gcgcagaTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTg                   <  1:323473/245‑1 (MQ=255)
                                                                                                                                                                                                               gcgcagaTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGACAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCACTAAAACTGGAAACAATATAATAACTGTGATCAATACCCTCGTGATAAAGGTTTAACGTCTTAATATTCATCTCCTGGCATATCCTTTCAAGAtttt                   >  2:212093/7‑250 (MQ=21)
                                                                                                                                                                                                                  cagaTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGACTTGATATTAATCTCCTGGCAGATCTTTTCAAGATTTgg                  >  1:32312/4‑249 (MQ=37)
                                                                                                                                                                                                                                gACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGGGGTAGGCAAAAAGCTCGCCACTAAAACAGGAGACAAAATAATAGCAGTGATCATAACCCTCGTGATAACGGATTAACGGCTTGATATTCATCTCCTGGCAGGTCTTTTAAAAATTTTGAGTCCGCAGCGGtttt  >  2:393610/1‑249 (MQ=25)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ATAGTGTTGCACAGCGCGCTCAGCGCCGAAGAGCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCT  >  NC_000913/377222‑377694

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