Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,791,148 C→T 5.6% G183G (GGC→GGT PP_3353 → sulfatase‑modifying protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,791,1480CT 5.6% 305.7 / 4.5 100G183G (GGC→GGTPP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base C (48/46);  new base T (2/4);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCG  >  NC_002947/3791001‑3791292
                                                                                                                                                   |                                                                                                                                                
ggACGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCcg                                                                                                                                                >  2:1609541/1‑150 (MQ=255)
ggACGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCcg                                                                                                                                                >  2:1086010/1‑150 (MQ=255)
 gACGGGCTCGGCCCCGGCTTGCGGTTTGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCcgc                                                                                                                                               <  2:1589273/150‑1 (MQ=255)
    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                                                                                                                         >  2:2111489/1‑82 (MQ=255)
                                                                                           cATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGc                                                                                                                         <  1:2111489/82‑1 (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                                       <  2:1484681/150‑1 (MQ=255)
                                                                                                         gATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGc                                       <  1:565830/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                                                                     <  1:859679/113‑1 (MQ=255)
                                                                                                                cTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTac                                                                     >  2:859679/1‑113 (MQ=255)
                                                                                                                cTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATAcc                                >  1:1994327/1‑150 (MQ=255)
                                                                                                                cTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATAcc                                <  2:1047606/150‑1 (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                                                                                >  2:2030459/1‑91 (MQ=255)
                                                                                                                           ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTc                                                                                <  1:2030459/91‑1 (MQ=255)
                                                                                                                           ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAg                                                                              <  2:1949449/93‑1 (MQ=255)
                                                                                                                           ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAg                                                                              >  1:1949449/1‑93 (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                <  2:435248/150‑1 (MQ=255)
                                                                                                                                ccTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTcc                <  1:2407808/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         <  2:1994327/150‑1 (MQ=255)
                                                                                                                                       ccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCg         >  1:1662925/1‑150 (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)
                                                                                                                                                   |                                                                                                                                                
GGACGTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCG  >  NC_002947/3791001‑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: