Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2560TG43.1% 16.3 / 18.5 58A19A (GCT→GCGPP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base T (25/4);  new base G (11/11);  total (43/15)
Fisher's exact test for biased strand distribution p-value = 1.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGG  >  NC_002947/3449136‑3449383
                                                                                                                        |                                                                                                                                
aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGcc                                                                                                                  >  2:15725/1‑135 (MQ=255)
aTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCGCTGGCGGCGCCCGGGcc                                                                                                                  >  2:680005/1‑135 (MQ=255)
         gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCGGCGCCCCGGCCGCCAGGctg                                                                                                         >  1:114518/1‑136 (MQ=255)
         gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGgtg                                                                                                         >  1:670836/1‑136 (MQ=255)
         gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCACCGCTGGCGGCGCCCGGGCAGCCgggggg                                                                                                         >  1:203691/1‑134 (MQ=255)
         gCAACCGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCGCAGCGCAGGCCGCGCCCGGCCAGccacggat                                                                                                         >  1:150799/1‑130 (MQ=255)
               gtgtGGAGCACCGCCAGGCACTCAACACGCGCGCGCGGGTGGTGGTGCTGTGGGCGATTGGGCCGCCGGCGTGGTCGCCGGGGCGATGGGGGGGGCCGCCCGGGATGGGGGCGCCCGGGCAGCCTGGGTGTGGCAg                                                                                                   <  2:442084/133‑1 (MQ=255)
                 cgggaggaGCCCAAGGCCCTCAACAAGCGCGCGCTCGAGGTGCTGCTGTGCGCGCGGGGGCCACCCGCCGGGGCGCGGTGGCGCTGCGGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGc                                                                                                 <  2:1040136/134‑1 (MQ=255)
                   ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCGGCGCCCGCGCAGCCAGGGTGTGGCAGGCCa                                                                                               >  1:461477/1‑136 (MQ=255)
                    gtggcccGCACGGCCCTATACACGCGCGCGCTCGGGGGGCGGCTGTGCGCGCGGGGGCCACCCGCCGGGCCGCCGGGGCGGTGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCaa                                                                                              <  2:712214/130‑1 (MQ=255)
                        aTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTTGTCGCCGTTGCGCTGCTTCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCAGGCCAAcccc                                                                                          >  1:673688/1‑136 (MQ=255)
                        aTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAAcccc                                                                                          >  1:368557/1‑136 (MQ=255)
                            cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCgata                                                                                      >  2:302378/1‑136 (MQ=255)
                            cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGGGGCCGGCCAACGCCgata                                                                                      >  1:207639/1‑136 (MQ=255)
                             aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGgccgccgccgaaaa                                                                                     >  1:384516/1‑131 (MQ=255)
                             aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCCGGCAGCCTGGGTGTGGCAGGCAAACCCCTATaa                                                                                     >  2:8123/1‑136 (MQ=255)
                                gCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTAGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgc                                                                                  <  2:97009/136‑1 (MQ=255)
                                    tctACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAg                                                                              >  2:258268/1‑136 (MQ=255)
                                       aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTc                                                                           >  2:794372/1‑136 (MQ=255)
                                        cAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCg                                                                          >  1:1080182/1‑136 (MQ=255)
                                                cgcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCAGGTCCCCCTTCCCCAGCAGCAGACCCCCACCACAATCCGCCCCCCCGCAGCCAGGGTGTTGCAGGCCACCCCCCTAAAGGCCCACCGCGCCGATCac                                                                  >  2:139199/1‑135 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAgg                                                                 >  1:259042/1‑136 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCCGGCCCACGCCTATAAGGCGCCGCTCGCTGGTCAgg                                                                 >  1:437474/1‑136 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCAATAAGGCGCAGCTCGCTGATCAgg                                                                 >  2:1037384/1‑136 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAgg                                                                 >  1:430169/1‑136 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCCCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACCCCGATAAGGCGCAGCTCGCTGATCAgg                                                                 >  2:36045/1‑136 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCCGGCAGCCGGGGTGTGGCAGGCCAACGCCAATAAGGCGCAGCTCGCTGATCAgc                                                                 >  2:22166/1‑135 (MQ=255)
                                                 gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCGGGCGGCGCCCCGGCAGCCGGGGTGTGGGAGGCCCACGCCAATAAGGCGCCGCTCGCTGATAAgg                                                                 >  2:78718/1‑136 (MQ=255)
                                                 gcTAGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGGGCAGCCTGGCTGTGGCAGGCCACCCCCGAAAAGGAGCAGCGCGCGGATCAgc                                                                 >  1:60104/1‑135 (MQ=255)
                                                   tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACCCCGATAAGGCGCAGCTCGCTGATCAGGcc                                                               >  1:459085/1‑136 (MQ=255)
                                                   tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCCCGCAGGCGGCGCCCGGGCCCCCTGGGTGTGGCAGGGCAACGCCTATAAGGCGCCCGTCGCTGATCAGccc                                                               >  2:404086/1‑136 (MQ=255)
                                                     gAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCAGATCAGGCCgc                                                             >  2:995510/1‑135 (MQ=255)
                                                              cAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGCGGTGTGGGAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCCATGagatacg                                                    >  1:548939/1‑130 (MQ=255)
                                                               aGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGGGTGGCAGGCCCCCGCCGATAAGGCGCAGCTCGCTGATCAGGCCGCTGACTTCGGc                                                   >  1:397365/1‑136 (MQ=255)
                                                                    gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCAACGCCGATAATGCGCAGCGCGCGGATCAGGCCGATGACTTCggggagga                                              >  2:902071/1‑130 (MQ=255)
                                                                    gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCACCCGGGGTGTGGCCGGCCCACCCCGATAAGGCGCCGCTCGCGGATCAGGCCGAGGGCGTCGGCCAGCa                                              >  2:960625/1‑136 (MQ=255)
                                                                    gcgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCAGCCCAGGCGGCGCCCCGGCATCCAGTGTGTGGCAGGCCAACGCCTAAAAGGCGCAGCTCGCTGATAAGGCCGCAGCCTTCGGCCAGCa                                              >  2:855405/1‑136 (MQ=255)
                                                                      gcTTGGGCCCCCCGCGGGGGCGCGGGGGCGCGGGGGCGGGCCGCCCGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                            <  2:563235/136‑1 (MQ=255)
                                                                      gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                            >  1:22139/1‑136 (MQ=255)
                                                                           ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGCCGGCCCCCCCCTATAAGGCGCAGCTCGCTGATCAGGCCGCAGACGTCGGCCAGCAGTTGGCg                                       >  2:990896/1‑136 (MQ=255)
                                                                                  cGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCAACCCCGATAAGGCGCAGCTCGCTGATCAGGCCGCTGACGTCGGCAAGCAGTTGGCGGGGAAgg                                >  1:240347/1‑136 (MQ=255)
                                                                                         tCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCCACGCCGATAAGGCGCCGCGCGCGGATCAGGCCGCTGACTTCGGCCAGCAGGTGGGGGGGAAGGGTcgcgc                         >  2:663013/1‑136 (MQ=255)
                                                                                         gcggggTGGCGGTGCTGGGGGCCGCGAGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                         <  2:627784/131‑1 (MQ=255)
                                                                                         gcGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                         <  1:1013991/135‑1 (MQ=255)
                                                                                         gcGCCGTGGCGCGGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                         <  1:1002966/135‑1 (MQ=255)
                                                                                              gTGGCGCGGCTGCCGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACg                    <  1:351045/136‑1 (MQ=255)
                                                                                                ggCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGGGGCCCACCCCGATAAGGCGCCGCTCGCGGGGGAGGGCGCTGGCGTCGGGAAGCAGTTGGGGGGGGAGGCTCGCGCGCACggg                  >  2:917696/1‑135 (MQ=255)
                                                                                                ggCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCCGGCCGCCGGGGGGTGGCAGGCCAACCCCGATAAAGCCCCCCTCGCGGATAAGGCCCCAGACTTCGGGACACAGGTGGCGGAGAAGGATCGCGCGCCCggg                  >  2:919423/1‑135 (MQ=255)
                                                                                                ggCGCTGCTGCTGGCCGCCAGCCCTGGCGGCGCCCCGGGCCGCGGGGGGTGGGCGGGCAACGCCAAAAAGGCGCCGCTCGGGGGTGAGGCCGCGGCCTAGGGGAAGCAGTTGGCGGAGCAGGAGAGCGCGCCCgcg                  >  2:452185/1‑134 (MQ=255)
                                                                                                 gcggcgctgctgCCCGCCCGCGGGGCCGGCGCCCGGGCAGCCTGGGTGTTGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTc                 <  1:754005/131‑1 (MQ=255)
                                                                                                  cgcTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCGATAAGGCGCAGCTCGCTGGTCAGGCCGATGACTTCGGCCAGCAGTTGGGGGTGAAGGGTCGCGCGCGCggggt                >  2:907716/1‑133 (MQ=255)
                                                                                                          ctgGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCGCGCTGATCAGGCCGCTGACTTCGGCAAGCAGTTGGCGGTGAAGGATCGCGCGCCCGGTGTGGAGCggg        >  2:1104532/1‑136 (MQ=255)
                                                                                                           gggCCGCCTGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCTGga       <  2:63522/135‑1 (MQ=255)
                                                                                                            ggCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCGGATCAGGCCGATGAGTTCGGCAAGCAGTTGGCGGAGGAGGATCGCGCGCACGGTCTGGAGCGGgag      >  1:1073162/1‑136 (MQ=255)
                                                                                                              ccgcccgcGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa    <  1:663013/136‑1 (MQ=255)
                                                                                                              ccgcccgcGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa    <  1:504696/136‑1 (MQ=255)
                                                                                                              ccgccAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa    <  1:1112921/136‑1 (MQ=255)
                                                                                                                cccAGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAgg  <  2:642787/135‑1 (MQ=255)
                                                                                                                        |                                                                                                                                
ATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGG  >  NC_002947/3449136‑3449383

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: