Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2560TG56.2% ‑4.1 / 25.2 53A19A (GCT→GCG) PP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base T (19/2);  new base G (16/11);  total (40/13)
Fisher's exact test for biased strand distribution p-value = 2.22e-02
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.

AAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAG  >  NC_002947/3449126‑3449382
                                                                                                                                  |                                                                                                                              
aaaCGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcgg                                                                                                                           >  1:114515/1‑136 (MQ=255)
                 tAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTggg                                                                                                          >  2:955915/1‑136 (MQ=255)
                   gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCAGGgtg                                                                                                        >  1:846538/1‑136 (MQ=255)
                   gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCAGGCGGCGCCCCGGCCGCCTGGgtg                                                                                                        >  2:1092473/1‑136 (MQ=255)
                   gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGgtg                                                                                                        >  1:858970/1‑136 (MQ=255)
                   gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCAGGgtg                                                                                                        >  2:932732/1‑136 (MQ=255)
                   gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGGGCAGCCTGGgtg                                                                                                        >  2:1192191/1‑136 (MQ=255)
                   gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTCCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGgtg                                                                                                        >  2:826833/1‑136 (MQ=255)
                   gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGAGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGCCCGCCCGCGCGGGCGCCGCCCGGGCACCCAggggg                                                                                                        >  1:1338416/1‑134 (MQ=255)
                         ggccggaggaGCGCCAGCCCCTCAACCCGCGCGGGGGCGAGGTGCGGCCGTGCGCGCGGGGGCCCCCCGCGGGGGCGCCGGGGCGGGGCGGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAg                                                                                                  <  2:1165835/132‑1 (MQ=255)
                          tCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCCCCCGGGCAGCCAGGGTGTGGCAgg                                                                                                 >  2:952635/1‑136 (MQ=255)
                                   ccccccggCCCTCTACCCGCGCGCGCTCGGGGGGCGGCCGGGCGCGCGGGGCCCCCCCGCCGGGGCGCGGGGGCGGGGCGGCGCCCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCt                                                                                        <  2:1054771/130‑1 (MQ=255)
                                      cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCAGGCGGCGCCCCGGCCGCCAGGGGTTGGCAGGCCCACGCCgata                                                                                     >  1:794339/1‑136 (MQ=255)
                                       aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCACCGCCTCTaa                                                                                    >  1:213284/1‑136 (MQ=255)
                                       aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCACGCGGCGCCCCCGCAGCCTGGGTGTGGCAGGCCCACGCCTATaa                                                                                    >  1:510039/1‑136 (MQ=255)
                                       aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGGAGGCCAACGCCGATaa                                                                                    >  2:16627/1‑136 (MQ=255)
                                                 aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTc                                                                          >  2:44688/1‑136 (MQ=255)
                                                        agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCACCGCCGATAAGGCGCAGCGCGCTgcta                                                                   >  1:291038/1‑133 (MQ=255)
                                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCGGATCAgg                                                                >  2:616580/1‑136 (MQ=255)
                                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGGGGCCGGCCCACCCCGAAAAGGCGCCGCGCGCGGAGCAgg                                                                >  1:1137109/1‑136 (MQ=255)
                                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCGCCCGGGCCGCCGGGGTGGGGCTGGCCAACCCCGATAAGGCGCCGCGCGCAGATGAgc                                                                >  1:1369002/1‑135 (MQ=255)
                                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCTGGCGGCGCCCCGGCAGCCGGGGTGTGGCAGGCCCACGCCAATAAGGCGCAGCTCGCGGAACAgg                                                                >  1:453151/1‑136 (MQ=255)
                                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCGCGCTGATCAgg                                                                >  1:625798/1‑136 (MQ=255)
                                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGATCAgg                                                                >  2:311452/1‑136 (MQ=255)
                                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGCGCAGCCTGGGTGTGGCAGGCCAACCCCTAAAATTCGCAGCACGCTCATCAgg                                                                >  1:72174/1‑136 (MQ=255)
                                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCCGCGCCCGGGCCGCCTGGGTGTTGAAGGCCCACCCCAATAAAGCGCAGCTCGCGGAACAgg                                                                >  1:808621/1‑136 (MQ=255)
                                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCAGGCGGCGCCCGCGCAGCCGGGGTGTGGCCGGCCCACCCCGATAAGGCGCCGCACGCAGATAAgg                                                                >  2:534982/1‑136 (MQ=255)
                                                           gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCGCGGGCGGCGCCCCGGCAGCCAGGGGGTGGCAGGCCAACGCCGAAAAAGCGCACCACGCAGAtaacg                                                                >  2:619864/1‑132 (MQ=255)
                                                             tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGGGTGGCAGGCCAACGCCTATAAGGCGCCGCTCGCGGAGCAGGcc                                                              >  1:749190/1‑136 (MQ=255)
                                                             tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCCCCTGGGTGTGGCACGCCAACGCCGATAAGGCGCAGCTCGCGGATCAGGcc                                                              >  2:1146963/1‑136 (MQ=255)
                                                                          gTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGGGGCCACCGCCGATAAGGCGCCGCGCGCGGATCAGGCCGGGGACTTgggga                                                 >  1:1327470/1‑134 (MQ=255)
                                                                              gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCCCTGGCGGCGCCCCGGCAGCCAGGGTGTGGCAGGCCAACCCCGATAAGGAGCAGCGCGCTGATCAGGCCGATGACATCGGCAAGCa                                             >  2:1085298/1‑136 (MQ=255)
                                                                                gcgggggCCCCCCGCGGGGGCGCCGGGGCGCGGCTGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                           <  1:1042242/132‑1 (MQ=255)
                                                                                gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCCGCGCGCGGCTCAGGCCGATGAATTAGGCGCGCAgg                                           >  1:791983/1‑135 (MQ=255)
                                                                                gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCcagcagg                                           >  2:974615/1‑135 (MQ=255)
                                                                                gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGGGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGAGGACTTCGGCAAGCAgg                                           >  1:123113/1‑135 (MQ=255)
                                                                                gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCAGGGGGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGATCCGGCCGCGGACGTCGGCAAGCCGt                                           >  2:998750/1‑136 (MQ=255)
                                                                                gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGGCCGGGCCGCCGCGGTGGGGCAGGCCAACGCCTAAAAGGAGCAGCGCGCGGATAACGCCGATGCCTTCGGCAAGctct                                           >  1:455690/1‑133 (MQ=255)
                                                                                gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCCACCCCGATAACGCGCAGCTCGCTGATCACGCCGAAGACGTACGCAAGCacg                                           >  2:1051086/1‑134 (MQ=255)
                                                                                  gggggCCCCCCGCCTGGGCGCCGTGGCGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTg                                         <  1:1312953/134‑1 (MQ=255)
                                                                                  gggggCCCCCCGCCGGGGCGCCGGGGCGCGGCGGCGGGCCGCCCGCGCGGGCCGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTg                                         <  1:181138/134‑1 (MQ=255)
                                                                                        cACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCCACGCCTATAAGGCGCAGCTCGCGGCTCAGGCCGCTGACTTCGGGAAGCAGTTGGGGgag                                   >  1:1322899/1‑136 (MQ=255)
                                                                                               cTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGGCCGGGCCGCCGGGGTGTGGGAGGGCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGAGTTAGGGCAGCAGTTGGCGGTGAAGGAtg                            >  2:98225/1‑135 (MQ=255)
                                                                                               cTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCACCGCCGATAAGGCGCCGCTCGCGGATCAGGCCGCTGGCTTCGGCAAGCAGTTGGCGGGGAAGGAtg                            >  2:87714/1‑135 (MQ=255)
                                                                                                  ggcGCCGGGGGGCGGCGGCGGGCCGCCAGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcg                         <  1:1092473/134‑1 (MQ=255)
                                                                                                   gggcggTGGCGGTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                        <  1:952635/131‑1 (MQ=255)
                                                                                                   gggcggTGCCGCTGCTGCGGGCCTCCTGCGCGGGGGGCGCCCGGGCAGCCTGGGTGTGGCAGGCGAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                        <  2:773935/131‑1 (MQ=255)
                                                                                                   gcggggTGGGGCTGCTGCTGGCCGCCTGGGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                        <  1:542049/131‑1 (MQ=255)
                                                                                                   gcgcgGTGGCGCTGCGGCGGGCCGCCTGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                        <  2:937377/135‑1 (MQ=255)
                                                                                                   gcgcgGTGGCGCTGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                        <  1:328861/135‑1 (MQ=255)
                                                                                                   gcGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                        <  1:1192191/135‑1 (MQ=255)
                                                                                                                        ccgcccgcGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa   <  2:1140526/136‑1 (MQ=255)
                                                                                                                         cgccCGCGCTGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATAAGGCCGATGACTTCGGCGAGCAGTTGGCGGAGCAGGATCGCGCGCCCGGTCTGGTGCGGGAGAAg  >  2:751414/1‑136 (MQ=255)
                                                                                                                                  |                                                                                                                              
AAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAG  >  NC_002947/3449126‑3449382

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