Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,613,032 A→T 5.4% intergenic (+26/‑2081) 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,0320AT 5.4% 640.7 / 17.7 198intergenic (+26/‑2081)hcpC‑II/PP_4084hemolysin‑coregulated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (92/95);  new base T (7/4);  total (99/99)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

CGACAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCACC  >  NC_002947/4612883‑4613181
                                                                                                                                                     |                                                                                                                                                     
  aCAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATg                                                                                                                                                     <  1:3863145/150‑1 (MQ=37)
  aCAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATg                                                                                                                                                     >  1:1855312/1‑150 (MQ=37)
       gACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgcccgc                                                                                                                                                >  2:1700340/1‑150 (MQ=37)
       gACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgcccgc                                                                                                                                                <  1:733226/150‑1 (MQ=37)
        aCAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgcccgcc                                                                                                                                               >  1:2610964/1‑150 (MQ=37)
         cAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGCGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgcccgccc                                                                                                                                              <  1:3219688/150‑1 (MQ=37)
          aaTGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGcccgcccg                                                                                                                                             >  1:2868622/1‑150 (MQ=37)
           aTGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCgg                                                                                                                                            >  2:4276992/1‑150 (MQ=37)
           aTGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCgg                                                                                                                                            >  2:1815370/1‑150 (MQ=37)
            tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCGGGCATGGCCTGCCCggg                                                                                                                                           <  1:439755/150‑1 (MQ=11)
            tGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCggg                                                                                                                                           <  1:4016248/150‑1 (MQ=37)
             gAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGc                                                                                                                                          >  2:720361/1‑150 (MQ=37)
             gAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGc                                                                                                                                          >  1:3812917/1‑150 (MQ=37)
              aaaaCTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCt                                                                                                                                         <  1:4286412/150‑1 (MQ=37)
               aaaCTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCtt                                                                                                                                        <  2:1735017/150‑1 (MQ=37)
                aaCTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCatt                                                                                                                                       <  1:2329339/150‑1 (MQ=37)
                 aCTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTtg                                                                                                                                      >  2:3541003/1‑150 (MQ=37)
                 aCTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTtg                                                                                                                                      <  2:2698670/150‑1 (MQ=37)
                    acacCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTgtgt                                                                                                                                   <  2:1767920/150‑1 (MQ=37)
                        ccAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGAt                                                                                                                               <  2:3995671/150‑1 (MQ=37)
                         cAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGAtt                                                                                                                              >  2:4212204/1‑150 (MQ=37)
                          aGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTc                                                                                                                             <  1:2277909/150‑1 (MQ=35)
                              ggTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGAc                                                                                                                         >  1:983780/1‑150 (MQ=34)
                                 tGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGAcggc                                                                                                                      >  1:2785719/1‑150 (MQ=34)
                                 tGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGAcggc                                                                                                                      <  1:2408063/150‑1 (MQ=34)
                                        tCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAg                                                                                                               >  1:2810065/1‑150 (MQ=34)
                                        tCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAg                                                                                                               >  1:3406065/1‑150 (MQ=34)
                                            tGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATc                                                                                                           <  1:953451/150‑1 (MQ=34)
                                            tGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATc                                                                                                           <  1:3463148/150‑1 (MQ=34)
                                              gCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCtt                                                                                                         >  2:2754561/1‑150 (MQ=34)
                                               cGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCttt                                                                                                        <  1:3176273/150‑1 (MQ=34)
                                                gTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCtttt                                                                                                       >  1:1717531/1‑150 (MQ=34)
                                                gTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTGTGTGTTGATTCAGACGGCGGACAAGGATCtttt                                                                                                       >  2:3055862/1‑150 (MQ=21)
                                                     gAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGcc                                                                                                  >  1:1592491/1‑150 (MQ=34)
                                                     gAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGcc                                                                                                  >  1:772263/1‑150 (MQ=34)
                                                      aaaGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGcccgcccg                                                                                                                                             >  1:3580192/1‑106 (MQ=35)
                                                      aaaGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCt                                                                                                 <  1:3154413/150‑1 (MQ=34)
                                                       aaGATCACGTGGGCTCACGTTGTAGCCAATACCGAAGGCTCTGTCGTCTGGCGCGCACCCGCCGCCTGCCCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTa                                                                                                <  1:252795/150‑1 (MQ=255)
                                                        aGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAg                                                                                               >  1:2696889/1‑150 (MQ=34)
                                                         gATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGc                                                                                              <  1:3702317/150‑1 (MQ=34)
                                                          aTCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGcc                                                                                             >  2:4087245/1‑150 (MQ=34)
                                                            cACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTGTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCaa                                                                                           >  1:1503712/1‑150 (MQ=21)
                                                                 gggCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGa                                                                                      <  2:2890628/150‑1 (MQ=34)
                                                                  ggCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGAt                                                                                     <  1:4144377/150‑1 (MQ=34)
                                                                  ggCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCATTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGAt                                                                                     <  2:2665845/150‑1 (MQ=21)
                                                                   gCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGAtt                                                                                    <  2:3414228/150‑1 (MQ=34)
                                                                     tCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTg                                                                                  >  2:2294275/1‑150 (MQ=34)
                                                                     tCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTg                                                                                  >  2:4478626/1‑150 (MQ=34)
                                                                      cACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgc                                                                                 <  2:4405062/150‑1 (MQ=34)
                                                                       aCGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgcg                                                                                >  2:3155845/1‑150 (MQ=34)
                                                                       aCGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgcg                                                                                >  2:1179623/1‑150 (MQ=34)
                                                                        cGATGTAGCCACTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAAGGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgcgc                                                                               <  2:2226689/150‑1 (MQ=21)
                                                                        cGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTgcgc                                                                               <  1:2368893/150‑1 (MQ=34)
                                                                          aTGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTa                                                                             <  2:1139602/150‑1 (MQ=34)
                                                                             tAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAg                                                                          >  2:742021/1‑150 (MQ=34)
                                                                             tAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCACCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAg                                                                          >  1:454560/1‑150 (MQ=21)
                                                                              aGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGc                                                                         >  1:2335762/1‑150 (MQ=34)
                                                                                 cAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCttt                                                                      <  1:1174680/150‑1 (MQ=34)
                                                                                  aaTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTg                                                                     >  2:547407/1‑150 (MQ=34)
                                                                                   aTACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGa                                                                    <  2:3392143/150‑1 (MQ=33)
                                                                                    tACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGCCAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAg                                                                   <  1:659820/150‑1 (MQ=17)
                                                                                    tACCGAAGGCTCTGAAGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCAGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAg                                                                   <  1:3974474/150‑1 (MQ=17)
                                                                                       cGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgc                                                                                                                                                    >  1:2636058/1‑66 (MQ=35)
                                                                                       cGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATgc                                                                                                                                                    <  2:2636058/66‑1 (MQ=35)
                                                                                       cGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCa                                                                >  2:1583571/1‑150 (MQ=32)
                                                                                       cGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCa                                                                >  1:4331919/1‑150 (MQ=32)
                                                                                       cGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCa                                                                >  2:4356810/1‑150 (MQ=32)
                                                                                             ctctGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGtt                                                          <  2:3348413/150‑1 (MQ=32)
                                                                                                tGACGACTGGCGCGTACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGc                                                       >  2:4159050/1‑150 (MQ=17)
                                                                                                   cgacTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCaaa                                                    <  2:1714242/150‑1 (MQ=32)
                                                                                                     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                                              >  1:2491907/1‑150 (MQ=33)
                                                                                                          cgcgcACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAAt                                             <  2:3850705/150‑1 (MQ=34)
                                                                                                           gcgcACCCGCCGCCTGACCCGTGCACATGCCCCGGGCCGGCCATGCCCGCCCGGGCTGTGTGTTGATTCAGACGGCGGACAAGGATCTTTCATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAAtt                                            >  2:3808414/1‑150 (MQ=255)
                                                                                                               aCCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGaa                                        <  2:4445705/150‑1 (MQ=34)
                                                                                                                    ccgccTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAgtggt                                   <  1:4075186/150‑1 (MQ=34)
                                                                                                                     cgccTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAgtggtg                                  >  2:2859033/1‑150 (MQ=34)
                                                                                                                        cTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCg                               >  2:3918537/1‑150 (MQ=34)
                                                                                                                           aCCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTc                            <  1:3262293/150‑1 (MQ=34)
                                                                                                                                  cAAATGCCCCGGGCAGGCCTTGCCCGCCCGGGCTTTGTGTTGATTCAGCCGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTcc                     <  2:3232832/150‑1 (MQ=17)
                                                                                                                                  cAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTcc                     <  1:3455770/150‑1 (MQ=32)
                                                                                                                                   acaaagcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCcg                    >  2:648843/7‑150 (MQ=11)
                                                                                                                                   aaaTGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCcg                    >  2:3253596/1‑150 (MQ=32)
                                                                                                                                    caatgcccGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCcgc                   >  1:2110464/6‑150 (MQ=17)
                                                                                                                                    aaTGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCcgc                   <  2:133541/150‑1 (MQ=32)
                                                                                                                                     aaagcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCcgcg                  >  1:1391461/5‑150 (MQ=11)
                                                                                                                                      tGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGa                 <  1:1689569/150‑1 (MQ=32)
                                                                                                                                       gCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAg                <  1:435609/150‑1 (MQ=32)
                                                                                                                                       agcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAg                >  1:2819565/3‑150 (MQ=11)
                                                                                                                                        gcccGGGCGGGCATGGCCTGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAgg               >  2:1468576/2‑150 (MQ=11)
                                                                                                                                         cccGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAggg              >  1:743058/1‑150 (MQ=32)
                                                                                                                                                aGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAgc       >  2:3621879/1‑150 (MQ=32)
                                                                                                                                                  gCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAgcgc     <  2:2611012/150‑1 (MQ=32)
                                                                                                                                                  gCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAgcgc     >  2:458442/1‑150 (MQ=32)
                                                                                                                                                    cATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCAc   >  2:1033858/1‑150 (MQ=32)
                                                                                                                                                     aTGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCAcc  <  1:1772621/150‑1 (MQ=32)
                                                                                                                                                     aTGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCAcc  >  2:2005363/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                     
CGACAAAGACAATGAAAACTACACCCAGCTGGTTGAAGTATCGCTGGCGTATCGAAAGATCACGTGGGCTCACGATGTAGCCAATACCGAAGGCTCTGACGACTGGCGCGCACCCGCCGCCTGACCCGTGCAAATGCCCCGGGCAGGCCATGCCCGCCCGGGCTTTGTGTTGATTCAGACGGCGGACAAGGATCTTTTATGCCTAGCCAATCTGATTTGCGCTACAGCTTTGAGGCACTCGTTGGCAAAGCTGAATTCGAAGTGGTGTCGTTCGAACTCCGCGAGGGCATCAGCGCACC  >  NC_002947/4612883‑4613181

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: