Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,306,728 A→G 8.0% intergenic (+16/+16) cysS → / ← PP_2906 cysteine‑‑tRNA ligase/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,306,7280AG 8.0% 490.4 / 21.8 162intergenic (+16/+16)cysS/PP_2906cysteine‑‑tRNA ligase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (76/73);  new base G (5/8);  total (81/81)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 165 total aligned reads displayed.

TCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCGG  >  NC_002947/3306582‑3306877
                                                                                                                                                  |                                                                                                                                                     
  cAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCa                                                                                                                                                  >  2:2099506/1‑150 (MQ=255)
  cAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCa                                                                                                                                                  <  2:363325/150‑1 (MQ=255)
     gcgcgcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGcc                                                                                                                                               >  1:2818986/1‑150 (MQ=255)
     acgcgcCTTCAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGcc                                                                                                                                               >  1:915945/2‑150 (MQ=255)
      cgcgcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGccc                                                                                                                                              <  2:2936896/150‑1 (MQ=255)
         gcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACTGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAg                                                                                                                                           <  1:3996025/150‑1 (MQ=255)
         gcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCAg                                                                                                                                           <  1:4299888/150‑1 (MQ=255)
         gcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCAg                                                                                                                                           <  2:1777564/150‑1 (MQ=255)
         gcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCAg                                                                                                                                           <  2:1692759/150‑1 (MQ=255)
          ccTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCAgt                                                                                                                                          <  1:4983631/150‑2 (MQ=255)
          ccTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGc                                                                                                                                          <  1:1651416/150‑1 (MQ=255)
          ccTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGc                                                                                                                                          >  1:2817472/1‑150 (MQ=255)
            tGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCtt                                                                                                                                        >  2:2554145/1‑150 (MQ=255)
                  cgcgTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTAttt                                                                                                                                  <  2:1137272/150‑1 (MQ=255)
                   gcgTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTg                                                                                                                                 >  2:205143/1‑150 (MQ=255)
                   gcgTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTg                                                                                                                                 <  1:2116954/150‑1 (MQ=255)
                   gcgTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTg                                                                                                                                 >  2:205121/1‑150 (MQ=255)
                     gTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCt                                                                                                                               <  2:3253764/150‑1 (MQ=255)
                       gCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTa                                                                                                                             <  2:4196800/150‑1 (MQ=255)
                         ggACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGa                                                                                                                           >  1:1955803/1‑150 (MQ=255)
                         ggACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGa                                                                                                                           >  2:2524258/1‑150 (MQ=255)
                          gACTAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGAc                                                                                                                          <  2:4518352/150‑1 (MQ=255)
                          gACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGAc                                                                                                                          <  2:1240881/150‑1 (MQ=255)
                          gACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGAc                                                                                                                          <  1:1141805/150‑1 (MQ=255)
                           aCAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGAcc                                                                                                                         >  1:3882078/1‑150 (MQ=255)
                            cAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGAccc                                                                                                                        <  1:302499/150‑1 (MQ=255)
                               gAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCa                                                                                                                     >  2:1656654/1‑150 (MQ=255)
                               gAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCa                                                                                                                     >  1:4694506/1‑150 (MQ=255)
                               gAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCa                                                                                                                     >  2:3995574/1‑150 (MQ=255)
                                  cTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGc                                                                                                                  <  2:4694506/150‑1 (MQ=255)
                                    gggCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCgg                                                                                                                <  2:3818527/150‑1 (MQ=255)
                                    gggCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCgg                                                                                                                <  1:81548/150‑1 (MQ=255)
                                       cTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCa                                                                                                             <  2:4110026/150‑1 (MQ=255)
                                       cTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCa                                                                                                             <  2:2546784/150‑1 (MQ=255)
                                       cTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCa                                                                                                             <  1:1288058/150‑1 (MQ=255)
                                        tGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAg                                                                                                            >  1:3571198/1‑150 (MQ=255)
                                        tGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAg                                                                                                            <  2:1705635/150‑1 (MQ=255)
                                            tCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCa                                                                                                        <  2:1556175/150‑1 (MQ=255)
                                             cTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCaa                                                                                                       <  2:5083129/150‑1 (MQ=255)
                                              tGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAAt                                                                                                      <  2:3040484/150‑1 (MQ=255)
                                              tGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAAt                                                                                                      <  1:245221/150‑1 (MQ=255)
                                               gACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTAttt                                                                                                                                  <  1:4449342/121‑1 (MQ=255)
                                                 ccGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATAcc                                                                                                   <  2:4348364/150‑1 (MQ=255)
                                                  cGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCa                                                                                                  <  1:3821259/150‑1 (MQ=255)
                                                  cGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCa                                                                                                  >  1:3009577/1‑150 (MQ=255)
                                                  cGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCa                                                                                                  <  1:3821116/150‑1 (MQ=255)
                                                         gcgATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGc                                                                                           >  2:4377401/1‑150 (MQ=255)
                                                           gATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGccc                                                                                         >  2:4831153/1‑150 (MQ=255)
                                                             tCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCtttt                                                                                                                                      >  1:4689949/1‑103 (MQ=255)
                                                             tCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCtttt                                                                                                                                      <  2:4689949/103‑1 (MQ=255)
                                                               aGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCcccacc                                                                                     <  1:1040858/150‑1 (MQ=255)
                                                                gCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCcccaccc                                                                                    >  1:3746319/1‑150 (MQ=255)
                                                                 cTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCcccaccca                                                                                   >  1:1466199/1‑150 (MQ=255)
                                                                   gACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAAc                                                                                 <  1:4902976/150‑1 (MQ=255)
                                                                     ccGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTc                                                                               <  1:4610099/150‑1 (MQ=255)
                                                                       gCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGc                                                                             >  1:2319527/1‑150 (MQ=255)
                                                                       gCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGc                                                                             >  1:1014594/1‑150 (MQ=255)
                                                                            gggTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTccc                                                                        <  1:1178262/150‑1 (MQ=255)
                                                                             ggTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTcccc                                                                       <  1:3235598/150‑1 (MQ=255)
                                                                                    tGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCACTTCCCCCAAtctc                                                                <  1:3871357/150‑1 (MQ=255)
                                                                                          aggaCAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGccc                                                          <  1:1847820/150‑1 (MQ=255)
                                                                                               aGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAAc                                                     >  1:1785154/1‑150 (MQ=255)
                                                                                                gCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAAcc                                                    >  2:3559973/1‑150 (MQ=255)
                                                                                                 cAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCa                                                   >  1:2212071/1‑150 (MQ=255)
                                                                                                  aaGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCAt                                                  >  1:587003/1‑150 (MQ=255)
                                                                                                   aGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATg                                                 >  1:4885107/1‑150 (MQ=255)
                                                                                                      gCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAg                                              >  1:5011701/1‑150 (MQ=255)
                                                                                                          gACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCCGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTcg                                          >  2:3800304/1‑150 (MQ=255)
                                                                                                                cTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATc                                    >  1:4412724/1‑150 (MQ=255)
                                                                                                                cTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATc                                    <  1:2650260/150‑1 (MQ=255)
                                                                                                                 tGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATcg                                   >  1:4894227/1‑150 (MQ=255)
                                                                                                                      tCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCgccgcc                              >  2:4067990/1‑150 (MQ=255)
                                                                                                                          gCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGa                          <  2:915945/150‑1 (MQ=255)
                                                                                                                          gCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGa                          >  2:370572/1‑150 (MQ=255)
                                                                                                                          gCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGa                          >  1:3473794/1‑150 (MQ=255)
                                                                                                                            tgactgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                        <  1:4128544/150‑1 (MQ=255)
                                                                                                                             gactgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCg                       >  1:4083139/1‑150 (MQ=255)
                                                                                                                             gactgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCg                       >  2:3779944/1‑150 (MQ=255)
                                                                                                                                tgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACa                    >  1:3692690/1‑150 (MQ=255)
                                                                                                                                tgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACa                    >  2:1480088/1‑150 (MQ=255)
                                                                                                                                 gaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACaa                   <  2:2944251/150‑1 (MQ=255)
                                                                                                                                 gaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACaa                   >  1:959974/1‑150 (MQ=255)
                                                                                                                                  aaGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAAc                  >  1:610340/1‑150 (MQ=255)
                                                                                                                                     tAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAAcccc               >  2:1334331/1‑150 (MQ=255)
                                                                                                                                      aaCTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCa              <  2:826648/150‑1 (MQ=255)
                                                                                                                                      aaCTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCa              >  1:2326715/1‑150 (MQ=255)
                                                                                                                                       gcTGGGGCTGCGGTGCAGCCCCAGTTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCa                             <  1:721104/133‑1 (MQ=255)
                                                                                                                                       gcTGGGGCTGCGGTGCAGCCCCAGTTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCa                             <  1:721205/133‑1 (MQ=255)
                                                                                                                                       gcTGGGGCTGCGGTGCAGCCCCAGTTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCa                             >  2:721104/2‑134 (MQ=255)
                                                                                                                                        cTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTTGTCGATCGTGGCGGCGATCGTCGAcctgcat            <  1:2099506/150‑6 (MQ=255)
                                                                                                                                        cTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGc            <  2:755166/150‑1 (MQ=255)
                                                                                                                                         tGGGGCTGCGGTGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGcc           >  2:2400796/1‑150 (MQ=255)
                                                                                                                                         tGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGcc           >  1:1514232/1‑150 (MQ=255)
                                                                                                                                         tGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGcc           <  2:2818986/150‑1 (MQ=255)
                                                                                                                                           gggCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCg         >  1:4681317/1‑150 (MQ=255)
                                                                                                                                             gCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGaa       >  1:2560805/1‑150 (MQ=255)
                                                                                                                                             gCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGaa       <  2:5122383/150‑1 (MQ=255)
                                                                                                                                               tGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGc     <  1:2809015/150‑1 (MQ=255)
                                                                                                                                                  aCCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCgg  >  1:409306/1‑150 (MQ=255)
                                                                                                                                                  aCCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCgg  >  2:2439980/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCGG  >  NC_002947/3306582‑3306877

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: