Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,613,027 A→G 5.1% intergenic (+21/‑2086) hcpC‑II → / → PP_4084 hemolysin‑coregulated protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,613,0270AG 5.1% 616.1 / 14.5 198intergenic (+21/‑2086)hcpC‑II/PP_4084hemolysin‑coregulated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (91/95);  new base G (7/3);  total (100/98)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 200 total aligned reads displayed.

CACTCGACAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGC  >  NC_002947/4612879‑4613176
                                                                                                                                                    |                                                                                                                                                     
   tCGACAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGcc                                                                                                                                                   >  1:2995177/1‑150 (MQ=37)
   tCGACAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGcc                                                                                                                                                   >  1:1882104/1‑150 (MQ=37)
    cGACAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCa                                                                                                                                                  >  2:3870392/1‑150 (MQ=37)
      aCAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATg                                                                                                                                                <  1:3863145/150‑1 (MQ=37)
      aCAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATg                                                                                                                                                >  1:1855312/1‑150 (MQ=37)
         aaGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgccc                                                                                                                                             >  2:1706937/1‑150 (MQ=37)
          aGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCACCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgcccg                                                                                                                                            >  1:22719/1‑150 (MQ=37)
           gACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgcccgc                                                                                                                                           <  1:733226/150‑1 (MQ=37)
            aCAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgcccgcc                                                                                                                                          >  1:2611012/1‑150 (MQ=37)
                tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCGGGCATGGCCTGCCCggg                                                                                                                                      <  1:439755/150‑1 (MQ=11)
                tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCGGGCATGGCCTGCCCggg                                                                                                                                      <  2:1731874/150‑1 (MQ=11)
                tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCggg                                                                                                                                      <  2:849318/150‑1 (MQ=37)
                tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCggg                                                                                                                                      <  1:4016248/150‑1 (MQ=37)
                tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCggg                                                                                                                                      >  1:3194844/1‑150 (MQ=37)
                tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCggg                                                                                                                                      >  1:2421912/1‑150 (MQ=37)
                 gAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGc                                                                                                                                     >  1:1209677/1‑150 (MQ=37)
                 gAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGc                                                                                                                                     >  1:3812917/1‑150 (MQ=37)
                 gAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGc                                                                                                                                     >  2:720361/1‑150 (MQ=37)
                  aaatCTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCt                                                                                                                                    >  2:1070312/1‑150 (MQ=37)
                  aaaaCTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCt                                                                                                                                    >  2:1451356/1‑150 (MQ=37)
                  aaaaCTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCt                                                                                                                                    <  1:4286412/150‑1 (MQ=37)
                   aaaCTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCtt                                                                                                                                   <  2:4436360/150‑1 (MQ=37)
                     aCTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTtg                                                                                                                                 <  2:2698670/150‑1 (MQ=37)
                      cTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTtgt                                                                                                                                <  2:3004713/150‑1 (MQ=37)
                                cTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCGAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAg                                                                                                                      >  1:1155523/1‑150 (MQ=21)
                                  ggTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGAc                                                                                                                    >  1:983780/1‑150 (MQ=34)
                                    ttGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGAcgg                                                                                                                  >  1:2812154/1‑150 (MQ=34)
                                     tGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGAcggc                                                                                                                 <  2:4283964/150‑1 (MQ=34)
                                     tGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGAcggc                                                                                                                 >  2:3448822/1‑150 (MQ=34)
                                            tCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAg                                                                                                          >  1:3406065/1‑150 (MQ=34)
                                             cGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAgg                                                                                                         <  1:2790780/150‑1 (MQ=34)
                                                tGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATc                                                                                                      <  1:953451/150‑1 (MQ=34)
                                                tGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATc                                                                                                      <  1:3463148/150‑1 (MQ=34)
                                                  gCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCtt                                                                                                    >  2:2754561/1‑150 (MQ=34)
                                                  gCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTGTGTGTTGATTCAGACGGCGGACAAGGATCtt                                                                                                    >  1:1750745/1‑150 (MQ=21)
                                                   cGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCttt                                                                                                   <  1:3176273/150‑1 (MQ=34)
                                                        cGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGc                                                                                              >  1:2840158/1‑150 (MQ=34)
                                                         gAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGcc                                                                                             >  1:772263/1‑150 (MQ=34)
                                                          aaaGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGcccgcccg                                                                                                                                        >  1:3580192/1‑106 (MQ=35)
                                                          aaaGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGcccgcccg                                                                                                                                        <  2:3580192/106‑1 (MQ=35)
                                                          aaaGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCt                                                                                            <  1:3154413/150‑1 (MQ=34)
                                                             gATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGc                                                                                         <  1:3702317/150‑1 (MQ=34)
                                                                cACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCaa                                                                                      <  2:1887832/150‑1 (MQ=34)
                                                                  cGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATc                                                                                    <  2:4396430/150‑1 (MQ=34)
                                                                   gTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCt                                                                                   >  1:2329867/1‑150 (MQ=34)
                                                                      ggCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGAt                                                                                <  1:4144377/150‑1 (MQ=34)
                                                                      ggCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCATTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGAt                                                                                <  2:2665845/150‑1 (MQ=21)
                                                                       gCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGAtt                                                                               <  1:2338979/150‑1 (MQ=34)
                                                                         tCACGATGTAGCCAATACCGCAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGAGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTg                                                                             <  2:371356/150‑1 (MQ=21)
                                                                         tCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCCGGCCATGCACGCCCGGGCGGTGTGGTGATTCAGACGGCGGACCAGGATCTTGTGTGCCTAGCCAATCTGATTTg                                                                             >  2:4478663/1‑150 (MQ=11)
                                                                         tCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTg                                                                             >  2:2294275/1‑150 (MQ=34)
                                                                         tCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTg                                                                             >  1:2730230/1‑150 (MQ=34)
                                                                         tCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTg                                                                             >  1:2590760/1‑150 (MQ=34)
                                                                           aCGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgcg                                                                           >  2:1179623/1‑150 (MQ=34)
                                                                            cGATGTAGCCACTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAAGGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgcgc                                                                          <  2:2226689/150‑1 (MQ=21)
                                                                            cGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgcgc                                                                          >  2:3383444/1‑150 (MQ=34)
                                                                             gATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCt                                                                         >  2:3250061/1‑150 (MQ=34)
                                                                              aTGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTa                                                                        <  2:1139602/150‑1 (MQ=34)
                                                                                 tAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCACCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAg                                                                     >  1:454560/1‑150 (MQ=21)
                                                                                      aaTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTg                                                                >  2:547407/1‑150 (MQ=34)
                                                                                      aaTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTg                                                                <  1:1458816/150‑1 (MQ=34)
                                                                                       aTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGa                                                               <  2:1738311/150‑1 (MQ=33)
                                                                                       aTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGa                                                               <  2:4017647/150‑1 (MQ=33)
                                                                                        tACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGCCAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAg                                                              <  1:659820/150‑1 (MQ=17)
                                                                                        tACCGAAGGCTCTGAAGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCAGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAg                                                              <  1:3974474/150‑1 (MQ=17)
                                                                                           cGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgc                                                                                                                                               <  2:2636058/66‑1 (MQ=35)
                                                                                           cGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCa                                                           >  2:1583571/1‑150 (MQ=32)
                                                                                           cGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCa                                                           >  1:4331919/1‑150 (MQ=32)
                                                                                                gCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGt                                                      <  1:2313264/150‑1 (MQ=32)
                                                                                                  tctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTg                                                    <  2:3338074/150‑1 (MQ=32)
                                                                                                  tctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTg                                                    <  2:3399549/150‑1 (MQ=32)
                                                                                                  tctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTg                                                    <  1:16691/150‑1 (MQ=32)
                                                                                                  tctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTg                                                    >  2:3056868/1‑150 (MQ=32)
                                                                                                  tctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTGTGAGGCACTCGTTg                                                    >  1:1529471/1‑150 (MQ=17)
                                                                                                  tctGACGACTGGAGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGACATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTg                                                    >  1:2654745/1‑150 (MQ=17)
                                                                                                        gacTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAg                                              <  2:1981584/150‑1 (MQ=32)
                                                                                                         acTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGc                                             <  1:1040950/150‑1 (MQ=32)
                                                                                                           tGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCGTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAATGCTg                                           >  2:4233453/1‑150 (MQ=17)
                                                                                                             gcgcgcACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGaa                                         <  2:2868622/150‑1 (MQ=33)
                                                                                                                 gcACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCg                                     <  2:1783880/150‑1 (MQ=34)
                                                                                                                        ccgccTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAgtggt                              <  1:4075186/150‑1 (MQ=34)
                                                                                                                            cTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCg                          >  2:3918537/1‑150 (MQ=34)
                                                                                                                              gACCCGTGCAAATGCCCCGGGCAGGCCATGCCCACCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATTCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGtt                        <  1:1706937/150‑1 (MQ=21)
                                                                                                                               aCCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTc                       <  1:3262293/150‑1 (MQ=34)
                                                                                                                                  cGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGaa                    <  2:729753/150‑1 (MQ=32)
                                                                                                                                  cGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGaa                    <  1:2501126/150‑1 (MQ=32)
                                                                                                                                  cGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGaa                    <  2:419449/150‑1 (MQ=32)
                                                                                                                                    tGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACt                  >  2:4449025/1‑150 (MQ=32)
                                                                                                                                     gCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTc                 >  1:238387/1‑150 (MQ=32)
                                                                                                                                      cAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTcc                <  2:4142536/150‑1 (MQ=32)
                                                                                                                                      cAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTcc                <  1:3455770/150‑1 (MQ=32)
                                                                                                                                       acaaagcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCcg               >  2:648843/7‑150 (MQ=11)
                                                                                                                                        caaagcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCcgc              >  1:1099238/6‑150 (MQ=11)
                                                                                                                                           gCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAg           <  1:435609/150‑1 (MQ=32)
                                                                                                                                           agcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAg           >  1:1996267/3‑150 (MQ=11)
                                                                                                                                            gcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAgg          >  2:1468576/2‑150 (MQ=11)
                                                                                                                                             cccGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAggg         >  2:2843051/1‑150 (MQ=32)
                                                                                                                                             cccGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAggg         >  1:743058/1‑150 (MQ=32)
                                                                                                                                              ccGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGc        >  2:319228/1‑150 (MQ=32)
                                                                                                                                                    aGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAgc  <  2:3825996/150‑1 (MQ=32)
                                                                                                                                                    |                                                                                                                                                     
CACTCGACAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGC  >  NC_002947/4612879‑4613176

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: