Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,574,229 G→C 9.0% G301R (GGC→CGC)  glgX → glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,574,2290GC 9.0% 190.9 / 3.8 62G301R (GGC→CGC) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base C (4/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGAAAT  >  NC_002947/4574081‑4574377
                                                                                                                                                     |                                                                                                                                                     
gCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGcgcg                                                                                                                                                       <  7:198000/149‑1 (MQ=255)
gCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGcgcg                                                                                                                                                       <  3:261020/149‑1 (MQ=255)
   cccGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCa                                                                                                                                                    >  2:319842/1‑149 (MQ=255)
    ccGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCAt                                                                                                                                                   >  1:149430/1‑149 (MQ=255)
    ccGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTACAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTAGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCag                                                                                                                                                   >  2:283241/1‑148 (MQ=255)
       cgcTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGa                                                                                                                                                <  8:289608/149‑1 (MQ=255)
        gcTACCTGGCCAGCGGCAAGATCGCCGACTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGa                                                                                                                                                >  1:248587/1‑148 (MQ=255)
            ccTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCCGCGCATCGACagg                                                                                                                                            <  7:314272/148‑3 (MQ=255)
              tGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGcc                                                                                                                                         >  4:152955/1‑149 (MQ=255)
              tGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCAACGACAACGcc                                                                                                                                         >  4:215597/1‑149 (MQ=255)
                gCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGcc                                                                                                                                         <  7:27235/147‑1 (MQ=255)
                gCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTc                                                                                                                                       <  8:39369/149‑1 (MQ=255)
                gCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTc                                                                                                                                       <  2:56949/149‑1 (MQ=255)
                 ccAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCg                                                                                                                                      <  6:350494/149‑1 (MQ=255)
                 ccAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCCGCGCATCGACa                                                                                                                                              >  7:256909/1‑141 (MQ=255)
                 ccAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCCGCGCATCGACa                                                                                                                                              <  8:256909/141‑1 (MQ=255)
                   aGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGta                                                                                                                                    <  8:58442/149‑1 (MQ=255)
                      ggCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGtac                                                                                                                                   >  1:337213/1‑148 (MQ=255)
                               gCCTAGTTCAATGTGTTTGTCGCGCACCTGAAGGAGGCGGGGCTTGAGGTTATCCTGGACGTGGTTTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGGGCGGAATCTAAAACGCCTCGTACTACCGCCTGAt                                                                                                                        <  4:223921/149‑1 (MQ=255)
                               gCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGAt                                                                                                                        <  8:219245/149‑1 (MQ=255)
                                  gAGTTCAAGGAGATGGGCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGTCGTGGGGTTCAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGcc                                                                                                                     <  6:115673/149‑1 (MQ=255)
                                     ttCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGga                                                                                                                  >  3:288418/1‑149 (MQ=255)
                                      tCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGga                                                                                                                  <  7:188242/148‑1 (MQ=255)
                                             gaTGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGAAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgc                                                                                                         >  4:287929/1‑149 (MQ=255)
                                             gaTGGGCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGACATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgc                                                                                                         <  8:188157/149‑1 (MQ=255)
                                                  tCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGctac                                                                                                    <  7:22503/149‑1 (MQ=255)
                                                     cgcACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTAGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGctacta                                                                                                  >  7:94782/1‑148 (MQ=255)
                                                       atcgtGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACAt                                                                                               <  7:56061/145‑1 (MQ=255)
                                                          cTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCaa                                                                                            <  3:287929/149‑1 (MQ=255)
                                                             cACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAAc                                                                                           <  4:288418/147‑1 (MQ=255)
                                                                gacgCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTc                                                                                      >  2:246121/1‑149 (MQ=255)
                                                                gacgCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACCAGCGCGGCCCGACCCTGCGGTTGCGCGGCATCGACAACGCCTCGTTCTACGGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTc                                                                                      <  6:345417/149‑1 (MQ=255)
                                                                  cgcgGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCTCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccg                                                                                    <  5:81148/149‑1 (MQ=255)
                                                                   gcgGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccgg                                                                                   >  7:354767/1‑149 (MQ=255)
                                                                     ggggCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTg                                                                                                                          >  2:338394/1‑109 (MQ=255)
                                                                     ggggCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTg                                                                                                                          <  1:338394/109‑1 (MQ=255)
                                                                      gggATTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCTGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggcac                                                                                <  3:137935/149‑1 (MQ=255)
                                                                        gCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGCGCGCGGCCCGACCCTGTAGAAGCGCCGCATCGACAACGGCTCGTACTACCGCCAGATGCCCGACGACAAGCGCTACTACATCAACGACTccggcaccg                                                                              >  1:308650/1‑149 (MQ=255)
                                                                              aGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCaa                                                                          <  1:55361/148‑1 (MQ=255)
                                                                               ggTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGCTGCCGGACGACAAGCGCTACTACATCAACAACTCCGGCACCGGCAaca                                                                        >  5:308845/1‑149 (MQ=255)
                                                                                gTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggc                                                                                  >  7:237837/1‑138 (MQ=255)
                                                                                gTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggc                                                                                  <  8:237837/138‑1 (MQ=255)
                                                                                   aTCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCTGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCt                                                                    <  1:319842/149‑1 (MQ=255)
                                                                                   aTCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCt                                                                    >  6:12408/1‑149 (MQ=255)
                                                                                          aCGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCt                                                              >  8:385318/1‑148 (MQ=255)
                                                                                          aCGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTg                                                             <  5:12408/149‑1 (MQ=255)
                                                                                           cGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                            <  8:42175/149‑1 (MQ=255)
                                                                                                ttgtagaCCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAc                                                       <  8:347203/143‑1 (MQ=255)
                                                                                                    cAACCACACCGCCGCAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCt                                                   <  8:294143/149‑1 (MQ=255)
                                                                                                     aaCCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTg                                                  >  1:213392/1‑149 (MQ=255)
                                                                                                              gccgACGGCAACCAGCGCGGCCCGACCCGGCCCTTCCGGGGCATCTAAAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCa                                         <  2:208989/149‑1 (MQ=255)
                                                                                                               ccgAAGGCAACGAGCGCGGCCCGACCTTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAg                                        >  4:193346/1‑149 (MQ=255)
                                                                                                                cgAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGt                                       >  2:420743/1‑149 (MQ=255)
                                                                                                                cgAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGt                                       <  8:354767/149‑1 (MQ=255)
                                                                                                                 gAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGtt                                      <  7:385318/149‑1 (MQ=255)
                                                                                                                       aaCGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCAc                                >  2:17167/1‑149 (MQ=255)
                                                                                                                                   ccGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgcttcg                    >  7:226352/1‑149 (MQ=255)
                                                                                                                                    cGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgcttcgc                   >  5:125757/1‑149 (MQ=255)
                                                                                                                                       cccTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTAc                >  1:251284/1‑149 (MQ=255)
                                                                                                                                         ttGTCGATGCCGCGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTg              >  5:245622/2‑149 (MQ=255)
                                                                                                                                          tGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTgg             >  2:12881/1‑149 (MQ=255)
                                                                                                                                          tGTCGATGCCGCGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTg              >  1:326150/1‑149 (MQ=255)
                                                                                                                                            tCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGc           >  8:40748/1‑149 (MQ=255)
                                                                                                                                                tGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGc       <  3:193346/149‑1 (MQ=255)
                                                                                                                                                  cgcgGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCt                                                   <  8:92109/103‑1 (MQ=255)
                                                                                                                                                  cgcgGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCt                                                   >  7:92109/1‑103 (MQ=255)
                                                                                                                                                   gcgGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCTGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACt               <  7:241484/138‑1 (MQ=255)
                                                                                                                                                   gcgGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACt               >  8:241484/1‑138 (MQ=255)
                                                                                                                                                     ggCATCTACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGAAAt  <  2:213392/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGAAAT  >  NC_002947/4574081‑4574377

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