Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,796,102 A→G 5.7% E185G (GAA→GGA) ‡ PP_3353 → sulfatase‑modifying protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,796,1020AG 5.7% 1246.9 / 38.0 386E185G (GAA→GGA) ‡PP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base A (189/175);  new base G (6/16);  total (195/191)
Fisher's exact test for biased strand distribution p-value = 2.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 399 total aligned reads displayed.

CGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATG  >  NC_002947/3795953‑3796251
                                                                                                                                                      |                                                                                                                                                     
cGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGa                                                                                                                                                       <  1:1654863/150‑1 (MQ=255)
 gTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGaa                                                                                                                                                      >  1:6801288/1‑150 (MQ=255)
 gTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGaa                                                                                                                                                      >  2:4459281/1‑150 (MQ=255)
       ggCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGccc                                                                                                                                                <  1:4472337/150‑1 (MQ=255)
         ccccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAg                                                                                                                                              >  2:5771230/1‑150 (MQ=255)
          cccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGc                                                                                                                                             >  1:753476/1‑150 (MQ=255)
          cccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGc                                                                                                                                             >  1:3758833/1‑150 (MQ=255)
           ccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGCg                                                                                                                                            <  2:6028142/150‑1 (MQ=255)
           ccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCg                                                                                                                                            >  1:4115892/1‑150 (MQ=255)
           ccGGCTGGCAGGTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGCg                                                                                                                                            <  2:1605181/150‑1 (MQ=255)
             ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCTCGCTGGGCAGTTCGCGGCCTGCCCAGCgg                                                                                                                                           <  1:2197246/150‑2 (MQ=255)
             ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCTCGCTGGGCAGTTCGCGGCCTGCCCAGCgg                                                                                                                                           <  1:6641045/150‑2 (MQ=255)
             ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGcagg                                                                                                                                          <  2:1878398/150‑4 (MQ=255)
                 ggCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCTCGCTGGGCAGTTCGCGGCCTGCCCAGCgggcata                                                                                                                                      <  2:3322757/150‑7 (MQ=39)
                 ggCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGgcgc                                                                                                                                      <  1:1824863/150‑1 (MQ=255)
                 ggCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGgcgc                                                                                                                                      <  2:745433/150‑1 (MQ=255)
                  gCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCa                                                                                                                                     >  2:1035815/1‑150 (MQ=255)
                    aGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTAGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGc                                                                                                                                   >  1:6557454/1‑150 (MQ=255)
                     gTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCt                                                                                                                                  <  1:5909108/150‑1 (MQ=255)
                       tCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTgg                                                                                                                                <  1:5144440/150‑1 (MQ=255)
                         gTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGaa                                                                                                                              >  2:6230944/1‑150 (MQ=255)
                           gCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATa                                                                                                                            >  1:3551903/1‑150 (MQ=255)
                             ccGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATAcc                                                                                                                          >  2:3687187/1‑149 (MQ=255)
                              cGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGc                                                                                                                         <  1:459446/150‑1 (MQ=255)
                              cGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGc                                                                                                                         <  2:7275021/150‑1 (MQ=255)
                                     aaCTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgcgc                                                                                                                  >  1:6381222/1‑150 (MQ=255)
                                     aaCTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgcgc                                                                                                                  >  1:8336822/1‑150 (MQ=255)
                                     aaCTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgcgc                                                                                                                  >  2:2795695/1‑150 (MQ=255)
                                       cTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcgg                                                                                                                <  1:5554976/150‑1 (MQ=255)
                                          gCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGcggcgg                                                                                                             >  2:5655126/1‑150 (MQ=255)
                                            gACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGcc                                                                                                           >  2:735735/1‑150 (MQ=255)
                                            gACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGcc                                                                                                           >  1:1173083/1‑150 (MQ=255)
                                                 ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAcc                                                                                                      >  1:3654399/1‑150 (MQ=255)
                                                  cagcccggCCAGGCTGCTGCCCGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCg                                                                                                     >  2:5205907/7‑150 (MQ=255)
                                                  cagcccggCCAGGCGGCTGCCCGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGACTACGCCATGCGCGGCGGCCTGACCg                                                                                                     >  2:910942/7‑150 (MQ=255)
                                                  cGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCg                                                                                                     >  1:1049677/1‑150 (MQ=255)
                                                   gCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGa                                                                                                    >  2:2189864/1‑150 (MQ=255)
                                                    cAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGAt                                                                                                   <  1:6223363/150‑1 (MQ=255)
                                                      aGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGc                                                                                                 <  1:8437554/150‑1 (MQ=255)
                                                       gggCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGcc                                                                                                >  1:2000799/1‑150 (MQ=255)
                                                       gggCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGcc                                                                                                >  2:4862242/1‑150 (MQ=255)
                                                         gCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGa                                                                                              <  2:1079421/150‑1 (MQ=255)
                                                             ggCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTc                                                                                          <  2:663687/150‑1 (MQ=255)
                                                              gcagcaGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCa                                                                                         >  2:6574297/1‑150 (MQ=255)
                                                              gcagcaGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCa                                                                                         <  2:3143456/150‑1 (MQ=255)
                                                                 gcagcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCt                                                                                      <  1:3951515/150‑1 (MQ=255)
                                                                   agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                                    <  2:8203827/150‑1 (MQ=255)
                                                                   agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                                    <  2:7034628/150‑1 (MQ=255)
                                                                    gcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTggg                                                                                   <  2:8332518/150‑1 (MQ=255)
                                                                              gCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGa                                                                         <  2:2995122/150‑1 (MQ=255)
                                                                              gCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGAGGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCATCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGa                                                                         <  1:8727968/150‑1 (MQ=255)
                                                                                tGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagc                                                                       >  2:5883848/1‑150 (MQ=255)
                                                                                  gACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcag                                                                     >  1:78915/1‑150 (MQ=255)
                                                                                  gACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcag                                                                     >  1:2703305/1‑150 (MQ=255)
                                                                                        cctCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAg                                                               <  2:8462523/148‑1 (MQ=255)
                                                                                        cATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAg                                                               <  2:802452/150‑1 (MQ=255)
                                                                                         aTCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAgg                                                              <  2:3482465/150‑1 (MQ=255)
                                                                                         aTCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAgg                                                              >  2:1234568/1‑150 (MQ=255)
                                                                                         aTCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAgg                                                              <  1:7114323/150‑1 (MQ=255)
                                                                                          tCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAggg                                                             <  1:1918049/150‑1 (MQ=255)
                                                                                                 gtgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTc                                                      >  2:430439/1‑150 (MQ=255)
                                                                                                  tgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCa                                                     >  1:3149362/1‑150 (MQ=255)
                                                                                                   gCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCAt                                                    >  1:227226/1‑150 (MQ=255)
                                                                                                   gCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCAt                                                    <  1:171186/150‑1 (MQ=255)
                                                                                                     aGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATgg                                                  <  1:4649610/150‑1 (MQ=255)
                                                                                                      gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                                 <  2:7711796/150‑1 (MQ=255)
                                                                                                      gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                                 <  2:7921104/150‑1 (MQ=255)
                                                                                                          gCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAAt                                             >  1:1285403/1‑150 (MQ=255)
                                                                                                             cTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATAcc                                          >  2:5051871/1‑150 (MQ=255)
                                                                                                              tGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCt                                         >  1:4030892/1‑150 (MQ=255)
                                                                                                              tGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCt                                         >  1:3425209/1‑150 (MQ=255)
                                                                                                               ggAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTg                                        >  1:2920485/1‑150 (MQ=255)
                                                                                                                gAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTgg                                       >  2:2521189/1‑150 (MQ=255)
                                                                                                                  aGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCa                                     <  1:8053252/150‑1 (MQ=255)
                                                                                                                    aTGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAgg                                   <  1:6064438/150‑1 (MQ=255)
                                                                                                                       cccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTc                                >  2:2225975/1‑150 (MQ=255)
                                                                                                                        ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCa                               >  2:6981354/1‑150 (MQ=255)
                                                                                                                        ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCa                               >  2:8647505/1‑150 (MQ=255)
                                                                                                                        ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCa                               >  1:8621854/1‑150 (MQ=255)
                                                                                                                        ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCa                               >  1:429055/1‑150 (MQ=255)
                                                                                                                        ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCa                               <  1:8319129/150‑1 (MQ=255)
                                                                                                                        ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCa                               >  1:7614990/1‑150 (MQ=255)
                                                                                                                           ggCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGtt                            <  1:2521189/150‑1 (MQ=255)
                                                                                                                            gCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTc                           <  1:3248440/150‑1 (MQ=255)
                                                                                                                             ccTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTcc                          <  2:2028489/150‑1 (MQ=255)
                                                                                                                              cTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTccc                         <  1:8217864/150‑1 (MQ=255)
                                                                                                                                   cctcGCTGGGCAGTTCGCGGCCTGCCCAGCGGGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTAt                     >  1:7434345/4‑150 (MQ=255)
                                                                                                                                   cccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTCCCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATc                    >  1:5444537/1‑150 (MQ=255)
                                                                                                                                    ctcGCTGGGCAGTTCGCGGCCTGCCCAGCAAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATc                    >  1:2104221/3‑150 (MQ=255)
                                                                                                                                         tGGGCAGTTCGCGGCCTGCCCAGCGGGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAAt               >  1:1464465/1‑150 (MQ=255)
                                                                                                                                         tGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAAt               >  1:8798074/1‑150 (MQ=255)
                                                                                                                                           ggCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAACACCTGGCAGGGTCAGTTCCCTTATCGCAATgc             >  2:2436414/1‑150 (MQ=255)
                                                                                                                                              aGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgcggc          <  2:2633256/150‑1 (MQ=255)
                                                                                                                                                 ccGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGTAGGGTCAGTTCCCTTATCGCAATGCGGCGaa       <  2:6874989/150‑1 (MQ=255)
                                                                                                                                                 ccGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGaa       <  2:5457382/150‑1 (MQ=255)
                                                                                                                                                 ccGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGaa       >  2:3274103/1‑150 (MQ=255)
                                                                                                                                                 ccGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGaa       <  1:7406013/150‑1 (MQ=255)
                                                                                                                                                  cgcgAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCGGTTCCCTTATCGCAATGCGGCGAAg      >  1:2831947/1‑150 (MQ=255)
                                                                                                                                                   gcgAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAgg     >  2:1449199/1‑150 (MQ=255)
                                                                                                                                                      aaCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATg  >  1:6461912/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATG  >  NC_002947/3795953‑3796251

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