Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,449,2510AC55.4% 14.6 / 18.7 60S18R (AGC→CGC) PP_5556bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base A (24/1);  new base C (25/6);  total (49/11)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  NC_002947/3449116‑3449381
                                                                                                                                       |                                                                                                                                  
ggCGATCACCACACGTATCAATGGTGGTAGCATGGTTGTGGAGGATCGAAAGGCCATCTACAAGCGAGCGCTCTTGGTGCTGCAGTGAGCGTTTGGGCAACCCGCTTGGTCGCCGTGGCGCTGCTGCTGGCCGCCa                                                                                                                                    <  2:1862465/136‑1 (MQ=255)
  cGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAgc                                                                                                                                  >  2:3186976/1‑136 (MQ=255)
   gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGccgcccgcc                                                                                                                                 >  2:305816/1‑135 (MQ=255)
   gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAgcg                                                                                                                                 >  1:2648848/1‑136 (MQ=255)
   gATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAgcg                                                                                                                                 >  1:409013/1‑136 (MQ=255)
         cAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTggcg                                                                                                                           >  1:1931202/1‑136 (MQ=255)
              gTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGcc                                                                                                                      >  2:484885/1‑136 (MQ=255)
                   aaTGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGc                                                                                                                 >  1:1303985/1‑136 (MQ=255)
                   aaTGGTGGTAGCAACGGTCTGGAGGATCGAAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAACCCTGGCGGCGCGCCGGc                                                                                                                 >  1:2351770/1‑136 (MQ=255)
                      ggtggtAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGc                                                                                                              >  1:2649798/1‑136 (MQ=255)
                             gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCTGGCGGCGGCCGCGCAGCCTggggg                                                                                                       >  1:2679755/1‑134 (MQ=255)
                             gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTgggag                                                                                                       >  1:2975419/1‑134 (MQ=255)
                             gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGgtg                                                                                                       >  1:678420/1‑136 (MQ=255)
                             gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCAGGCCGCGCCCGGGCCGCCTGGgtg                                                                                                       >  2:653773/1‑136 (MQ=255)
                             gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCCGGCCGCCGGGgtg                                                                                                       >  1:208214/1‑136 (MQ=255)
                             gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCCGGCAGCCTGGgtg                                                                                                       >  1:840030/1‑136 (MQ=255)
                             gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCCGGCAGCCTGGgtg                                                                                                       >  1:964320/1‑136 (MQ=255)
                                         aggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGGCCCCGCACGCGGCGCCCCCGCCGCCAGCGTGTGGCAGGCCCAc                                                                                           >  1:1124861/1‑136 (MQ=255)
                                            cgggCAGGGCCCTTTACTCGCGCGCGCGGGTGTTGCGGCTGTGCGGGGTTGGCCCACCCGCCGGGCCGGGGTGGGGGTGCTGCGGCCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGcc                                                                                        <  2:110233/133‑1 (MQ=255)
                                                cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCCCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGGCAAAGCCtgta                                                                                    >  1:1951925/1‑136 (MQ=255)
                                                 aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGGCCGCCCGGGCGGCGCGCGGGCCGCCGGGGGGGGGGAGGCCAACCCCCATaa                                                                                   >  1:429987/1‑136 (MQ=255)
                                                 aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATaa                                                                                   >  1:2080962/1‑136 (MQ=255)
                                                    cccccCTACCCGCGCGCGCGCGCGGGGCTGCCGGGCGGGCGGGGCCCCCCCGCCGGGGCGCGGGGGCGCGGCGGCGGCCCCCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgc                                                                                <  2:3074879/135‑1 (MQ=255)
                                                     cccTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCGATAAGgcg                                                                               >  2:2659087/1‑136 (MQ=255)
                                                     cccTATCCACGCGAGCGCGCGCGGTTCTGGAGTGCGGGCTTGGCCCACCCGCCTGGGCGGGGTGGCGGTGGGGGTGGCCGCCTGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGgcg                                                                               <  2:107526/136‑1 (MQ=255)
                                                         ctccccgcgtgcgCGCGTGGGGCTGCCGTGTGCGCGGGGGCCCCCCGCCGGGTCGCGGTGGCGCGGCGGCGGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGc                                                                           <  1:453704/130‑1 (MQ=255)
                                                           aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCCCGCGGGCGGCGCCCGGGCCGCCGGGGGGGGGCCGGCCCACCCCGAAAAGGcgccgcgc                                                                         >  2:2654950/1‑134 (MQ=255)
                                                           aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCCACCCCGATAAGGCGCAGcac                                                                         >  2:2567882/1‑134 (MQ=255)
                                                           aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCCACGCCTATAAGGCGCAGCTc                                                                         >  1:3119474/1‑136 (MQ=255)
                                                              cgcgcgcgcTCGTGGTGCGGCCGGGCGCGGGGGGGCCCCGCGCGTGGGCGCGGGGGGGCGGCGGCGGGCCGCCGGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCt                                                                      <  1:1526731/135‑1 (MQ=255)
                                                                  agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCGCGCGGAta                                                                  >  2:158107/1‑135 (MQ=255)
                                                                  agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCCGGCAGCCGGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCGGATc                                                                  >  1:2789719/1‑136 (MQ=255)
                                                                  agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCACCGCCGATAAGGCGCCGCTCGCTGATc                                                                  >  1:2441937/1‑136 (MQ=255)
                                                                  agcgCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCCCCCGCGCGGGCGGCCCCCCGGCCGCCACGGGGTGGCAGGCCCCCCCCAAGAAAGCGCAGCGAGCTGagc                                                                  >  2:3085978/1‑134 (MQ=255)
                                                                     gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCGCTGGCGGGGCCCGGGCCGCCTGGGTGTTGCAGGCCTACGCCTTTATGTCGCAGCTCGCGGATCAgg                                                               >  1:1181616/1‑136 (MQ=255)
                                                                     gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGCAGGCCACCGCCGATAAGGCGCAGCTCGCGGATCAgg                                                               >  1:1213854/1‑136 (MQ=255)
                                                                     gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCGGGCGGCGCCCGGGCCGCCGGGGGGGGGCAGGCCCACGCCGATAAGGCGCCGCGCGCGGCTCAgg                                                               >  2:2545392/1‑136 (MQ=255)
                                                                     gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGGGCAGCCTGGGTGGGGCAGGCCCACGCCTATAAGGCGCAGCTCGCGGATGAgg                                                               >  2:1607692/1‑136 (MQ=255)
                                                                     gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCAGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACCCCAATAAGGCGCACCACGCAGATCAgg                                                               >  1:2347012/1‑136 (MQ=255)
                                                                       tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGcc                                                             >  1:2244157/1‑136 (MQ=255)
                                                                       tCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACCCCGAAAAAGCGCAGCTCGCGGCTCAGGcc                                                             >  1:1230840/1‑136 (MQ=255)
                                                                             ggctggAGTGCGCGCGTGCGCCCCCCGCGTGGGCGCGTGGGCGCGGCGGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGAc                                                       <  2:2057893/135‑1 (MQ=255)
                                                                                 gcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCACCGCAGGCGGCGCCCGGGCAGCCTGGGTGTGGCCGGCCAACCCCAAAAAGGCGCACCTCGCGGATCAGGCCGAAGACATCg                                                   >  2:585675/1‑136 (MQ=255)
                                                                                 gcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCCACCCCGATAAAGCGCAGCGCGCGGATGAGGCCGCTGACttag                                                   >  2:503227/1‑134 (MQ=255)
                                                                                         cgcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAg                                           >  1:2944583/1‑136 (MQ=255)
                                                                                          gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGGCCGGGCCGCCGGGGTGTGGCCGGCCCACGCCTATAAGGCGCCGCTCGCGGATCAGGCCGATGACATCGGCAAGCAGt                                          >  1:1791165/1‑136 (MQ=255)
                                                                                          gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                          >  1:3051216/1‑136 (MQ=255)
                                                                                          gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGGGTGGCCGGCCACCGCCGATAAGGCGCCGCGCGCGGATCAGGCCGCGGACTTCGGCAAGCAgg                                          >  1:1554613/1‑135 (MQ=255)
                                                                                          gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCCCGCCGATAAGGCGCAGCGCGCTGATAAGGCCGCTGACGTCGGCAAGCAGt                                          >  2:2595836/1‑136 (MQ=255)
                                                                                          gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCGGGCCCCCGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGt                                          >  1:821804/1‑136 (MQ=255)
                                                                                          gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCCACGCCTATAAGGCGCAGCGCGCGGATCAGGCCGCTGACTTCGGCAAGCAGt                                          >  1:1841458/1‑136 (MQ=255)
                                                                                            ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCAGGCGGCGCCCCGGCCGCCTGGGTGTGGGAGGCCCACGCCGAGAAGGCGCAGCTCGCTGATCAGGCCGCTGACTTCGGCAAGCAGTTg                                        >  1:2536868/1‑136 (MQ=255)
                                                                                            ttGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCAGATCAGGCCGATGACATCGGCAAGCAGttt                                        >  2:1966611/1‑135 (MQ=255)
                                                                                                 ccACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCTGGGGGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGCTGACTTCGGCAAGCAGTTGGCggg                                   >  1:729681/1‑135 (MQ=255)
                                                                                                      cGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAgg                              >  2:2818265/1‑136 (MQ=255)
                                                                                                             ggggggTGGGGCTGGTGCGGGCCGCCGGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  2:1053834/131‑1 (MQ=255)
                                                                                                             gggCCGGGGGGCGGCGGCGGGCCCCCTGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  2:1155900/134‑1 (MQ=255)
                                                                                                             gcggggTGGCGCGGGGGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  1:397954/131‑1 (MQ=255)
                                                                                                             gcgcgGTGGCGCGGCGGCGGGCCGCCCGCGGGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc                       <  1:508275/135‑1 (MQ=255)
                                                                                                                      cgcTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCTGGGTGTGGCAGGCCAACGCCGATAAGGCGCCGCTCGCTGATCAGGCCGCTGACTTCGGGAAGCAGGTGGCGGAGCAGGATCGCGCGCGCgggcg              >  1:1031221/1‑133 (MQ=255)
                                                                                                                      cgcTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCGGGCCCACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACgggcg              >  1:2749744/1‑133 (MQ=255)
                                                                                                                             gcgggCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCgg       <  2:2648848/133‑1 (MQ=255)
                                                                                                                                  ccgcccgcGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGaa  <  1:158107/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
GGCGATCACCAAACGTATCAATGGTGGTAGCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAA  >  NC_002947/3449116‑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: ‑