Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913288,3820AC39.5% 78.3 / 12.2 44T532P (ACC→CCC) yagHCP4‑6 prophage; putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base A (20/6);  new base C (0/17);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 2.42e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGC‑GCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGC  >  NC_000913/288154‑288615
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         
ccggaccgggCGGAGGGCGTCGGGCTGGGGGTGGACGTGGTTACGCGGGTCTCCCGCTACAGCTACTGGTTGGGTGGGGGGAGGGGGCACACCGTGCGGGTGACGTAGGGGGGGGGGAAGGTGTCGGCGGATTACATGGGGGGGC‑GCGGCTTCTTCACCGGCGGGTTTGGGGGCCTGCACTGCGAGGACATCAGGGGCGCGGGGTGCTACGGGGATTTCGACTACTCCCCCTACGAGc                                                                                                                                                                                                                                  <  2:2034555/229‑1 (MQ=255)
 cGGAACAGGGGGGGGGCGTCGGGTTGCGGGTGGAGGTGGATCCGCTGTTCTACCGCTAAAGCTACTCGTTTGATGGCGAGAGGGGGCCCCCCGGGCCGGGGACGTAGGGGGCGGGGAAGCTGTCGGACGACTACATCGGCGGGC‑GCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCGGCGGCGACGGCTGCTACGCGGACTTCGACTCCTTCCCCTCCGAGCCg                                                                                                                                                                                                                                <  2:875057/239‑1 (MQ=255)
                  gTCTGGCTGGGGGTGGGGGTGGTTACGCTGGTCTACGGCTACAGCTACTGGTTGGATGGCGAGAGGGGGCACACCGGGCGGGTGAGGTAGGAGGCGTGGAGGCGGTGGGACGACTACATCGGGGGGC‑GCGGCTTCTTCCCCGGCGGGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGc                                                                                                                                                                                                                      <  2:2253106/232‑1 (MQ=255)
                                      ggATCCGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGCCTACATCGGCGGGC‑GCGGCTTCTTCCCCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGc                                                                                                                                                                                                                      <  1:1330647/212‑1 (MQ=255)
                                                                                                 gggTGAGTTATGGGGCGTGGAGCTTGCCGGACGACTCCA‑CGGGGGGCGGGGGCTTCTTCACCGGGGCGTTTGTGGGCCGGCCCTGGGAGGACATCGGGGGGGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCa                                                                                                                                <  2:1101052/238‑1 (MQ=255)
                                                                                                  ggTGACGTAGGGGGCGGGGAAGCTGCGGGACCACTACATCGGGGGGC‑GCGTCTTCTGCCCCGGCGGGTTGGGGGGCCTGCACTGCGAGGACACCGGCGGCGACGGCTGCTACGGGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCa                                                                                                                                <  2:285030/238‑1 (MQ=255)
                                                                                                       cGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGC‑GCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCt                                                                                                                                                                                                                     <  1:1779360/148‑1 (MQ=255)
                                                                                                                                    tACATCGGCGGGC‑GCGGCTTCTTCCCCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGCCTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCgag                                                                                                                                                                                          <  1:312946/146‑1 (MQ=255)
                                                                                                                                            cgggc‑gCGGCTTCTTCCCCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGa                                                                                         <  1:1149373/235‑1 (MQ=255)
                                                                                                                                              ggc‑gCGGCTTCTCCCCCGGCGGGTTTGGGGGCCTGCCCGGCGGGGACATCGGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCgg                                                                                    <  1:106644/238‑1 (MQ=255)
                                                                                                                                                c‑gcgGCTTCTTCCCCGGCGCTTTTGGGGGCCTGCACCGCGAGGACAGCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGa                                                                                  <  1:2296115/238‑1 (MQ=255)
                                                                                                                                                        cctccccgggCGCTTTGGGGGGCCTGCCCGGCGAGGCCACCGGCGGCGCCGGCTGCTACGCGGACTTCGCCTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCgg                                                                                    <  2:1824279/221‑1 (MQ=255)
                                                                                                                                                              ccGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGCCGACAGGCTCATGGCGGCGATGACGCGGGACTCGTGGTTCcac                                                                     >  1:1121279/1‑238 (MQ=255)
                                                                                                                                                                ggCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACAcc                                                                  >  1:7523/1‑239 (MQ=255)
                                                                                                                                                                  cgcgTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGGCGCGGGACTCGTGGTTCCACACCg                                                                 >  2:209176/1‑238 (MQ=255)
                                                                                                                                                                  cgcgTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGt                                                                          >  1:1824279/1‑229 (MQ=255)
                                                                                                                                                                    cgTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGc                                                               >  2:358733/1‑238 (MQ=255)
                                                                                                                                                                      tttGTGGGCCTGCCCGGCGAGGACATCAGCGGGGACGGCTGTTCCGCGGCCTTCGACTACTCCACCTCCGAGCCGGTCTAACGGCTCAgg                                                                                                                                                                                                                 <  2:397248/90‑1 (MQ=255)
                                                                                                                                                                                tGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGcc                                                                                                                                                                                                       >  1:397248/1‑90 (MQ=255)
                                                                                                                                                                                       cGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCCCCTACGAGCCGGTCTAACGGCTCAgg                                                                                                                                                                                                                 <  2:1064461/73‑1 (MQ=255)
                                                                                                                                                                                       cGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAgg                                                                                                                                                                                                                 >  1:1064461/1‑73 (MQ=255)
                                                                                                                                                                                       cGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGc                                            >  1:1041803/1‑238 (MQ=255)
                                                                                                                                                                                       cGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAACTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACCGGCTCAGGGCGGAGCTGACGCGGGACGCGGGGTTCCGCACCCGCACAGCCCCGCCGGGCTCGc                                            >  2:909896/1‑238 (MQ=255)
                                                                                                                                                                                         aGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGCCTCCTCCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGcc                                                                                                                                                                           <  2:1378994/109‑1 (MQ=255)
                                                                                                                                                                                            acacCAGCGGCGACGGCGGCTCCGCGGATTTCGACTACTTCCCCTACGAGCCGGTCTAACgg                                                                                                                                                                                                                       <  2:611793/62‑1 (MQ=255)
                                                                                                                                                                                                 aGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGa                                                                                                                         >  1:1497749/1‑151 (MQ=255)
                                                                                                                                                                                                 aGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCCCCACCCAGCGCACCCCCTGCTCGtt                                  >  2:1552963/1‑238 (MQ=255)
                                                                                                                                                                                                   cggcgACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAgct                                                                                                                                                                 >  1:1378994/1‑109 (MQ=255)
                                                                                                                                                                                                    ggcgACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCgg                                                                                                                                                                                                             >  1:611793/1‑64 (MQ=255)
                                                                                                                                                                                                        ggggCTGCAGCGGGGACTGAGACGTCCTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCg                           <  1:2320761/236‑1 (MQ=255)
                                                                                                                                                                                                        ccgctgcctccgcgGCCTCCGCCTCCTCCCCCTCCGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCg                           <  1:1731522/228‑1 (MQ=255)
                                                                                                                                                                                                        aCGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTc                            >  1:1405853/1‑237 (MQ=255)
                                                                                                                                                                                                          gactgctACGCGGCCTTCGACTCCTCCCCCAACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGc                                                                                                                                                 <  2:1643534/116‑1 (MQ=255)
                                                                                                                                                                                                                    cgGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCcgcg                                                                                                                                       >  1:1643534/1‑118 (MQ=255)
                                                                                                                                                                                                                      gACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGaaa                                                                                                                                                      >  1:116319/1‑101 (MQ=255)
                                                                                                                                                                                                                      gACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGaaa                                                                                                                                                      <  2:116319/101‑1 (MQ=255)
                                                                                                                                                                                                                      gACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTa             >  2:272034/1‑238 (MQ=255)
                                                                                                                                                                                                                                cctTCGCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTcgc  <  1:597739/238‑1 (MQ=255)
                                                                                                                                                                                                                                 ctTCCCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCgg                                                                                                                                                                             <  2:43065/67‑1 (MQ=255)
                                                                                                                                                                                                                                 ctTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCgg                                                                                                                                                                             >  1:43065/1‑67 (MQ=255)
                                                                                                                                                                                                                                 ctTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGa                                                                                                             >  1:1672413/1‑131 (MQ=255)
                                                                                                                                                                                                                                 ctTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGa                                                                                                             <  2:1672413/131‑1 (MQ=255)
                                                                                                                                                                                                                                   tCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGa                                                                                                             <  1:1591620/129‑1 (MQ=255)
                                                                                                                                                                                                                                   tCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGa                                                                                                             >  2:1591620/1‑129 (MQ=255)
                                                                                                                                                                                                                                     aCCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAggg                 >  1:1868813/1‑219 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         
CCGGAACATGCGGAGAGCGTCTGGCTGCGGGTGGACGTGGATACGCTGGTCTACCGCTACAGCTACTCGTTTGATGGCGAGACGTGGCACACCGTGCCGGTGACGTATGAGGCGTGGAAGCTGTCGGACGACTACATCGGCGGGC‑GCGGCTTCTTCACCGGCGCGTTTGTGGGCCTGCACTGCGAGGACATCAGCGGCGACGGCTGCTACGCGGACTTCGACTACTTCACCTACGAGCCGGTCTAACGGCTCAGGCCGGGTAGCCCAGCGCGCGCGAGAGCGCGAGCCCGGCCTGCTGAAGCTGCTCGCGGAAATTAGCCAGCTCCGCGTCGTCCACGCGGGAGGTCAGCGTCGACAGGCTCAGGGCGGCGATGACGCGGGACTCGTGGTTCCACACCGGCACCGCCACGCAGCGCACGCCCTGCTCGTTCTCTTCGCTGTCCAGGGCGTAGCCTTGCTCGC  >  NC_000913/288154‑288615

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑