Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,258,5750AG48.6% ‑6.4 / 53.2 37intergenic (‑62/‑216)pth/ychHpeptidyl‑tRNA hydrolase/DUF2583 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (10/9);  new base G (8/10);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGTCGCGTTGCGGCGTATTCAGCACCGGGGTTCGCCAGGCCGACAATCAATTTAATCGTCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATAAGGGTTGTAATTGTGATCACGCCCGC  >  NC_000913/1258452‑1258703
                                                                                                                           |                                                                                                                                
tGTCGCGTTGCGGCGTATTCAGCACCGGGGTTCGCCAGGCCGACAATCAATTTAATCGTCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTAt                                                                                                                   <  2:79627/139‑1 (MQ=255)
 gTCGCGTTGCGGCGTATTCAGCACCGGGGTTCGCCAGGCCGACAATCAATTTAATCGTCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATc                                                                                                                  <  2:42458/139‑1 (MQ=255)
           ggcgTATTCAGCACCGGGGTTCGCCAGGCCGACAATCAATTTAATCGTCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGAt                                                                                                        <  1:161374/139‑1 (MQ=255)
               tATTCAGCACCGGGGTTCGCCAGGCCGACAATCAATTTAATCGTCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACg                                                                                                    <  1:32856/139‑1 (MQ=255)
                         cGGGGTTCGCCAGGCCGACAATCAATTTAATCGTCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGc                                                                                          >  2:453773/1‑139 (MQ=255)
                              ttCGCCAGGCCGACAATCAATTTAATCGTCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAg                                                                                     <  1:453773/139‑1 (MQ=255)
                                                aatTTAATCGTCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCa                                                                   >  1:110451/1‑139 (MQ=255)
                                                aatTTAATCGTCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCa                                                                   <  2:110451/139‑1 (MQ=255)
                                                          tCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGc                                                                                                                       <  2:405584/77‑1 (MQ=255)
                                                          tCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGc                                                                                                                       >  1:405584/1‑77 (MQ=255)
                                                          tCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAg                                                         >  2:353996/1‑139 (MQ=255)
                                                                                  tACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGa                                                                                                         >  2:149136/1‑67 (MQ=255)
                                                                                  tACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGa                                                                                                         <  1:149136/67‑1 (MQ=255)
                                                                                      tAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAg                                                                 <  2:543880/103‑1 (MQ=255)
                                                                                      tAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAg                                                                 >  1:543880/1‑103 (MQ=255)
                                                                                      tAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGa                                                                                                         <  1:178707/63‑1 (MQ=255)
                                                                                      tAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGa                                                                                                         >  2:178707/1‑63 (MQ=255)
                                                                                       aaCTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTg                                                                                   >  1:433114/1‑84 (MQ=255)
                                                                                       aaCTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTg                                                                                   <  2:433114/84‑1 (MQ=255)
                                                                                        aCTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAg                                                                                     >  1:382098/1‑81 (MQ=255)
                                                                                        aCTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAg                                                                                     <  2:382098/81‑1 (MQ=255)
                                                                                          tGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATAAgg                         >  2:335905/1‑139 (MQ=255)
                                                                                            gcgcGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAAc                                                                                                                          <  2:171515/40‑1 (MQ=255)
                                                                                            gcgcGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAAc                                                                                                                          >  1:171515/1‑40 (MQ=255)
                                                                                              gcgTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATAAGGGTTg                     >  2:215131/1‑139 (MQ=255)
                                                                                               cgTAGTTTACTGGTTGCGGCCTCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTg                                                                                   >  2:347641/1‑76 (MQ=255)
                                                                                               cgTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGTAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATAAGGGtt                      <  1:27644/137‑1 (MQ=255)
                                                                                               cgTAGTTTACTGGTTGCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATAAGGGtt                      >  2:27644/1‑137 (MQ=255)
                                                                                               cgTAGTTTACTGGTTGCGGCCACGCTTGACAAAAAACTGCGTAACAAATGCAGATAACGTAATAATTGCCTGAGTg                                                                                   <  1:347641/76‑1 (MQ=255)
                                                                                                 tAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAG‑CAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATAAGGGTTGTAAt                 >  1:499221/1‑139 (MQ=255)
                                                                                                           gTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATaa                           <  2:473445/120‑1 (MQ=255)
                                                                                                           gTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATaa                           >  1:473445/1‑120 (MQ=255)
                                                                                                              gCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGc                                                                                                            >  2:191010/1‑36 (MQ=255)
                                                                                                              gCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGc                                                                                                            <  1:191010/36‑1 (MQ=255)
                                                                                                              gCGGCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAAGCGTTTATTTGTAAAGTTTTGTTGAAATAAGGGTTGTAATTGTGATCACGcc     >  2:384693/1‑139 (MQ=255)
                                                                                                              gCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATAAGGGTTGTAATTGTGATCACGcc     <  1:215131/139‑1 (MQ=255)
                                                                                                                 gCCCCGCTTGGCAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATAAGGGTTGTAATTGTGATCACGCCCGc  <  2:348322/139‑1 (MQ=255)
                                                                                                                           |                                                                                                                                
TGTCGCGTTGCGGCGTATTCAGCACCGGGGTTCGCCAGGCCGACAATCAATTTAATCGTCACGTTTTTTTGTCCTGAGTGTGTACATAACTGGCGCGTAGTTTACTGGTTGCGGCCCCGCTTGACAAAAAACTGCGTATCAAATGCAGATAACGTAATAATTGCCTGAGTGGACTATTAGAAAGTCAAGGTGTTCAGGCGTTTATTTGTAAAGTTTTGTTGAAATAAGGGTTGTAATTGTGATCACGCCCGC  >  NC_000913/1258452‑1258703

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