Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,574,228 C→G 9.3% R300R (CGC→CGG) ‡ glgX → glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,574,2280CG 9.3% 183.4 / 7.0 62R300R (CGC→CGG) ‡glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base G (3/3);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CGCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGC  >  NC_002947/4574080‑4574372
                                                                                                                                                     |                                                                                                                                                 
cGCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATgcgc                                                                                                                                                   <  3:58783/149‑1 (MQ=255)
cGCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATgcgc                                                                                                                                                   <  3:411674/149‑1 (MQ=255)
 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                                                                             <  8:237837/138‑1 (MQ=255)
                                                                                 gTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggc                                                                             >  7:237837/1‑138 (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)
                                                                                                                                                     |                                                                                                                                                 
CGCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGC  >  NC_002947/4574080‑4574372

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: