Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,793,076 A→C 5.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,793,0760AC 5.5% 657.9 / 19.3 193E185D (GAA→GAC) ‡PP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base A (94/87);  new base C (5/6);  total (99/94)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 197 total aligned reads displayed.

GTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATGG  >  NC_002947/3792927‑3793225
                                                                                                                                                     |                                                                                                                                                     
gTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGaa                                                                                                                                                       <  2:4141151/150‑1 (MQ=255)
  gCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACt                                                                                                                                                     <  2:1569710/150‑1 (MQ=255)
   cTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTg                                                                                                                                                    >  2:817391/1‑150 (MQ=255)
   cTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTg                                                                                                                                                    <  2:3890425/150‑1 (MQ=255)
   cTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTg                                                                                                                                                    >  1:2729915/1‑150 (MQ=255)
   cTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTg                                                                                                                                                    >  2:4007806/1‑150 (MQ=255)
    tCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGc                                                                                                                                                   <  1:3130269/150‑1 (MQ=255)
     cGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGcc                                                                                                                                                  <  2:3859027/150‑1 (MQ=255)
       gCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCa                                                                                                                                                <  2:2943766/150‑1 (MQ=255)
       gCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCa                                                                                                                                                <  2:3213006/150‑1 (MQ=255)
       gCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCa                                                                                                                                                <  2:4163598/150‑1 (MQ=255)
        ccccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAg                                                                                                                                               >  1:835514/1‑150 (MQ=255)
         cccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGc                                                                                                                                              <  1:3559856/150‑1 (MQ=255)
         cccGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGc                                                                                                                                              >  2:279787/1‑150 (MQ=255)
           cGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGCgg                                                                                                                                            <  2:4136869/150‑2 (MQ=255)
           cGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGa                                                                                                                                            <  2:3467558/150‑1 (MQ=255)
           cGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGa                                                                                                                                            <  1:1533319/150‑1 (MQ=255)
           cGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGa                                                                                                                                            <  1:3298354/150‑1 (MQ=255)
           cGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGa                                                                                                                                            <  1:1008674/150‑1 (MQ=255)
                 gCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCa                                                                                                                                      <  1:89335/150‑1 (MQ=255)
                   aGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCGTCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGc                                                                                                                                    <  1:931904/150‑1 (MQ=255)
                            ccGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACg                                                                                                                           >  2:4051197/1‑150 (MQ=255)
                            ccGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACg                                                                                                                           >  1:1993872/1‑150 (MQ=255)
                                 gccagttGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATg                                                                                                                      >  1:3967550/7‑150 (MQ=255)
                                  ccAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATgc                                                                                                                     >  2:1151641/1‑150 (MQ=255)
                                      cTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcgg                                                                                                                 <  1:2125536/150‑1 (MQ=255)
                                       tGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcggc                                                                                                                >  1:3937475/1‑150 (MQ=255)
                                       tGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcggc                                                                                                                >  2:380588/1‑150 (MQ=255)
                                        ggCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcggcg                                                                                                               <  2:2631617/150‑1 (MQ=255)
                                          cGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCggcggc                                                                                                             <  2:193581/150‑1 (MQ=255)
                                           gACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGcc                                                                                                            >  2:3998296/1‑150 (MQ=255)
                                              acCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATAGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGa                                                                                                         >  2:2026211/1‑150 (MQ=255)
                                               cccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAc                                                                                                        <  1:3883801/150‑1 (MQ=255)
                                               cccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAc                                                                                                        <  1:1794454/150‑1 (MQ=255)
                                                ccGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGAcc                                                                                                       >  2:1245089/1‑150 (MQ=255)
                                                 cGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGTTCGCGGCCTGCCCAGCg                                                                                                                                             >  1:2958171/1‑111 (MQ=39)
                                                 cGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGACTACGCCATGCGCGGCGGCCTGACCg                                                                                                      >  2:2405832/1‑150 (MQ=255)
                                                  gCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGa                                                                                                     >  2:2263035/1‑150 (MQ=255)
                                                  gCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGa                                                                                                     <  1:3635535/150‑1 (MQ=255)
                                                    aaGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATg                                                                                                   >  2:2069903/1‑150 (MQ=255)
                                                     aGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGc                                                                                                  <  2:2341177/150‑1 (MQ=255)
                                                      gggCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGcc                                                                                                 <  1:4291584/150‑1 (MQ=255)
                                                          cGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACt                                                                                             >  1:4214515/1‑150 (MQ=255)
                                                            ggCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTc                                                                                           <  1:976646/150‑1 (MQ=255)
                                                              cagcagCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAg                                                                                         <  2:685365/150‑1 (MQ=255)
                                                                gcagcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCt                                                                                       >  1:783949/1‑150 (MQ=255)
                                                                  agcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTgg                                                                                     <  2:2560788/150‑1 (MQ=255)
                                                                   gcCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTggg                                                                                    <  2:3302771/150‑1 (MQ=255)
                                                                     cTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGt                                                                                  <  2:2492750/150‑1 (MQ=255)
                                                                       ggCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTac                                                                                <  2:2114322/150‑1 (MQ=255)
                                                                           gggCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTaccacc                                                                            <  1:1022960/150‑1 (MQ=255)
                                                                              cTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGag                                                                         >  1:62147/1‑150 (MQ=255)
                                                                                ggACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagca                                                                       <  2:580189/150‑1 (MQ=255)
                                                                                ggACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagca                                                                       <  1:4438653/150‑1 (MQ=255)
                                                                                 gACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGagcag                                                                      >  1:1518469/1‑150 (MQ=255)
                                                                                      ccATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCaa                                                                 <  1:2480182/150‑1 (MQ=255)
                                                                                       cctCCGGTGGCGCAGCTCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAg                                                                <  2:3132463/148‑1 (MQ=255)
                                                                                       cATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAg                                                                >  2:4049479/1‑150 (MQ=255)
                                                                                       cATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAg                                                                >  2:958296/1‑150 (MQ=255)
                                                                                       cATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAg                                                                >  1:2181582/1‑150 (MQ=255)
                                                                                        ctCCGGTGGTGCAGATCGCCCTGGACGAGGCCCAGGCCGATGCCCGCTGGGCAGGCCGCGAACGGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAgg                                                               <  2:3560889/149‑1 (MQ=255)
                                                                                            ggtggtGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCACGGGCaa                                                           >  2:3639838/1‑150 (MQ=255)
                                                                                             gtggtgCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAg                                                          >  2:3076560/1‑150 (MQ=255)
                                                                                                  gCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCAt                                                     <  1:1489187/150‑1 (MQ=255)
                                                                                                  gCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCAt                                                     >  1:3052321/1‑150 (MQ=255)
                                                                                                  gCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCAt                                                     >  1:2775881/1‑150 (MQ=255)
                                                                                                    agagCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATgg                                                   >  2:4510697/1‑150 (MQ=255)
                                                                                                      aTCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGcc                                                 >  1:1750250/1‑150 (MQ=255)
                                                                                                      aTCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCCGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGcc                                                 >  2:2477308/1‑150 (MQ=255)
                                                                                                       tCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCa                                                >  2:2677973/1‑150 (MQ=255)
                                                                                                       tCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCa                                                >  2:2678145/1‑150 (MQ=255)
                                                                                                       tCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCa                                                >  2:3610428/1‑150 (MQ=255)
                                                                                                         gCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAAt                                              >  2:696723/1‑150 (MQ=255)
                                                                                                             tGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCt                                          >  2:1835687/1‑150 (MQ=255)
                                                                                                              ggAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTg                                         >  1:1637992/1‑150 (MQ=255)
                                                                                                               gAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCgcggcg                                                                                                               >  1:4141530/1‑79 (MQ=255)
                                                                                                               gAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTgg                                        <  1:4007806/150‑1 (MQ=255)
                                                                                                                atgATGCCCTGGCCTATGGCCTCTGGGCTGGCCGAGATCTGCCCTGCGAGGCGCAGCTGGAATACGCCTTGCGCGGCGGACTGACGGATGGCGTCTTCAGCTGGGGTACCACCGAGGAGCCCAAGGGCTAGCTCATGGCCAATACCTGGc                                       <  2:1750250/148‑1 (MQ=255)
                                                                                                                  gATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAg                                     >  1:2585361/1‑150 (MQ=255)
                                                                                                                  cctGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAg                                     >  2:3349916/3‑150 (MQ=255)
                                                                                                                      cccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCCCATGGCCAATACCTGGCAGGGTc                                 <  2:2722572/150‑1 (MQ=255)
                                                                                                                       ccAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCa                                >  2:1999408/1‑150 (MQ=255)
                                                                                                                         aGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGt                              >  1:3322872/1‑150 (MQ=255)
                                                                                                                          ggCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGtt                             <  1:1817218/150‑1 (MQ=255)
                                                                                                                             cTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTccc                          <  2:3967550/150‑1 (MQ=255)
                                                                                                                               aTGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCtt                        >  2:2168246/1‑150 (MQ=255)
                                                                                                                                 gCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTAt                      <  2:4238938/150‑1 (MQ=255)
                                                                                                                                   ccGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCg                    >  1:2646770/1‑150 (MQ=255)
                                                                                                                                      cTGGGCAGTTCGCGGCCTGCCCAGCGGGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCaa                 >  2:1785677/1‑150 (MQ=255)
                                                                                                                                       tGGGCAGTTCGCGGCCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAAt                >  1:2383644/1‑150 (MQ=255)
                                                                                                                                           cAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgcgg            <  1:4095466/150‑1 (MQ=255)
                                                                                                                                            aGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgcggc           >  1:1408131/1‑150 (MQ=255)
                                                                                                                                            aGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATgcggc           >  1:385738/1‑150 (MQ=255)
                                                                                                                                              gCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCAGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGa         >  2:3361178/1‑150 (MQ=255)
                                                                                                                                               ccGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGaa        >  2:1392365/1‑150 (MQ=255)
                                                                                                                                               ccGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGaa        >  2:1533721/1‑150 (MQ=255)
                                                                                                                                                cgcgAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAg       >  1:2128886/1‑150 (MQ=255)
                                                                                                                                                 gcgAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAgg      >  1:3693149/1‑150 (MQ=255)
                                                                                                                                                    aaCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATg   >  2:1376937/1‑150 (MQ=255)
                                                                                                                                                    aaCTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATg   <  2:64358/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGCTCGGCCCCGGCTGGCAGTTCGTGCCCGGCGCCAACTGGCGACACCCGCAAGGGCCGGGCAGCAGCCTGGCCGGGCTGGACAACCATCCGGTGGTGCAGATCGCCCTGGAAGATGCCCAGGCCTATGCCCGCTGGGCAGGCCGCGAACTGCCCAGCGAGGCGCAGCTGGAATACGCCATGCGCGGCGGCCTGACCGATGCCGACTTCAGCTGGGGTACCACCGAGCAGCCCAAGGGCAAGCTCATGGCCAATACCTGGCAGGGTCAGTTCCCTTATCGCAATGCGGCGAAGGATGG  >  NC_002947/3792927‑3793225

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