Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2510AC55.6% 1.3 / 10.0 37S18R (AGC→CGC) PP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base A (13/3);  new base C (14/6);  total (27/10)
Fisher's exact test for biased strand distribution p-value = 7.00e-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.

CGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  NC_002947/3449118‑3449381
                                                                                                                                     |                                                                                                                                  
cGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGccgcccgc                                                                                                                                  >  1:334986/1‑136 (MQ=255)
                     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)
                                                                                                                                     |                                                                                                                                  
CGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  NC_002947/3449118‑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: ‑