Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,613,037 C→T 5.5% intergenic (+31/‑2076) 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,0370CT 5.5% 560.8 / 15.3 182intergenic (+31/‑2076)hcpC‑II/PP_4084hemolysin‑coregulated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (80/91);  new base T (7/3);  total (87/95)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 201 total aligned reads displayed.

AAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCACCCTTT  >  NC_002947/4612888‑4613185
                                                                                                                                                     |                                                                                                                                                    
aaGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgccc                                                                                                                                                      >  2:1706937/1‑150 (MQ=37)
 aGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCACCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgcccg                                                                                                                                                     >  1:22719/1‑150 (MQ=37)
  gACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgcccgc                                                                                                                                                    >  2:1700337/1‑150 (MQ=37)
   aCAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgcccgcc                                                                                                                                                   >  1:2610964/1‑150 (MQ=37)
      aTGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCGGGCATGGCCTGCCCgg                                                                                                                                                <  1:713660/150‑1 (MQ=11)
      aTGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCgg                                                                                                                                                >  2:4276992/1‑150 (MQ=37)
       tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCGGGCATGGCCTGCCCggg                                                                                                                                               <  2:1731874/150‑1 (MQ=11)
       tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCggg                                                                                                                                               >  1:3194844/1‑150 (MQ=37)
       tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCggg                                                                                                                                               >  1:2421912/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)
             cTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTtgt                                                                                                                                         <  2:3004713/150‑1 (MQ=37)
              tACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTtgtg                                                                                                                                        <  1:1832510/150‑1 (MQ=37)
               acccccAGCGGGGTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTgtgt                                                                                                                                       <  2:2732913/147‑1 (MQ=255)
                  cccAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGa                                                                                                                                    <  1:4472544/150‑1 (MQ=37)
                  cccAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGa                                                                                                                                    <  1:4228275/150‑1 (MQ=37)
                   ccAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGAt                                                                                                                                   <  1:73424/150‑1 (MQ=37)
                   ccAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGAt                                                                                                                                   <  2:3995671/150‑1 (MQ=37)
                   ccAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGAt                                                                                                                                   >  1:4294435/1‑150 (MQ=37)
                    cAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGAtt                                                                                                                                  >  2:4212204/1‑150 (MQ=37)
                           ttGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGAcgg                                                                                                                           >  1:2812154/1‑150 (MQ=34)
                            tGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGAcggc                                                                                                                          >  2:3448822/1‑150 (MQ=34)
                                    cGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAgg                                                                                                                  <  1:2790780/150‑1 (MQ=34)
                                       tGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATc                                                                                                               <  1:3463148/150‑1 (MQ=34)
                                         gCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTGTGTGTTGATTCAGACGGCGGACAAGGATCtt                                                                                                             >  1:1750745/1‑150 (MQ=21)
                                            tATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTa                                                                                                          <  1:3318681/150‑1 (MQ=34)
                                               cGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGc                                                                                                       >  1:2840158/1‑150 (MQ=34)
                                                 aaaGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGcccgcccg                                                                                                                                                 <  2:3580192/106‑1 (MQ=35)
                                                 aaaGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGcccgcccg                                                                                                                                                 >  1:3580192/1‑106 (MQ=35)
                                                 aaaGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCt                                                                                                     <  1:3154499/150‑1 (MQ=34)
                                                 aaaGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCt                                                                                                     >  1:4377507/1‑150 (MQ=34)
                                                  aaGATCACGTGGGCTCACGTTGTAGCCAATACCGAAGGCTCTGTCGTCTGGCGCGCACCCGCCGCCTGCCCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTa                                                                                                    <  1:252795/150‑1 (MQ=255)
                                                    gATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGc                                                                                                  <  1:3702317/150‑1 (MQ=34)
                                                     aTCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGcc                                                                                                 >  2:540037/1‑150 (MQ=34)
                                                     aTCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGcc                                                                                                 >  2:4087245/1‑150 (MQ=34)
                                                      tCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCa                                                                                                >  2:1156070/1‑150 (MQ=34)
                                                          gTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCt                                                                                            >  1:2329867/1‑150 (MQ=34)
                                                             ggCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGAt                                                                                         <  1:798518/150‑1 (MQ=34)
                                                              gCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGAtt                                                                                        <  1:2338979/150‑1 (MQ=34)
                                                                tCACGATGTAGCCAATACCGCAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGAGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTg                                                                                      <  2:371356/150‑1 (MQ=21)
                                                                tCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTg                                                                                      >  2:2294275/1‑150 (MQ=34)
                                                                tCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTg                                                                                      >  2:4478626/1‑150 (MQ=34)
                                                                tCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTg                                                                                      >  1:2730230/1‑150 (MQ=34)
                                                                 cACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgc                                                                                     <  2:4405062/150‑1 (MQ=34)
                                                                   cGATGTAGCCACTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAAGGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgcgc                                                                                   <  2:2226689/150‑1 (MQ=21)
                                                                   cGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgcgc                                                                                   >  2:3383444/1‑150 (MQ=34)
                                                                    gATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCt                                                                                  >  2:3250061/1‑150 (MQ=34)
                                                                        tAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAg                                                                              >  2:742021/1‑150 (MQ=34)
                                                                           ccAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCtt                                                                           >  2:2690808/1‑150 (MQ=34)
                                                                            cAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCttt                                                                          <  1:1174680/150‑1 (MQ=34)
                                                                            cAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCttt                                                                          >  2:1012590/1‑150 (MQ=34)
                                                                             aaTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTg                                                                         <  1:4063087/150‑1 (MQ=34)
                                                                             aaTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTg                                                                         <  1:1458816/150‑1 (MQ=34)
                                                                              aTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGa                                                                        <  2:1738311/150‑1 (MQ=33)
                                                                                aCCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAgg                                                                      <  1:3870392/150‑1 (MQ=32)
                                                                                  cGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCa                                                                    >  2:4356810/1‑150 (MQ=32)
                                                                                       gCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGt                                                               <  1:4332331/150‑1 (MQ=32)
                                                                                       gCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGt                                                               <  1:2313264/150‑1 (MQ=32)
                                                                                         tctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTg                                                             <  2:3399549/150‑1 (MQ=32)
                                                                                         tctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTg                                                             <  2:3338074/150‑1 (MQ=32)
                                                                                         tctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTg                                                             >  2:3056868/1‑150 (MQ=32)
                                                                                         tctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTg                                                             <  1:16691/150‑1 (MQ=32)
                                                                                         tctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTGTGAGGCACTCGTTg                                                             >  1:1529471/1‑150 (MQ=17)
                                                                                         tctGACGACTGGAGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGACATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTg                                                             >  1:2654745/1‑150 (MQ=17)
                                                                                           tGACGACTGGCGCGTACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGc                                                           >  2:4159050/1‑150 (MQ=17)
                                                                                                acTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGc                                                      <  1:1040950/150‑1 (MQ=32)
                                                                                                  tGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCttt                                                                                                            >  2:1945094/1‑94 (MQ=33)
                                                                                                  tGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCttt                                                                                                            <  1:1945094/94‑1 (MQ=33)
                                                                                                    gcgcgcACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGaa                                                  <  1:1137114/150‑1 (MQ=33)
                                                                                                    gcgcgcACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGaa                                                  <  2:2868622/150‑1 (MQ=33)
                                                                                                       cgcACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTc                                               <  2:2648381/150‑1 (MQ=34)
                                                                                                       cgcACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTc                                               <  1:954578/150‑1 (MQ=34)
                                                                                                        gcACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCg                                              >  1:4142536/1‑150 (MQ=34)
                                                                                                        gcACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCg                                              <  2:1783880/150‑1 (MQ=34)
                                                                                                          aCCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGaa                                            >  1:543751/1‑150 (MQ=34)
                                                                                                          aCCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGaa                                            <  2:4445705/150‑1 (MQ=34)
                                                                                                          aCCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGaa                                            <  1:764536/150‑1 (MQ=34)
                                                                                                            ccgccgCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAgt                                          >  1:3251392/1‑150 (MQ=34)
                                                                                                                   cTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCg                                   >  2:3918537/1‑150 (MQ=34)
                                                                                                                     gACCCGTGCAAATGCCCCGGGCAGGCCATGCCCACCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATTCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGtt                                 <  1:1706937/150‑1 (MQ=21)
                                                                                                                         cGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGaa                             <  1:2501126/150‑1 (MQ=32)
                                                                                                                            gCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTc                          >  1:238387/1‑150 (MQ=32)
                                                                                                                             cAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGCGTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGGTCGAACTcg                         >  1:3098139/1‑149 (MQ=12)
                                                                                                                             cAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTcc                         <  1:3455770/150‑1 (MQ=32)
                                                                                                                              acaaagcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCcg                        >  2:648843/7‑150 (MQ=11)
                                                                                                                               caaagcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCcgc                       >  1:3993646/6‑150 (MQ=11)
                                                                                                                               aaTGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCcgc                       <  2:1518903/150‑1 (MQ=32)
                                                                                                                                  agcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAg                    >  1:1996267/3‑150 (MQ=11)
                                                                                                                                    cccGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAggg                  >  2:2843051/1‑150 (MQ=32)
                                                                                                                                     ccGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGc                 >  2:319228/1‑150 (MQ=32)
                                                                                                                                           aGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAgc           <  2:3825996/150‑1 (MQ=32)
                                                                                                                                             gCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAgcgc         <  2:2611012/150‑1 (MQ=32)
                                                                                                                                             gCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAgcgc         >  2:458442/1‑150 (MQ=32)
                                                                                                                                             gCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAgcgc         >  2:1182786/1‑150 (MQ=32)
                                                                                                                                                aTGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCAcc      >  2:2005363/1‑150 (MQ=32)
                                                                                                                                                    ccgcccgGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCACCCttt  <  2:2110464/150‑1 (MQ=32)
                                                                                                                                                    ccgcccgGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCACCCttt  >  2:1838674/1‑150 (MQ=32)
                                                                                                                                                    ccgcccgGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCACCCttt  >  1:1381639/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                    
AAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCACCCTTT  >  NC_002947/4612888‑4613185

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: