Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,579,221 C→G 5.4% 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,579,2210CG 5.4% 849.5 / 16.6 267R300R (CGC→CGGglgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base C (139/109);  new base G (4/13);  total (144/123)
Fisher's exact test for biased strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

CCGCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGAAAT  >  NC_002947/4579072‑4579370
                                                                                                                                                     |                                                                                                                                                     
ccGCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACCAGCGCGGCCCGACCCTGTCGATgcgc                                                                                                                                                       <  2:194487/150‑1 (MQ=255)
      ccGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATc                                                                                                                                                 <  1:3297774/150‑1 (MQ=255)
       cgcgCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGcgcgctcgtt                                                                                                                                                <  2:3400999/150‑7 (MQ=255)
        gcgcTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGa                                                                                                                                               >  1:3376703/1‑150 (MQ=255)
        gcgcTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCACACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGa                                                                                                                                               <  1:3508951/150‑1 (MQ=255)
          gcTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGTCCGACCCTGTCGATGCGCGGCATCGACa                                                                                                                                             >  2:7176675/1‑150 (MQ=255)
           cTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCCGCGCATCGACag                                                                                                                                            <  2:4951253/150‑2 (MQ=255)
           cTACCTGGCAAGCGGAAATATCGCCTAGTTCAAGGAGTTGGTCGCGCACCTGAACGACGGGGGGATTTTGGTTATCCTGGACGTGGTCTACAACCACACCGCCTAAGGAAACGAGCGGGGCCCGACCCTGTGTATGCCGCGCATCGACag                                                                                                                                            <  1:6128647/150‑2 (MQ=255)
             aCCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACg                                                                                                                                          >  2:4160332/1‑150 (MQ=255)
               cTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGcc                                                                                                                                        >  2:285010/1‑150 (MQ=255)
               cTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGcc                                                                                                                                        <  1:6543842/150‑1 (MQ=255)
                 ggCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGTTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGAAACGAGCGCGGCCCGACCCTTTCTATGCGCGGCATCGACAACGCCTc                                                                                                                                      <  2:2533800/150‑1 (MQ=255)
                 ggCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTc                                                                                                                                      >  1:5299917/1‑150 (MQ=255)
                 ggCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTc                                                                                                                                      >  1:1066992/1‑150 (MQ=255)
                  gCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCg                                                                                                                                     <  2:4874722/150‑1 (MQ=255)
                   ccAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGt                                                                                                                                    >  1:569262/1‑150 (MQ=255)
                    cAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGta                                                                                                                                   <  2:5094366/150‑1 (MQ=255)
                           aaGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGc                                                                                                                            >  1:2871872/1‑150 (MQ=255)
                                   cGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGcc                                                                                                                    <  2:424900/150‑1 (MQ=255)
                                    gAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCg                                                                                                                   >  2:4041286/1‑150 (MQ=255)
                                       ttCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgac                                                                                                                >  2:5175873/1‑150 (MQ=255)
                                          aaGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacgac                                                                                                             >  2:3290008/1‑150 (MQ=255)
                                          aaGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacgac                                                                                                             >  2:169557/1‑150 (MQ=255)
                                          aaGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacgac                                                                                                             >  1:4906048/1‑150 (MQ=255)
                                          aaGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacgac                                                                                                             >  2:4894682/1‑150 (MQ=255)
                                          aaGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacgac                                                                                                             >  1:393360/1‑150 (MQ=255)
                                          aaGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacgac                                                                                                             >  1:3829145/1‑150 (MQ=255)
                                            ggAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCCCACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCTTCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACaa                                                                                                           <  2:3927589/150‑1 (MQ=255)
                                            ggAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACaa                                                                                                           >  1:5856032/1‑150 (MQ=255)
                                              cccgcgTCGTGCAGGTGCGCGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgc                                                                                                         >  1:4663203/6‑150 (MQ=255)
                                              agaTGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgc                                                                                                         <  1:6802983/150‑1 (MQ=255)
                                               taTGGTCGCGCACCTGCACGAGGGGGGGGTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCTATGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgcg                                                                                                        <  1:3156952/149‑1 (MQ=255)
                                                aTGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgcgc                                                                                                       >  1:678230/1‑150 (MQ=255)
                                                 tGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGAGCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCGGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGct                                                                                                      <  1:1020331/150‑1 (MQ=255)
                                                   gTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACACCCACACCGCCGACGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGctac                                                                                                    <  1:545485/150‑1 (MQ=255)
                                                   gTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGctac                                                                                                    <  2:5609881/150‑1 (MQ=255)
                                                    tCGCGCACCTGCACGACGCGGGGCTTGAGGTGAGCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgac                                                                                                                <  2:4128341/137‑1 (MQ=255)
                                                     cgcgCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGctacta                                                                                                  <  2:569262/150‑1 (MQ=255)
                                                            cTGCACGAGGGGGGGCTTGTGGTGATCCTGGAGGTGTTCTACAACCACACCGCCGATGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCTACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAAc                                                                                           <  1:6443463/150‑1 (MQ=255)
                                                            cTGCACGACGGGGGGTTTGAGGTGATCGTGGAGGGGGTCTACAACCACACGGCCGAAGGCAACGCGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAAc                                                                                           <  1:869229/150‑1 (MQ=255)
                                                                  gacgGGGGGTTTGTGTTGATCCTGGAGGTGGTCTAAAACCACACCGCCTATGGAAAGGCGCGCGGCCCGACCATTTCTATGCGGGGAATCTACAACGCCTCTTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTcc                                                                                     <  1:1208452/150‑1 (MQ=255)
                                                                  gacgCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTcc                                                                                     >  2:6064869/1‑150 (MQ=255)
                                                                   ccgcgGGGCTTGTGGTTATCCTGGAGGTGGTCTACAACCACACCGCCGATGGCAACGCGCGCGGCCCGACCCTGTGGATGCGCGTCATGGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccg                                                                                    <  2:5306860/149‑1 (MQ=255)
                                                                     gcgGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGGTGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggc                                                                                  >  2:6177472/1‑150 (MQ=255)
                                                                          gCTTGTGGTGATCCTGGACGTGGTCTACAACCACACCGCCGATGGCAACGAGCGCGGCCCGACCCTGTCTATGCGCGGCTTCGACAACGCATCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggcaccgg                                                                             <  2:2039845/150‑1 (MQ=255)
                                                                             tGAGGTTATCCTGGACGTGGAAAAAAACCACACCCCCCAAGGAAACGACCCCGCCCCCACCCTGTCGATGCGCGCAATCCAAAACGCCTCGTACTACCCCCTGATGCCGTACGACAACCGCTACTACAACAACGACTCCGGCACCGGCaa                                                                          <  1:7241457/150‑1 (MQ=255)
                                                                                    aTCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTg                                                                   <  2:5465324/150‑1 (MQ=255)
                                                                                     tCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTgg                                                                  >  1:4562490/1‑150 (MQ=255)
                                                                                          gACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTg                                                             <  2:6366869/150‑1 (MQ=255)
                                                                                           cggTGGTCTAAAACAACACCGCCTAAGGAAACGCGCGGGGCCCGACCCTGTGTATGGGGGGATTCCACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                            <  1:2697240/148‑1 (MQ=255)
                                                                                           aCGTGTTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCTACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                            <  2:4700685/150‑1 (MQ=255)
                                                                                           aCGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACACGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                            <  2:6827806/150‑1 (MQ=255)
                                                                                             gTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCg                                                                                                                                     >  1:73366/1‑75 (MQ=255)
                                                                                             gTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCg                                                                                                                                     <  2:73366/75‑1 (MQ=255)
                                                                                               ggTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCa                                                        >  2:30591/1‑150 (MQ=255)
                                                                                                gTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAc                                                       >  2:896192/1‑150 (MQ=255)
                                                                                                  cTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAccc                                                     >  1:3656191/1‑150 (MQ=255)
                                                                                                  cTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAccc                                                     <  2:2872460/150‑1 (MQ=255)
                                                                                                  aaaaaaCCACACCGCGGACGGAAACGCGCGGGCCCCGCCCCTGTGGAGGCGCGGAATCGACAACGCCTCTTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAccc                                                     <  1:808182/146‑1 (MQ=255)
                                                                                                   tACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAcccc                                                    >  2:2257593/1‑150 (MQ=255)
                                                                                                   tACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAcccc                                                    >  1:6877771/1‑150 (MQ=255)
                                                                                                      aaCCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTCCTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGGCCTGAGCCCCCCCTGc                                                 >  2:3693620/1‑150 (MQ=255)
                                                                                                       aCCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCg                                                >  1:4132284/1‑150 (MQ=255)
                                                                                                       aCCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCg                                                >  1:5218726/1‑150 (MQ=255)
                                                                                                        ccACACCGCCGATGGCAACGCGCGCGGCCCCACCCTGGCGAGGCGGGCCATCCACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGt                                               <  1:2779389/150‑1 (MQ=255)
                                                                                                        ccACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGt                                               >  1:600042/1‑150 (MQ=255)
                                                                                                         cacaCCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACa                                                                                                                                             >  2:1251031/1‑55 (MQ=255)
                                                                                                         cacaCCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtg                                              >  1:2447553/1‑150 (MQ=255)
                                                                                                            acCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACCAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCGGCGtgctg                                           >  2:1495169/1‑150 (MQ=255)
                                                                                                            acCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtgctg                                           >  2:1456297/1‑150 (MQ=255)
                                                                                                            acCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtgctg                                           >  1:7261497/1‑150 (MQ=255)
                                                                                                            acCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtgctg                                           <  2:5689141/150‑1 (MQ=255)
                                                                                                                ccgAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGt                                       >  1:4216303/1‑150 (MQ=255)
                                                                                                                 cgACGGCAACGAGCGCGGCCCGACCCTGTCGAGGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGtt                                      <  2:1760715/150‑1 (MQ=255)
                                                                                                                  gAAGGCAACGCGCGCGGCCCGACCCTGTCGAGGCGCGGCATCGACACCGCCCCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTg                                     <  1:6058987/150‑1 (MQ=255)
                                                                                                                  gAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTg                                     >  2:5279848/1‑150 (MQ=255)
                                                                                                                     ggCAACGAGCGCGGCCCGACCGTGTCGTTGCGGGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTc                                  <  1:7132475/150‑1 (MQ=255)
                                                                                                                       cAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCAc                                >  2:2340952/1‑150 (MQ=255)
                                                                                                                       cAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCAc                                >  1:231270/1‑150 (MQ=255)
                                                                                                                        aaCGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCGCCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCAcc                               >  2:2191773/1‑150 (MQ=255)
                                                                                                                        aaCGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCAcc                               >  1:3109134/1‑150 (MQ=255)
                                                                                                                        aaCGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCAcc                               >  2:1217975/1‑150 (MQ=255)
                                                                                                                        aaCGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCAcc                               >  1:355854/1‑150 (MQ=255)
                                                                                                                           gAGCGCGGCCCGACCCGGTCGTTGCGCGGCAGCGACAACGCCTCGTACTACCGCCTGTTGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGAc                            <  1:4806842/150‑1 (MQ=255)
                                                                                                                               gcgGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgc                        >  2:4831531/1‑150 (MQ=255)
                                                                                                                                 ggCCCGCCCCGGTCGAGGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACACGCGCTACTCCCGCCACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgctt                      <  1:1911421/150‑1 (MQ=255)
                                                                                                                                 ggCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgctt                      >  2:6599449/1‑150 (MQ=255)
                                                                                                                                  gCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgcttc                     >  1:4932971/1‑150 (MQ=255)
                                                                                                                                    ccGACCCTTTCGATGCGCGGAATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgcttcgc                   <  1:1710437/150‑1 (MQ=255)
                                                                                                                                     cGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgcttcgct                  >  2:632759/1‑150 (MQ=255)
                                                                                                                                       aCCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTAc                >  2:7269456/1‑150 (MQ=255)
                                                                                                                                         ccTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTg              <  2:5790992/150‑1 (MQ=255)
                                                                                                                                          ttGTCGATGCCGCGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTgg             >  2:7065272/2‑150 (MQ=255)
                                                                                                                                          cTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCa                                         <  2:1304951/122‑1 (MQ=255)
                                                                                                                                          cTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCa                                         >  1:1304951/1‑122 (MQ=255)
                                                                                                                                            ggctaggcgcGGAATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGc           <  1:2018152/144‑1 (MQ=255)
                                                                                                                                              ctaTGCGCGGAATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCg         <  2:6244939/148‑1 (MQ=255)
                                                                                                                                                aggcgcGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGc       <  2:4627772/148‑1 (MQ=255)
                                                                                                                                                aTGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGc       >  2:774322/1‑150 (MQ=255)
                                                                                                                                                    gcgGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGaaa   >  2:5030525/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGAAAT  >  NC_002947/4579072‑4579370

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