Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2710AC50.0% 28.1 / 14.8 41A24A (GCA→GCC) PP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base A (8/12);  new base C (20/0);  total (29/12)
Fisher's exact test for biased strand distribution p-value = 4.51e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGAT  >  NC_002947/3449139‑3449403
                                                                                                                                    |                                                                                                                                    
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)
                                                                                                                               cGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCg    >  1:68568/1‑136 (MQ=255)
                                                                                                                                gggCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGa   >  2:347476/1‑136 (MQ=255)
                                                                                                                                 ggCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGAt  <  2:331218/136‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                    
GTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGAT  >  NC_002947/3449139‑3449403

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