Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,791,154 A→C 5.1% 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 5.1% 287.4 / 3.2 96E185D (GAA→GAC) ‡PP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base A (45/46);  new base C (2/3);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCG  >  NC_002947/3791005‑3791292
                                                                                                                                                     |                                                                                                                                          
gTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGaa                                                                                                                                            >  2:1544657/1‑150 (MQ=255)
       gCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCa                                                                                                                                     <  1:1544657/150‑1 (MQ=255)
        ccccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAg                                                                                                                                    >  1:965973/1‑150 (MQ=255)
         cccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGc                                                                                                                                   <  2:2248124/150‑1 (MQ=255)
          ccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCg                                                                                                                                  >  1:772197/1‑150 (MQ=255)
          ccGACTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCg                                                                                                                                  >  1:936855/1‑150 (MQ=255)
            ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAGGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAg                                                                                                                                <  2:965973/150‑1 (MQ=255)
            ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAg                                                                                                                                <  1:1220689/150‑1 (MQ=255)
            ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAg                                                                                                                                <  2:223316/150‑1 (MQ=255)
            ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGACTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAg                                                                                                                                    >  2:1927004/1‑146 (MQ=255)
            ggCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGACTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAg                                                                                                                                    <  1:1927004/146‑1 (MQ=255)
             gCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGTCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAgg                                                                                                                               >  2:2229158/1‑150 (MQ=255)
               tGGCAGTTCGTTCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCTTGGCTCTTCCGATCTTCCATCCGGTGGTGCAGTTCGCACCCGAAGAAGCAAAGCCCCATCCCCGGTGGGCAGTTCGCGGCCTGCCCAGCgggcat                                                                                                                             <  1:550881/150‑6 (MQ=255)
                            ccGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACg                                                                                                                >  2:2407808/1‑150 (MQ=255)
                             cGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGc                                                                                                               >  1:435248/1‑150 (MQ=255)
                                    aaCTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgcgc                                                                                                        >  1:2052729/1‑150 (MQ=255)
                                    aaCTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgcgc                                                                                                        <  1:991085/150‑1 (MQ=255)
                                       tGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcggc                                                                                                     >  2:624266/1‑150 (MQ=255)
                                       tGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACACGTTGCgcggc                                                                                                     >  1:1308800/1‑150 (MQ=255)
                                       tGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGAAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcggc                                                                                                     >  2:1497195/1‑150 (MQ=255)
                                          cGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCggcggc                                                                                                  <  1:1086010/150‑1 (MQ=255)
                                           gACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAGGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGcc                                                                                                 <  2:533846/150‑1 (MQ=255)
                                           gACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGcc                                                                                                 >  2:2146732/1‑150 (MQ=255)
                                              acCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGa                                                                                              <  1:1840770/150‑1 (MQ=255)
                                               cccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAc                                                                                             <  2:1938670/150‑1 (MQ=255)
                                                ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGaa                                                                                                                                            >  2:774035/1‑102 (MQ=255)
                                                ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGaa                                                                                                                                            <  1:774035/102‑1 (MQ=255)
                                                ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAcc                                                                                            >  1:2083177/1‑150 (MQ=255)
                                                     aGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGc                                                                                       <  1:2146732/150‑1 (MQ=255)
                                                      gggCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGcc                                                                                      >  2:1948015/1‑150 (MQ=255)
                                                             atctcagcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCCACCGCGGCGGCCTGACCGATGCCGACTTCa                                                                               <  1:1270548/146‑1 (MQ=255)
                                                              cagcagCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAg                                                                              >  1:414551/1‑150 (MQ=255)
                                                                gcagcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCt                                                                            >  2:2088039/1‑150 (MQ=255)
                                                                gcagcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCt                                                                            <  1:2322096/150‑1 (MQ=255)
                                                                  agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                          <  2:1251193/150‑1 (MQ=255)
                                                                  agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCACCGAGCGGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAgagat                                                                          >  2:1270548/1‑146 (MQ=255)
                                                                   gcCTAGCCGCGCTGGACAACCATCCGGTGGTGAAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAAGTGCGCAGCGCGGCGCGGATGGAATACCCCGTGCGCGGCGGTTTGACCGATGCCGACTTCAGCTggg                                                                         <  1:1025530/150‑1 (MQ=255)
                                                                    ccTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgggg                                                                        >  2:1825050/1‑150 (MQ=255)
                                                                         ccGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTacca                                                                   <  1:440635/150‑1 (MQ=255)
                                                                           gggCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTaccacc                                                                 <  1:392246/150‑1 (MQ=255)
                                                                              cTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTaccacc                                                                 >  1:186896/1‑147 (MQ=255)
                                                                              cTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTaccacc                                                                 <  2:186896/147‑1 (MQ=255)
                                                                              cTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGag                                                              <  2:772197/150‑1 (MQ=255)
                                                                                 gACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcag                                                           <  1:624266/150‑1 (MQ=255)
                                                                                 gACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcag                                                           >  1:749863/1‑150 (MQ=255)
                                                                                 gACAACCAGCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcag                                                           <  1:1482401/150‑1 (MQ=255)
                                                                                    aaCCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGccc                                                        >  1:671679/1‑150 (MQ=255)
                                                                                      ccATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCaa                                                      >  1:1047606/1‑150 (MQ=255)
                                                                                      ccATCCGGTGATGCAGATCTCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAccccca                                                      <  2:1803519/150‑3 (MQ=255)
                                                                                       cATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGc                                                                                                                         <  1:2111489/82‑1 (MQ=255)
                                                                                       cATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGc                                                                                                                         >  2:2111489/1‑82 (MQ=255)
                                                                                         tCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAggg                                                   >  1:1135755/1‑150 (MQ=255)
                                                                                             gtggtgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAg                                               >  2:1981506/1‑150 (MQ=255)
                                                                                              ctcttccgATCTCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGTCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGc                                              <  1:1353078/143‑1 (MQ=255)
                                                                                               ggtgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTTGTCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCt                                             >  1:1260244/1‑150 (MQ=255)
                                                                                                gtgCAGATCGCCCTGAAAGATGCCCAGACCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTTAGCTGGGGCCCCACCGCGCAGCCCCAGGGCAAGCTc                                            >  2:313597/1‑150 (MQ=255)
                                                                                                    aGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATgg                                        <  2:1712083/150‑1 (MQ=255)
                                                                                                     gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                       <  1:565830/150‑1 (MQ=255)
                                                                                                     gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                       <  2:1484681/150‑1 (MQ=255)
                                                                                                     gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                       >  2:166548/1‑150 (MQ=255)
                                                                                                      aTCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGcc                                      >  1:1073304/1‑150 (MQ=255)
                                                                                                          cccTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGAACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATa                                  <  1:1948015/150‑1 (MQ=255)
                                                                                                           ccTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATAc                                 >  2:59339/1‑150 (MQ=255)
                                                                                                            cTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTac                                                                     >  2:859679/1‑113 (MQ=255)
                                                                                                            cTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTac                                                                     <  1:859679/113‑1 (MQ=255)
                                                                                                            cTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATAcc                                <  2:1047606/150‑1 (MQ=255)
                                                                                                            cTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATAcc                                >  1:1994327/1‑150 (MQ=255)
                                                                                                             tGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCCGAACGGTATCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCt                               <  2:933715/150‑1 (MQ=255)
                                                                                                               gAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAg                                                                                                                                >  1:1494059/1‑51 (MQ=255)
                                                                                                               gAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAg                                                                                                                                <  2:1494059/51‑1 (MQ=255)
                                                                                                                   aTGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGTGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAgg                         <  2:344813/150‑1 (MQ=255)
                                                                                                                   aTGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAgg                         <  1:2088039/150‑1 (MQ=255)
                                                                                                                       ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTc                                                                                <  1:2030459/91‑1 (MQ=255)
                                                                                                                       ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTc                                                                                >  2:2030459/1‑91 (MQ=255)
                                                                                                                       ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAg                                                                              >  1:1949449/1‑93 (MQ=255)
                                                                                                                       ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAg                                                                              <  2:1949449/93‑1 (MQ=255)
                                                                                                                       ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTggg                                                                         >  1:1038816/1‑98 (MQ=255)
                                                                                                                       ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTggg                                                                         <  2:1038816/98‑1 (MQ=255)
                                                                                                                       ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCa                     >  1:2204551/1‑150 (MQ=255)
                                                                                                                          ggCCTATGCCCGCTGGGCAGGCCGCGAACTGCCAAGCGAGGCGCAGCTGGAATACGCCATGCGCGCCGGCCTGCCCGATGCTGACTGCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGtt                  <  2:1548113/150‑1 (MQ=255)
                                                                                                                            ccTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTcc                <  1:2407808/150‑1 (MQ=255)
                                                                                                                            ccTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTcc                <  2:435248/150‑1 (MQ=255)
                                                                                                                                 gCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTAt           >  2:770598/1‑150 (MQ=255)
                                                                                                                                  ctcGCTGGGCAGTTCGCGGCCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCt                               <  2:1756774/127‑1 (MQ=255)
                                                                                                                                  ctcGCTGGGCAGTTCGCGGCCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCt                               >  1:1756774/3‑129 (MQ=255)
                                                                                                                                  cccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATc          <  2:2204551/150‑1 (MQ=255)
                                                                                                                                   ccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTccc               <  1:2399498/144‑1 (MQ=255)
                                                                                                                                   ccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTccc               >  2:2399498/1‑144 (MQ=255)
                                                                                                                                   ccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTTCCTTATCg         >  1:911784/1‑150 (MQ=255)
                                                                                                                                   ccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCg         <  1:2229158/150‑1 (MQ=255)
                                                                                                                                   ccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCg         >  1:1662925/1‑150 (MQ=255)
                                                                                                                                   ccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCg         <  2:1994327/150‑1 (MQ=255)
                                                                                                                                     gCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCTGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCa       >  2:1560077/1‑150 (MQ=255)
                                                                                                                                      cTGGGCAGTTCGCGGCCTGGCCAGCGAGGCGAAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCaa      >  1:1665026/1‑150 (MQ=255)
                                                                                                                                       tGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAAt     >  1:1224550/1‑150 (MQ=255)
                                                                                                                                         ggCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgc   >  1:2015296/1‑150 (MQ=255)
                                                                                                                                          gCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgcg  <  1:166548/150‑1 (MQ=255)
                                                                                                                                                cgcgAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGt                       >  1:1773367/1‑123 (MQ=255)
                                                                                                                                                cgcgAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGt                       <  2:1773367/123‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
GTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCG  >  NC_002947/3791005‑3791292

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: