Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,758,669 A→G 6.0% intergenic (+368/‑320) glnG → / → PP_5050 two‑component system DNA‑binding response regulator GlnL/GlnG/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,758,6690AG 6.0% 442.8 / 4.3 160intergenic (+368/‑320)glnG/PP_5050two‑component system DNA‑binding response regulator GlnL/GlnG/transposase
Reads supporting (aligned to +/- strand):  ref base A (68/75);  new base G (8/2);  total (78/82)
Fisher's exact test for biased strand distribution p-value = 5.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

GCCTACCACGGCGGTGGGCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCC  >  NC_002947/5758520‑5758818
                                                                                                                                                      |                                                                                                                                                     
  cttccATGGCGGTGGGCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGACGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGaca                                                                                                                                                     <  1:1007494/147‑1 (MQ=255)
  cTACCACGGCGGTGGGCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGGTGaca                                                                                                                                                     >  2:193520/1‑150 (MQ=255)
   ttccccggcggTGGTCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGacac                                                                                                                                                    <  2:3734400/145‑1 (MQ=255)
   tACCTCGGCGGTGGTCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGacac                                                                                                                                                    <  2:2282203/150‑1 (MQ=255)
   gtcttacgcggTGGCCACATATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGacac                                                                                                                                                    <  2:656523/143‑1 (MQ=255)
    aCCACGGCGGTGGGCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTAACACGGCGTGGacacc                                                                                                                                                   >  2:1878771/1‑149 (MQ=255)
      cACGGCGGTGGGCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCCCTCCAACCACCGCGTGGAcacaca                                                                                                                                                 >  1:574104/1‑150 (MQ=255)
        cggcggTGGGCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATc                                                                                                                                               <  1:1276985/150‑1 (MQ=255)
           cggTGGGCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATcgcg                                                                                                                                            <  1:4152236/150‑1 (MQ=255)
           cggTGGGCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATcgcg                                                                                                                                            <  1:2061013/150‑1 (MQ=255)
            ggTGGCCACAGATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCgg                                                                                                                                           <  2:530497/150‑1 (MQ=255)
              tGGGCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGgtc                                                                                                                                         <  1:193462/150‑1 (MQ=255)
                 gCACATATCGGGGTAGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGgtcgtc                                                                                                                                      <  2:2457794/150‑1 (MQ=255)
                      tatCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTAACACCGCGTTGACACACATCGCGGTCGTCGTAgg                                                                                                                                 >  1:3727708/1‑150 (MQ=255)
                         ccgggTCGTCGTAGGAGCGTGTTCACCCGCGAAGAGGCCGCGCCTTCTCGCCAATCTCAACCTGACTTTCAGCGTTGCCGCCTTCACATTTTTCTCTGGTAATCCCGCTCCTAACACGGCGTGGACACACATCGCGGTCGTCGTAGGAGc                                                                                                                              <  2:3043237/148‑1 (MQ=255)
                         cGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGc                                                                                                                              <  1:3371294/150‑1 (MQ=255)
                             gtcgtcGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGt                                                                                                                          <  2:2575774/150‑1 (MQ=255)
                                    tAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTGCTACCCCGGCGTTGACACAAATGGGGGTGGTCCTAGGTGGGGGTTGACCcg                                                                                                                   >  2:918865/1‑150 (MQ=14)
                                     aGGAGCGGGTTCACCCGCGTAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCcgc                                                                                                                  <  1:2999279/150‑1 (MQ=255)
                                        aGCGGGTTCACCCGCGACGAGGCCGCGCCGGCGGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACCCATCGCTGTCGTCGTAGGAGCGGGTTCACCCGCGaa                                                                                                               <  2:183369/150‑1 (MQ=17)
                                        aGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGaa                                                                                                               <  1:439126/150‑1 (MQ=255)
                                          cGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCCCCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAaga                                                                                                             <  2:1714691/150‑1 (MQ=255)
                                                cACCCGCGACGAGGCCGCGCCTGCTGGCCAATCTCCACCGGCCGTTCAGCGTTGCGGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCgg                                                                                                       <  2:2965824/150‑1 (MQ=255)
                                                  cccGCGACGAGGCCGCGCCTGGTGGCCAACCTCCCCCTGCCGTTCAGGGTTGCCGCCGCCACGGTATGCGCGGTTAACCCCGCTCCTCCCCCGGCGTGGACACACATCGGGGTCGTCGTAGGAGCGGGTTCACCGGCGAAGAGGCCGGGc                                                                                                     <  1:186830/150‑1 (MQ=17)
                                                  cccGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACGCATCGCGGTCGTCGGAGGAGCGGGTTCACCCGCGAAGAGGCCGGGc                                                                                                     >  2:3376680/1‑150 (MQ=255)
                                                   ccGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCCCGGCGGGGCCACACAGCGCGGGCGGCGGAGGAGCGGGGTCCCCCGCGGAGAGGCCGGGcc                                                                                                    >  2:158260/1‑150 (MQ=17)
                                                   ccGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGcc                                                                                                    >  2:1413977/1‑150 (MQ=255)
                                                   ccGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGGTCACCCGCGCAGAGGCCGGGcc                                                                                                    >  2:2865234/1‑150 (MQ=255)
                                                        aaGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCctgctg                                                                                               <  2:4278929/150‑1 (MQ=255)
                                                          gagGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCCCGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCCGGCCAGCTGGc                                                                                             >  1:2148911/1‑150 (MQ=255)
                                                          gagGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCCCGGCGGGGGCACACATCGCGGGCGGCGTAGGAGCCGGTTCACCCGGGGAGAGGGCCGGCCAGCAGGc                                                                                             >  1:607588/1‑150 (MQ=14)
                                                             gCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTCAACCCGCTCCTACCACGGCGGGGGCACCCAGCGCGGGCGGCGTAGGTGCGGGGTCACCCGCGCAGAGGCCCGGCCGGCTGGCCaa                                                                                          >  2:946146/1‑150 (MQ=14)
                                                             gCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACCGCGGGGGCCCCCATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGCAGAGGCCGGGCCAGCAGGCCaa                                                                                          >  1:3274645/1‑150 (MQ=14)
                                                                  gcCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTcc                                                                                     <  1:2590390/150‑1 (MQ=255)
                                                                          gCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACCCATCGCGGTCGGCGTAGGAGCGGGTTCACCCGCGGAGAGGCCGGGCCGGCTGGCCAAACTCCACCTGCCg                                                                             >  1:3635892/1‑150 (MQ=17)
                                                                          gCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCAAGGCGTGGACACACATGGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCg                                                                             >  2:3755355/1‑150 (MQ=255)
                                                                            cAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGGTCACCCGCGAAGAGGCCGGGCCAGCTGGCCAATCACCACCTGCCGtt                                                                           >  2:615676/1‑150 (MQ=255)
                                                                               tctcCACCTGCCGTTCTGCGTTGCCGCCGCCAGCGTTTTCGCGGGTAAACCCGCTCCTACCACGGGGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAg                                                                        <  2:3846583/150‑1 (MQ=255)
                                                                                           gTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGccac                                                            <  2:2031454/150‑1 (MQ=255)
                                                                                           gTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGccac                                                            <  2:759103/150‑1 (MQ=255)
                                                                                           gTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTAACACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTAACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTACACCAGCCCTTAAGCGTTGCCGccac                                                            >  1:1388338/1‑150 (MQ=17)
                                                                                            ttCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGGCGTCGTAGGTGCGGGTTCACCCGCGAAGAGGCCGGGTCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGccacc                                                           >  2:796902/1‑150 (MQ=255)
                                                                                            ttCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCCTAACACGGCGTGGGCACCCAGCGCGGTCGGCGGAGGCGCGGGTTCACCCGCCCAGAGGCCGGGCCGGCTGGCCCAAATCCACCGGCCGGTCAGCGGTGCCGccacc                                                           >  1:1360463/1‑150 (MQ=14)
                                                                                             tCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGccacca                                                          <  2:2977710/150‑1 (MQ=255)
                                                                                                  gTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGGCACACCTCGCGGGCGTCGTTGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCGGCTGGCCAATCTCCACCGGCCGGTCAGCGTTGCCGCCACCACCgga                                                     >  1:975303/1‑148 (MQ=14)
                                                                                                  gTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTa                                                     <  2:2814936/150‑1 (MQ=255)
                                                                                                    tttacccgccACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTAtt                                                   >  2:3244102/6‑150 (MQ=255)
                                                                                                     gccgccgccACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATAGCGGTCGTCGTAGGGTCGGGTTCACCCGCCACGAGGCCGGGCCTGCTGGCCAATCTCCACCGGCCGTTCAGCGTTGCCGCCACCCCCgggtgc                                                  >  2:3240381/1‑145 (MQ=17)
                                                                                                         ccgccACCGTATTCGCGGGTAAACCCGCTCCTACCCCGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCgg                                              <  2:2914048/150‑1 (MQ=255)
                                                                                                         ccgccACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTTGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCgg                                              <  1:2886356/150‑1 (MQ=255)
                                                                                                           gccACCGTATTCGCGGGTAAAACCGCGCCTACCACCGCGTGGACACACAACGCGGTCGTGGTAGGAGCGGGTTCACCCCCCAAGAGGCCGGGCCAGCTGGCCAAAATCCACCAGCCGTTCAGCGTTGCCGCCACCACCCAAATAGCGGGt                                            >  1:224628/1‑150 (MQ=14)
                                                                                                              aCCGTATTCGCGGGTAAACCCGCTCCTACCCCGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTaaa                                         <  1:963839/150‑1 (MQ=255)
                                                                                                               ccGTATTCGGGGGTAACCCGCCTCGTACCCGGGCGTTGACACACATCGCGGTCTTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAAc                                        <  1:2305299/150‑1 (MQ=14)
                                                                                                                 gTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGACCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAAccc                                      >  2:2370017/1‑150 (MQ=255)
                                                                                                                    atcccgGGGAAACCCGGCCCCCACCCCGGGGTGGCCACACATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGc                                    <  2:3095973/146‑1 (MQ=11)
                                                                                                                     gcgcgGGTAAACCCGCTCCTCCCCCGGGGTTGACACACATGTCGGTGGTAGTTGTCTCGTGTCGTGCGGCGTAGAGGGCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCtt                                  <  1:1945747/149‑1 (MQ=9)
                                                                                                                     acgggggTAAACCCGCCCGGCCCACGGGGTGGCCACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGACGAGGCCGGGCCTTCGGGCCAAACTCCCCCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCtt                                  <  1:2222903/146‑1 (MQ=14)
                                                                                                                       gcgGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGt                                >  2:308668/1‑150 (MQ=255)
                                                                                                                       gcgGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGt                                >  1:1766284/1‑150 (MQ=255)
                                                                                                                       gcgGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGt                                >  2:4260196/1‑150 (MQ=255)
                                                                                                                       gcgGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGt                                >  2:4475062/1‑150 (MQ=255)
                                                                                                                       gcgGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGt                                >  1:1356848/1‑150 (MQ=255)
                                                                                                                         gggTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTAt                              >  2:1222802/1‑150 (MQ=255)
                                                                                                                         ggcaaaccccGCCCCTCCCCCGCCGAGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCCCCCGCGCAGAGGCCGCGCCAGCTGGACAATCTCCACCTGCCGGTCCGCGTTGCCCCCTCCACCGTATTCGCGGGTAAACCCGCTTGTAt                              <  2:3752175/143‑1 (MQ=14)
                                                                                                                          ggTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATg                             <  2:574104/150‑1 (MQ=255)
                                                                                                                          ggTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCGTGTATg                             >  1:357585/1‑150 (MQ=255)
                                                                                                                           gTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATgg                            >  1:3105294/1‑150 (MQ=255)
                                                                                                                            tAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGt                           >  2:1159895/1‑150 (MQ=255)
                                                                                                                             aaaCCCGCTCCTACCCCGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTa                          <  2:3659377/150‑1 (MQ=255)
                                                                                                                             aaaCCCGCTCCTACCACGGCGTGGTCACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTa                          <  2:1464693/150‑1 (MQ=255)
                                                                                                                             aaaCCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTa                          <  2:4039736/150‑1 (MQ=255)
                                                                                                                             aaaCCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTa                          <  2:1514478/150‑1 (MQ=255)
                                                                                                                                cccGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGa                       >  2:613096/1‑150 (MQ=255)
                                                                                                                                cccGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGa                       >  1:3862422/1‑150 (MQ=255)
                                                                                                                                cccGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGa                       >  2:2257950/1‑150 (MQ=255)
                                                                                                                                cccGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGa                       >  2:2257966/1‑150 (MQ=255)
                                                                                                                                 ccGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGAc                      >  2:4007185/1‑150 (MQ=255)
                                                                                                                                    gCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACtt                    <  1:2909556/150‑1 (MQ=255)
                                                                                                                                       ccTACCACGGGGTGGACACAGATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGaa                 <  2:4222443/150‑1 (MQ=255)
                                                                                                                                        cTACCACGGGGTGGACACACATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCTAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAg                <  1:3109974/150‑1 (MQ=255)
                                                                                                                                        cTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAg                >  2:3138381/1‑150 (MQ=255)
                                                                                                                                        cTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAg                >  1:2996176/1‑150 (MQ=255)
                                                                                                                                        cTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAg                >  2:4173715/1‑150 (MQ=255)
                                                                                                                                          aCCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAggg              <  1:544267/150‑1 (MQ=255)
                                                                                                                                           ccACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGgga             <  1:1008961/150‑1 (MQ=255)
                                                                                                                                           ccACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGgga             >  2:106273/1‑150 (MQ=255)
                                                                                                                                           ccACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGgga             <  1:193520/150‑1 (MQ=255)
                                                                                                                                            cACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggag            >  1:2019799/1‑150 (MQ=255)
                                                                                                                                             aCGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggagg           >  2:3593773/1‑150 (MQ=255)
                                                                                                                                             aCGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggagg           >  2:3593339/1‑150 (MQ=255)
                                                                                                                                              cGGCGTGGACACACATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggagga          <  2:3253962/150‑1 (MQ=255)
                                                                                                                                              cGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggagga          >  1:1602620/1‑150 (MQ=255)
                                                                                                                                              cGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggagga          <  2:835040/150‑1 (MQ=255)
                                                                                                                                                   tGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTcgc                                                >  2:4154235/1‑107 (MQ=40)
                                                                                                                                                   tGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTcgc                                                <  1:4154235/107‑1 (MQ=40)
                                                                                                                                                   tGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGaggaggagg     <  1:158260/150‑1 (MQ=255)
                                                                                                                                                   tGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGaggaggagg     <  1:4482428/150‑1 (MQ=255)
                                                                                                                                                   tGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGaggaggagg     >  1:3652795/1‑150 (MQ=255)
                                                                                                                                                     gACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGc   >  2:289722/1‑150 (MQ=255)
                                                                                                                                                      acacacATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGcc  >  1:3191867/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCCTACCACGGCGGTGGGCACATATCGGGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGCGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACCACGGCGTGGACACACATCGCGGTCGTCGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGGCCTGCTGGCCAATCTCCACCTGCCGTTCAGCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCC  >  NC_002947/5758520‑5758818

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