Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2560TG46.0% 6.4 / 14.3 38A19A (GCT→GCG) PP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base T (18/2);  new base G (8/9);  total (27/11)
Fisher's exact test for biased strand distribution p-value = 9.74e-03
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.

GTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  NC_002947/3449139‑3449381
                                                                                                                     |                                                                                                                             
gtgctAGCATGGGGCTGGAGGACCGCATGGCCCTCTACCCCGGCGCGCTCGTGGTGCTGCTGTGAGCGGTTGCGCCACCCGCGTGGGCGGGGGGGCGGTGGTGCTGGGCGCCCGCGCGGGCGGCGCCCGGGCAGcc                                                                                                             <  2:853342/136‑1 (MQ=255)
      gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGgtg                                                                                                       >  1:790111/1‑136 (MQ=255)
      gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGgtg                                                                                                       >  1:428536/1‑136 (MQ=255)
                ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCCGGCAGCCGGGGTGTGGCAGGccc                                                                                             >  2:591249/1‑135 (MQ=255)
                     aTCGCAACGCCCTCTACACGCGCGCGCGCGCGGTGCGGCCGGGCGCGCGGGGGCCCCCCGCCGGGGCGCCGGGGCGCTGCTGCGGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGcc                                                                                        <  1:809/136‑1 (MQ=255)
                         cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCaata                                                                                    >  1:662294/1‑136 (MQ=255)
                         cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGCAGGCCAACGCCgatt                                                                                    >  1:615387/1‑135 (MQ=255)
                         cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCCGCGCCCGGCCCCCCTGGGGGTGGCAGGCCAACCCCgata                                                                                    >  1:331218/1‑136 (MQ=255)
                         cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCACCCCCgata                                                                                    >  1:959862/1‑136 (MQ=255)
                          aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCAGCGCTGGCGGCGCCCGGGCGGCCGGGGTGTGGCAGGCCAACGCCTaaaa                                                                                   >  1:436000/1‑136 (MQ=255)
                          aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCCGCCGGGGTGTGGCAGGCCAACGccgaaaa                                                                                   >  2:623053/1‑131 (MQ=255)
                          aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCGCCCCGGCCCCCTCGGTGTGGCAGGCCACCGCCAATaa                                                                                   >  1:55233/1‑136 (MQ=255)
                          aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCTATaa                                                                                   >  2:785645/1‑136 (MQ=255)
                            ggCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCTATTAgg                                                                                 >  2:843111/1‑136 (MQ=255)
                               ccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcgc                                                                              >  1:56870/1‑136 (MQ=255)
                                  ctACAAGCGAGCGCGGGAGGTGGTGCTGGGAGCGCTGGGGCCAGCCGGCGGGGCGCGGGGGCGCTGCGGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGc                                                                           <  2:902976/136‑1 (MQ=255)
                                    aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGGAGGCCCACGCCGATAAGGCGCCGCTc                                                                         >  2:574596/1‑136 (MQ=255)
                                       cgcgagcgCTCGCGGTGCGGCAGTGAGCGGTTGGGCCACCCGCCTGGGCGGGGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                      <  1:246105/135‑1 (MQ=255)
                                              gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCCCCGCCGTGTGGCAGGCCCACGCCTATAAGGCCCAGCTCGCGGATCacg                                                               >  2:960771/1‑134 (MQ=255)
                                              gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGCCCGCCAGGGTGTGGCAGGCCAACCCCTAAAAGGCGCACCTCGCTGATCAgg                                                               >  1:956381/1‑136 (MQ=255)
                                              gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCCGCGCCCCCGCCGCCAGGGTGTGGCACGCCCACGCCAAAAAGGCGCAGCACGCTGATCAgg                                                               >  2:465685/1‑136 (MQ=255)
                                                tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGcc                                                             >  1:902976/1‑136 (MQ=255)
                                                tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCCACGCCTATAAGGCGCAGCTCGCAGGTCAGGcc                                                             >  2:60230/1‑136 (MQ=255)
                                                  gtggTGCTGTCGTGCGGGCGTGGCCCCCCCGCCGGGGCGGGGTGGCGCTGCGGGTGGCCGCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGa                                                           <  1:364159/134‑1 (MQ=255)
                                                                   gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCACGCCGATAAGGCGCCGCTCGCGGATCAGGCCGAGGACGTCGGCCAGCAGt                                          >  1:767582/1‑136 (MQ=255)
                                                                   gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGAAGACTTCGGCAAGCAGt                                          >  1:32433/1‑136 (MQ=255)
                                                                   gcTGGGGCCCCCCGGCGGGGCGGGGTGGCGCTGCTGCGGCCCGCCCGGGCGGGGGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                          <  2:15285/136‑1 (MQ=255)
                                                                     ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCCGGCCCACGCCTATAAGGCGCAGCTCGCGGATCAGGGCGATGACTTCGGCCAGCAGGTg                                        >  2:577812/1‑136 (MQ=255)
                                                                        ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTAGGCAAGCAGTTGGCg                                     >  2:250511/1‑136 (MQ=255)
                                                                        ggCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCTATAATGCGCAGCTCGCTGATCAGGCCGATGACTTCGGGAAGCAGTTGGCg                                     >  2:255464/1‑136 (MQ=255)
                                                                                gCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCTGGCGGCGCCCGGGCAGCCAGGGGGTGGCGGGCCCACGCCGATAAGGCGCAGCAAGCGGGTCAGGCCGCTGACTTCGGCCAGGAGTTGGCGGGGGAggc                             >  2:63555/1‑135 (MQ=255)
                                                                                    ggTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCCACGCCGATAAGGCGCCGCGCGCGGATCAGGCCGATGACTTCGGCGAGCAGGTGGCGGAGAAGGCTcgc                         >  1:260378/1‑136 (MQ=255)
                                                                                      gcggggTGGCGCGGCGGCGGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  2:634805/131‑1 (MQ=255)
                                                                                      gcgcgGTGGCGGGGCGGCTGGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  2:203327/135‑1 (MQ=255)
                                                                                      gcGCCGTGGCGCTGCTGCTGGCCGCCTGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  1:716669/135‑1 (MQ=255)
                                                                                              gcgcTGCGGCTGGCCGCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTc               <  1:499971/136‑1 (MQ=255)
                                                                                                       ctgGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCGGATCAGGCCGCTGACTTCGGGAAGCAGTTGGGGGAGAAGGATCGCGCGCCCGGTCTGGTGCggg      >  2:222666/1‑136 (MQ=255)
                                                                                                           ccgcccgcGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa  <  1:299173/136‑1 (MQ=255)
                                                                                                                     |                                                                                                                             
GTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  NC_002947/3449139‑3449381

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: ‑