Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,305,086 C→G 6.5% intergenic (+693/+22) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,305,0860CG 6.5% 465.8 / 15.8 136intergenic (+693/+22)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (67/60);  new base G (5/4);  total (72/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

Only 100 of 138 total aligned reads displayed.

GCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGCCGCTACATCACCGTTCAGCTTGCGTCCATAAGACGGAACGATCGCTTTCAACGTAGCCTGCCATTGCGGGCTGGAAACACGATCGCCAAATACTTTTTCCAGCAGATTCAA  >  NC_000913/2304822‑2305356
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                              
     agcgTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTgg                                                                                                                                                                                                                                                                  <  2:100425/273‑1 (MQ=255)
         tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCCCAACACCGCTTTGCCGTTCCCCACCCCGTAGGCTGGAACAGTCGGTTACGCCGCATCCCACATCAAACTCCCGAGCCGGTTGCTTGATGCGACGCTGCCGCGTCTTATCATCTTACGCCACAACCGACGTACGCTGCCGGATGTGTCGTac                                                                                                                                                                                                                                                             >  4:202497/1‑274 (MQ=11)
                  gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCTTCCAACCTCTAACGCCCGAGCCGGTTGCCTGATGCGATTCTGGCCCGTCTTCTCAGCATCCCCCACATCCGGCCTACCATCCCGGGTGTGGCGTATCCttccct                                                                                                                                                                                                                                                     >  4:100424/1‑270 (MQ=14)
                       cacCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAat                                                                                                                                                                                                                                               >  2:635290/1‑275 (MQ=255)
                                    tGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCActtct                                                                                                                                                                                                                                  >  2:662257/1‑275 (MQ=255)
                                                  ggTCGGGTAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGtt                                                                                                                                                                                                                     <  2:380258/274‑1 (MQ=255)
                                                        aTAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAAc                                                                                                                                                                                                          >  1:73579/1‑279 (MQ=255)
                                                                                         cccAGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGctgcctg                                                                                                                                                                             <  2:250690/275‑1 (MQ=255)
                                                                                              gCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGc                                                                                                                                                                         <  2:545258/274‑1 (MQ=255)
                                                                                                                          tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGcc                                                                                                                                                                                                                                                                          >  6:859290/1‑149 (MQ=255)
                                                                                                                             tATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGa                                                                                                                                                                                                                                                                       <  7:662836/149‑1 (MQ=255)
                                                                                                                             tATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGa                                                                                                                                                                                                                                                                       <  6:1536038/149‑1 (MQ=255)
                                                                                                                              aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                                                                                                                                                      <  5:964875/149‑1 (MQ=255)
                                                                                                                              aTCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCCGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                                                                                                                                                      <  5:1280314/149‑1 (MQ=255)
                                                                                                                                cAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtgt                                                                                                                                                                                                                                                                    <  6:1039594/149‑1 (MQ=255)
                                                                                                                                  ggCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTgg                                                                                                                                                                                                                                                                  <  5:859290/149‑1 (MQ=255)
                                                                                                                                  ggCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTgg                                                                                                                                                                                                                                                                  <  6:691837/149‑1 (MQ=255)
                                                                                                                                         acCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGa                                                                                                                                                                                                                                                                       <  5:399581/137‑1 (MQ=255)
                                                                                                                                         acCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATGGTATGCCGGa                                                                                                                                                                                                                                                                       >  6:399581/1‑137 (MQ=255)
                                                                                                                                             cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTAc                                                                                                                                                                                                                                                       >  5:1166618/1‑149 (MQ=255)
                                                                                                                                             cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTCCTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTcgcg                                                                                                                          >  2:335241/1‑274 (MQ=255)
                                                                                                                                             cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTa                                                                                                                                                                 >  2:359899/1‑235 (MQ=255)
                                                                                                                                                tgAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGc                                                                                                                                                    >  2:299331/1‑245 (MQ=255)
                                                                                                                                                    gTGCTCCACCCCGTAGGTTGGAGAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACtt                                                                                                                                                <  4:299328/245‑1 (MQ=255)
                                                                                                                                                     gtctcccaccctgaGGGTGGTTAAGGGGTGTACGCCGTATCCCAACTCTAAACCCCGGGCCGTTTGCCTGATGGAGCGCTGGAGCGTCTTTTCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGc                                                                                                                 <  4:662252/262‑1 (MQ=255)
                                                                                                                                                      gCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGc                                                                                                                 <  2:585511/274‑1 (MQ=255)
                                                                                                                                                       cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATTCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTTCAACGCAATATCCGCCACTTCTTTTTTAACGGGTTGCGGCTTCAACTCCGGTTTCGGCGTACTATCTGCGGCCTGCGGCCTTTCGTTGTTCAACCCCAGCACTCTGCTGGTTTGCTTCAACCCGCGTTCTTTGcca                                                                                                               >  4:585507/1‑273 (MQ=37)
                                                                                                                                                       cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCCCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGCTGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTTCGGTTTCAACTGCTGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTCGTTCAACCCCAGTACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGc                                                                                                                 >  2:520305/1‑273 (MQ=255)
                                                                                                                                                        tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATc                                                                                                                                                                                                                                            >  8:573642/1‑149 (MQ=255)
                                                                                                                                                          cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                                                                                                                                                      >  2:771001/1‑121 (MQ=255)
                                                                                                                                                          cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                                                                                                                                                      >  4:697988/1‑121 (MQ=255)
                                                                                                                                                          cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCACCCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGCCGc                                                                                                              >  2:245017/1‑273 (MQ=255)
                                                                                                                                                           aCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGc                                                                                                                                                                                                                                         <  6:1166618/149‑1 (MQ=255)
                                                                                                                                                              ccGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTgg                                                                                                                                                                                                                                                                  <  2:697993/121‑1 (MQ=255)
                                                                                                                                                              ccGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTgg                                                                                                                                                                                                                                                                  <  4:770995/121‑1 (MQ=255)
                                                                                                                                                               cGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAtgtg                                                                                                                                                                                                                                                                   >  2:345763/1‑119 (MQ=255)
                                                                                                                                                                gTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCActt                                                                                                                                                                                                                                    >  6:1716985/1‑149 (MQ=255)
                                                                                                                                                                   ggTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGt                                                                                                                                                                                                                                                               <  4:345759/119‑1 (MQ=255)
                                                                                                                                                                     ttGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCt                                                                                                                                                                        <  2:677269/204‑1 (MQ=255)
                                                                                                                                                                      tGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGc                                                                                                                                                                                                                                         >  4:373497/1‑138 (MQ=255)
                                                                                                                                                                        gATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGc                                                                                                                                                                                                                                                                 >  6:541643/1‑112 (MQ=255)
                                                                                                                                                                        gATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGc                                                                                                                                                                                                                                                                 <  5:541643/112‑1 (MQ=255)
                                                                                                                                                                         aTAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGaa                                                                                                                                                                                                                           >  7:808736/1‑149 (MQ=255)
                                                                                                                                                                           aaGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAata                                                                                                                                                                                                                                              >  1:103235/1‑128 (MQ=255)
                                                                                                                                                                           aaGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAata                                                                                                                                                                                                                                              <  3:103235/128‑1 (MQ=255)
                                                                                                                                                                             ggCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACggg                                                                                                                                                                                                                       >  6:1651964/1‑149 (MQ=255)
                                                                                                                                                                                 gTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGc                                                                                                                                                                                                                   >  6:1533388/1‑149 (MQ=255)
                                                                                                                                                                                      gccgcATCCAACATCTAACGCCCGAGCCGGTTGACTGATGCGACGCTGGCGCGTCTTGTCAGCATACGCCACATCAGTCATACCATGCCGGATGTGGCGTATCCTTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGctt                                                                                                                                                                                                              >  5:387867/1‑149 (MQ=38)
                                                                                                                                                                                       ccgcATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGCCGCTACATCACCGTTCAGCTTGCGTCCATAAGa                                                                                <  2:202500/274‑1 (MQ=255)
                                                                                                                                                                                              cAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGAt                                                                                                                                                                                                                                                                      <  7:363972/85‑1 (MQ=255)
                                                                                                                                                                                              cAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCg                                                                                                                                                                                                      >  7:55517/1‑149 (MQ=255)
                                                                                                                                                                                                  aTCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGa                                                                                                                                                                                                                                                                       <  8:720926/80‑1 (MQ=255)
                                                                                                                                                                                                  aTCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGa                                                                                                                                                                                                                                                                       >  7:720926/1‑80 (MQ=255)
                                                                                                                                                                                                  aTCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGttt                                                                                                                                                                                                  >  5:784883/1‑149 (MQ=255)
                                                                                                                                                                                                       aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTgg                                                                                                                                                                                                                                                                  <  1:257938/80‑1 (MQ=255)
                                                                                                                                                                                                       aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTgg                                                                                                                                                                                                                                                                  >  3:257938/1‑80 (MQ=255)
                                                                                                                                                                                                       aCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCg                                                                                                                                                                                             >  7:1532790/1‑149 (MQ=255)
                                                                                                                                                                                                         gCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTg                                                                                                                                           <  2:114453/197‑1 (MQ=255)
                                                                                                                                                                                                         ccccGAGCCGGTGGCCAGATCCGACGCTGGCGCGACTTATCAGCATACGCCACATCCGGCATACCATGCCTGTGGTGGCGTATCATTACAACGCAATATCCGCCACGTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGCCGCTACATCACCGTTCAGCTTGCGTCCATAAGACGGAACGATCGCTTTCAAc                                                             <  4:258829/274‑1 (MQ=255)
                                                                                                                                                                                                          cccGGGCCGGTTGCCGGATGCGACGCGGGCGCGTCTTATCAGCATACGCCCCATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTGCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGCCGCTACATCACCGTTCAGCTTGCGTCCATAAGACGGAACGATCGCTTTCAACg                                                            <  4:635285/275‑1 (MQ=255)
                                                                                                                                                                                                             gAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTAt                                                                                                                                                                                       >  8:536017/1‑149 (MQ=255)
                                                                                                                                                                                                                ccGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATctg                                                                                                                                                                                    <  4:785380/149‑1 (MQ=255)
                                                                                                                                                                                                                 cGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATctgc                                                                                                                                                                                   <  5:1716985/149‑1 (MQ=255)
                                                                                                                                                                                                                 cGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATctgc                                                                                                                                                                                   <  5:1651964/149‑1 (MQ=255)
                                                                                                                                                                                                                     tGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGctgcc                                                                                                                                                                               <  7:1422907/149‑1 (MQ=255)
                                                                                                                                                                                                                       ccTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGctgcctg                                                                                                                                                                             >  6:578175/1‑149 (MQ=255)
                                                                                                                                                                                                                          gATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCgg                                                                                                                                                                          <  8:808736/149‑1 (MQ=255)
                                                                                                                                                                                                                             gCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCtt                                                                                                                                                                       <  7:536017/149‑1 (MQ=255)
                                                                                                                                                                                                                               gACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGt                                                                                                                                                                     <  8:188793/149‑1 (MQ=255)
                                                                                                                                                                                                                                   cTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTa                                                                                                                                                                 <  5:876337/149‑1 (MQ=255)
                                                                                                                                                                                                                                        gcgTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCa                                                                                                                                                            >  7:1018205/1‑149 (MQ=255)
                                                                                                                                                                                                                                             ttATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAAcccc                                                                                                                                                       >  5:993734/1‑149 (MQ=255)
                                                                                                                                                                                                                                              tATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCa                                                                                                                                                      <  5:1533388/149‑1 (MQ=255)
                                                                                                                                                                                                                                                tCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGc                                                                                                                                                    >  8:90944/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     aTACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTc                                                                                                                                               >  8:1121888/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      tACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCg                                                                                                                                              >  8:648801/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      tACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGc                                                                                                                                             <  1:229629/150‑1 (MQ=255)
                                                                                                                                                                                                                                                      tACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGc                                                                                                                 <  1:139459/178‑1 (MQ=255)
                                                                                                                                                                                                                                                      tACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGc                                                                                                                 >  3:139459/1‑178 (MQ=255)
                                                                                                                                                                                                                                                      tACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACGTCGc                                                                                                                                             >  3:229629/1‑150 (MQ=255)
                                                                                                                                                                                                                                                          ccACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAAcc                                                                                                                                                         >  6:908769/1‑134 (MQ=255)
                                                                                                                                                                                                                                                          ccACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTgg                                                                                                                                          >  5:935963/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           cacaTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATctgc                                                                                                                                                                                   >  5:679675/1‑107 (MQ=255)
                                                                                                                                                                                                                                                            acaTCCGGCATGGTATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCACCCCCAGCa                                                                                                                                                   >  8:1148814/1‑138 (MQ=255)
                                                                                                                                                                                                                                                            acaTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTAt                                                                                                                                                                                       <  7:969310/102‑1 (MQ=255)
                                                                                                                                                                                                                                                            acaTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTc                                                                                                                                                                    <  6:1389004/121‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTCCGGCATGGTATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCttt                                                                                                                                                                                                                                >  6:800140/1‑59 (MQ=21)
                                                                                                                                                                                                                                                              aTCCGGCATGGTATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCttt                                                                                                                                                                                                                                <  5:800140/59‑1 (MQ=21)
                                                                                                                                                                                                                                                              aTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAg                                                                                                                                                     <  1:143829/134‑1 (MQ=255)
                                                                                                                                                                                                                                                              aTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTa                                                                                                                                      >  8:313766/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              aTCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGCCGCTACATCACCGTTCAGCTTGCGTCCATAAGACGGAACGATCGCTTTCAACGTAGCCTGCCATTGCGGGCTGGAAACACGATCGCCAAATACTTTTTCcagcag        >  2:386607/1‑275 (MQ=255)
                                                                                                                                                                                                                                                               tCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTAc                                                                                                                                     >  7:652543/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               tCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTAc                                                                                                                                     >  8:729117/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               tCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTAc                                                                                                                                     >  5:731156/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                ccGGCATGGTATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACt                                                                                                                                    >  5:669939/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  ggCATACCATGCCGGATGTGGCATATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGCCGCTACATCACCGTTCAGCTTGc                                                                                          <  2:611374/189‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCaa                                                                                                                                <  8:1552734/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGCCGCTACATCACCGTTCAGCTTGCGTCCATAAGACGGAACGATCGCTTTCAACGTAGCCTGCCATTGCGGGCTGGAAACACGATCGCCAAATACTTTTTCCAGCAGATTCaa  <  2:315022/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aCCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTc                                                                                                                             >  8:1145623/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        ccATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTcg                                                                                                                            >  6:1609990/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                              
GCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGGCATACCATGCCGGATGTGGCGTATCATTACAACGCAATATCCGCCACTTCTTTTTGAACGGGTTGCGGTTTCAACTGCGGTTTCGGCGTACTATCTGCTGCCTGCGGCTTGTCGTAGTTCAACCCCAGCACTTCGCTGGTGTACTGCAACTCGCGTTCTGTTGCCGCTACATCACCGTTCAGCTTGCGTCCATAAGACGGAACGATCGCTTTCAACGTAGCCTGCCATTGCGGGCTGGAAACACGATCGCCAAATACTTTTTCCAGCAGATTCAA  >  NC_000913/2304822‑2305356

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: