Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,791,154 A→C 7.5% E185D (GAA→GAC) ‡ PP_3353 → sulfatase‑modifying protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,791,1540AC 7.5% 374.0 / 16.2 137E185D (GAA→GAC) ‡PP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base A (62/64);  new base C (5/6);  total (67/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

GTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATG  >  NC_002947/3791005‑3791302
                                                                                                                                                     |                                                                                                                                                    
gTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACAAGCCATTAATATTCATTTCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGaa                                                                                                                                                      <  2:3050301/150‑1 (MQ=255)
gTGCTCGGACCCGGCTGGCAGTTAGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGaa                                                                                                                                                      >  2:2340940/1‑150 (MQ=255)
 tGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAAc                                                                                                                                                     >  2:1950011/1‑150 (MQ=255)
   cTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTg                                                                                                                                                   <  1:120425/150‑1 (MQ=255)
    tCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGc                                                                                                                                                  <  1:731004/150‑1 (MQ=255)
         cccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGc                                                                                                                                             <  2:328977/150‑1 (MQ=255)
         cccGGCTGGCAGTTCGGGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCGGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCGCGGCCTGTGCCCGCTGGGCAGGCCGCGACCTGCCCAGc                                                                                                                                             >  2:1402909/1‑150 (MQ=255)
             gCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGCgggc                                                                                                                                         >  2:2085268/1‑147 (MQ=255)
              cTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGgc                                                                                                                                        <  1:3053810/150‑1 (MQ=255)
              cTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGgc                                                                                                                                        >  1:2923759/1‑150 (MQ=255)
                 gCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGCgggcatag                                                                                                                                     <  1:2085268/150‑8 (MQ=255)
                 gCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCa                                                                                                                                     <  2:2902841/150‑1 (MQ=255)
                 gCAGTTCGAGCCCGGCGCCAACTGGCGACACCCGCAAGGGACGGGGTGACTGGAGTTCAGACGTGTGATCTTTCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCTGCTGGGCAGTTCGCGGCCTGCCCAGCgggcatag                                                                                                                                     <  1:361193/150‑8 (MQ=255)
                  cAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAg                                                                                                                                    >  1:1602092/1‑150 (MQ=255)
                   aGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGc                                                                                                                                   <  2:669136/150‑1 (MQ=255)
                   aGTTCGTGCCCGGCGCCAACTGGCGAAACCCGCACGGGCGGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGc                                                                                                                                   <  2:2639679/150‑1 (MQ=255)
                    gTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCt                                                                                                                                  <  1:513939/150‑1 (MQ=255)
                     ttCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTg                                                                                                                                 >  1:641634/1‑150 (MQ=255)
                     ttCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTg                                                                                                                                 <  1:961021/150‑1 (MQ=255)
                            ccGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACg                                                                                                                          >  2:2931680/1‑150 (MQ=255)
                                    aaCTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGCg                                                                                                                                            >  2:1689083/1‑124 (MQ=255)
                                    aaCTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGCg                                                                                                                                            <  1:1689083/124‑1 (MQ=255)
                                    aaCTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgcgc                                                                                                                  >  1:2102027/1‑150 (MQ=255)
                                    aaCTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgcgc                                                                                                                  >  2:2087519/1‑150 (MQ=255)
                                         gCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCAt                                                                                                                      <  1:1976230/141‑1 (MQ=255)
                                         gCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCAt                                                                                                                      >  2:1976230/1‑141 (MQ=255)
                                           gACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGcc                                                                                                           >  1:1996882/1‑150 (MQ=255)
                                              acCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGa                                                                                                        >  1:2437366/1‑150 (MQ=255)
                                               cccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAc                                                                                                       <  1:675409/150‑1 (MQ=255)
                                                ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAg                                                                                                                                          >  1:2961489/1‑114 (MQ=255)
                                                ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAg                                                                                                                                          <  2:2961489/114‑1 (MQ=255)
                                                ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAcc                                                                                                      >  2:2176267/1‑150 (MQ=255)
                                                ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAcc                                                                                                      >  2:1561923/1‑150 (MQ=255)
                                                  gCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGa                                                                                                    >  1:2505703/1‑150 (MQ=255)
                                                   cAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGAt                                                                                                   <  1:1561923/150‑1 (MQ=255)
                                                    aaGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATg                                                                                                  >  1:26173/1‑150 (MQ=255)
                                                     aGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGc                                                                                                 >  1:1933121/1‑150 (MQ=255)
                                                          cGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACt                                                                                            >  1:1037793/1‑150 (MQ=255)
                                                           gggCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGaa                                                                                                                                                      <  1:567148/91‑1 (MQ=255)
                                                                gcagcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAg                                                                                                                                          >  2:496629/1‑98 (MQ=255)
                                                                gcagcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGc                                                                                                                                   <  1:2941774/105‑1 (MQ=255)
                                                                 cagcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTg                                                                                     <  2:887/150‑1 (MQ=255)
                                                                  agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                                    <  2:1163102/150‑1 (MQ=255)
                                                                  agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                                    >  2:2674638/1‑150 (MQ=255)
                                                                  agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                                    <  2:2852757/150‑1 (MQ=255)
                                                                         tgcgcgtGGACTACCATCCGGTGGTGCAGAACGCCCTGGAAGATGCCCAGGCCTATGCCCGCTTGGCACGCTGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTacca                                                                             <  1:1134202/144‑1 (MQ=255)
                                                                             gCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGa                                                                         <  1:1054267/150‑1 (MQ=255)
                                                                                 gACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGcc                                                                                                <  1:1062375/123‑1 (MQ=255)
                                                                                 gACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGcc                                                                                                >  2:1062375/1‑123 (MQ=255)
                                                                                 gACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcag                                                                     >  2:89960/1‑150 (MQ=255)
                                                                                  aCAACCATCCGGTGGTGCAGATCGCCGTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCATCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcagc                                                                    >  2:1217160/1‑150 (MQ=255)
                                                                                  aCAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcagc                                                                    <  1:89960/150‑1 (MQ=255)
                                                                                  aCAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcagc                                                                    <  1:1950011/150‑1 (MQ=255)
                                                                                      ccATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAg                                                               >  2:1727443/1‑150 (MQ=255)
                                                                                        aTCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAgg                                                              <  1:2087519/150‑1 (MQ=255)
                                                                                             gtggtgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGcc                                                                                                >  1:823658/1‑111 (MQ=255)
                                                                                             gtggtgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGcc                                                                                                <  2:823658/111‑1 (MQ=255)
                                                                                                gtgCAGATCTCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGAACCCATCTGGGGGACCACCGAGCAGCCCAAGGGCAAGCTc                                                      >  1:2941238/1‑150 (MQ=255)
                                                                                                gtgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTc                                                      >  1:593195/1‑150 (MQ=255)
                                                                                                gtgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTc                                                      >  1:39230/1‑150 (MQ=255)
                                                                                                   cataTCGCCCCGGAATACCCCCAGGCCTAGGCCCTCTGGCCAGGCCGCGAACTGCCCAGCGAGGCGCATCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATg                                                   <  2:641634/147‑1 (MQ=255)
                                                                                                     taTCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                                 <  1:2977381/149‑1 (MQ=255)
                                                                                                     gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCt                                                                                                          <  2:2542215/93‑1 (MQ=255)
                                                                                                     gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCt                                                                                                          >  1:2542215/1‑93 (MQ=255)
                                                                                                     gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                                 >  1:2324263/1‑150 (MQ=255)
                                                                                                     gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                                 <  2:319479/150‑1 (MQ=255)
                                                                                                     gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                                 <  2:521363/150‑1 (MQ=255)
                                                                                                          cccTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATa                                            >  1:2146618/1‑150 (MQ=255)
                                                                                                           ccTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATAc                                           >  2:1605307/1‑150 (MQ=255)
                                                                                                            cTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATAcc                                          <  2:1788266/150‑1 (MQ=255)
                                                                                                             tGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCGGCTg                                                                                                                                 >  1:177648/1‑62 (MQ=255)
                                                                                                              ggAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTg                                        >  2:67729/1‑150 (MQ=255)
                                                                                                               gAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTCCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTgg                                       <  2:2102027/150‑1 (MQ=255)
                                                                                                                  gATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAg                                    <  2:1933121/150‑1 (MQ=255)
                                                                                                                  gATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAg                                    <  1:799108/150‑1 (MQ=255)
                                                                                                                    tGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAggg                                  >  1:1660100/1‑150 (MQ=255)
                                                                                                                     ccccAGGCCTCTGCCCGCTGGACACGCGGACCCATACCCAGCGAGGCCCAGCTCGAATAAGCCATGCGCCCCCGCCTGCCTGATTCCTACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGt                                 <  2:2941238/149‑1 (MQ=255)
                                                                                                                       ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAgg                                                              >  1:385764/1‑119 (MQ=255)
                                                                                                                       ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCa                               >  2:1793012/1‑150 (MQ=255)
                                                                                                                         aGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGt                             >  1:1445598/1‑150 (MQ=255)
                                                                                                                         aGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGt                             >  1:134622/1‑150 (MQ=255)
                                                                                                                             cTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCAt                                                    <  1:2099402/123‑1 (MQ=255)
                                                                                                                             cTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCAt                                                    >  2:2099402/1‑123 (MQ=255)
                                                                                                                             cTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACAGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTccc                         <  1:2198309/150‑1 (MQ=255)
                                                                                                                                 gCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTAt                     <  2:2437366/150‑1 (MQ=255)
                                                                                                                                 cctcGCTGGGCAGTTCGCGGCCTGCCCAGCGGGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCt                                         <  1:671784/127‑1 (MQ=38)
                                                                                                                                  ctcGCTGGGCAGTTCGCGGCCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATc                    >  1:1770172/3‑150 (MQ=255)
                                                                                                                                   ccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCg                   >  1:2911877/1‑150 (MQ=255)
                                                                                                                                    cGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGc                  >  2:1697229/1‑150 (MQ=255)
                                                                                                                                         ggCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgc             <  1:1697229/150‑1 (MQ=255)
                                                                                                                                          gCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcagc                                                                    <  1:314364/94‑1 (MQ=255)
                                                                                                                                           cAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCCTGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTCGCAGGGTCAGTTCCCTTATCGCAATgcgg           <  1:2186088/150‑1 (MQ=255)
                                                                                                                                              gCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGa        >  2:37282/1‑150 (MQ=255)
                                                                                                                                               ccGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGaa       >  1:2514503/1‑150 (MQ=255)
                                                                                                                                               ccGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGaa       <  2:1445598/150‑1 (MQ=255)
                                                                                                                                               ccGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGaa       <  2:1037793/150‑1 (MQ=255)
                                                                                                                                                cgcgAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAg      >  2:2165666/1‑150 (MQ=255)
                                                                                                                                                 gcgAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCACTTCTCGAAATGCGACGAAgg     >  2:2474152/1‑150 (MQ=255)
                                                                                                                                                    aaCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATg  <  2:2324263/150‑1 (MQ=255)
                                                                                                                                                    aaCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACGTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATg  <  1:2340940/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATG  >  NC_002947/3791005‑3791302

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