Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,121,602 T→G 14.6% intergenic (+65/+66) PP_4505 → / ← PP_4506 sensor histidine kinase/nitroreductase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,121,6020TG14.6% 298.8 / 43.6 124intergenic (+65/+66)PP_4505/PP_4506sensor histidine kinase/nitroreductase family protein
Reads supporting (aligned to +/- strand):  ref base T (47/58);  new base G (4/14);  total (52/72)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

Only 100 of 132 total aligned reads displayed.

CGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAG  >  NC_002947/5121453‑5121751
                                                                                                                                                     |                                                                                                                                                     
 ggTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgtt                                                                                                                                                      <  1:345111/150‑1 (MQ=255)
   tAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgttgt                                                                                                                                                    <  1:3317826/150‑1 (MQ=255)
   tAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgttgt                                                                                                                                                    <  2:853470/150‑1 (MQ=255)
    aGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTa                                                                                                                                                   <  2:1738830/150‑1 (MQ=255)
       tgtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtgccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGa                                                                                                                                                <  1:2154978‑M2/24‑1 (MQ=255)
       tgtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtgccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGa                                                                                                                                                <  1:2154937‑M2/24‑1 (MQ=255)
       tgtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGa                                                                                                                                                <  1:49391‑M2/24‑1 (MQ=255)
        gtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtaaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAg                                                                                                                                               >  2:1308995‑M2/126‑150 (MQ=255)
         tctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtgccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGc                                                                                                                                              <  1:4209709‑M2/26‑1 (MQ=255)
         tctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGc                                                                                                                                              <  1:2238959‑M2/26‑1 (MQ=255)
         tctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtaaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGc                                                                                                                                              <  1:1308995‑M2/26‑1 (MQ=255)
         gATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGc                                                                                                                                              <  1:4491405/150‑1 (MQ=255)
           tggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtgccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCgg                                                                                                                                            <  1:2012618‑M2/28‑1 (MQ=255)
           tggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCgg                                                                                                                                            <  2:4885672‑M2/28‑1 (MQ=255)
           tggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCgg                                                                                                                                            <  2:2590247‑M2/28‑1 (MQ=255)
           tggcgttggtcgtgagatgggttcagcggccatcgagcagtacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCgg                                                                                                                                            <  2:4885730‑M2/28‑1 (MQ=255)
           tGCTGGAGCATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTCCAGGGGATGTGGATGCTGTTGTAGGAGCgg                                                                                                                                            <  2:4949432/150‑1 (MQ=255)
           tGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCgg                                                                                                                                            >  2:1514908/1‑150 (MQ=255)
           ccgcggacggcgccggtgctggaaacggctgatttcgttagtagctggggtgagtaacctggccttgggctgattcatgcatcaaggcatccGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCgg                                                                                                                                            >  1:3243862‑M1/93‑150 (MQ=255)
            gCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCggg                                                                                                                                           <  1:4606260/150‑1 (MQ=255)
              tGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGc                                                                                                                                              >  1:3787805/1‑145 (MQ=255)
              tGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGc                                                                                                                                              <  2:3787805/145‑1 (MQ=255)
               ggAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTc                                                                                                                                        <  1:2230653/150‑1 (MQ=255)
                gagaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCa                                                                                                                                       >  2:3371884/1‑150 (MQ=255)
                gagaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCa                                                                                                                                       >  2:148566/1‑150 (MQ=255)
                  gaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCAcc                                                                                                                                     <  1:2508027/150‑1 (MQ=255)
                     cgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtaaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCgg                                                                                                                                            >  1:610606‑M2/113‑140 (MQ=255)
                     cgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtaaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCgg                                                                                                                                            <  2:610606‑M2/28‑1 (MQ=255)
                     cgccggtgctggaaacggctgatttcgttagtagctggggtgagtaacctggccttgggctgattcatgcatcaaggcatccGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgc                                                                                                                                  >  1:1709973‑M1/83‑150 (MQ=255)
                      catcCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgcg                                                                                                                                 <  2:4930956/150‑1 (MQ=255)
                          gtgctggaaacggctgatttcgttagtagctggggtgagtaacctggccttgggctgattcatgcatcaaggcatccGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                             >  1:4592209‑M1/78‑150 (MQ=255)
                          caccGGGTGGGTTTGTGGCCAGGCTGGAGTGGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                             <  1:4804084/150‑1 (MQ=255)
                          cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAATCCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                             >  2:3282126/1‑150 (MQ=255)
                          cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCAc                                                                                                                                      <  2:95446/141‑1 (MQ=255)
                          cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                             >  1:862695/1‑150 (MQ=255)
                          cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                             >  1:863442/1‑150 (MQ=255)
                          cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTGAACCcgctccta                                                                                                                             <  1:4879364/150‑6 (MQ=255)
                             gggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCGGGTTTAcc                                                                                                                                     <  1:1176531‑M2/35‑1 (MQ=255)
                                    gtTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGc                                                                                                                   <  2:4366002/150‑1 (MQ=255)
                                      ttGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCt                                                                                                                 <  2:4476929/150‑1 (MQ=255)
                                        gtgGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTg                                                                                                               <  2:1037230/150‑1 (MQ=255)
                                        gtgGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGGTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTg                                                                                                               <  2:1037241/150‑1 (MQ=255)
                                          ggCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGAt                                                                                                             >  2:3796950/1‑150 (MQ=255)
                                          ggCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGAt                                                                                                             <  2:990432/150‑1 (MQ=255)
                                             cAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCa                                                                                                          >  2:1163457/1‑150 (MQ=255)
                                              cggCTGGTGGTGCCCTTGGCTGGCGCGGCAGGCAGTTGAGGGGGTGCAACGGCGCCCGGCCTCTGCGCGGGTAACCCCGCTCCTACCGGGGTTGTGGTTGCTGTTGTAGGAGCGGGTACACCCGCGAAGAGGCCGGATGCCTTGATGCAt                                                                                                         <  1:2543809/149‑1 (MQ=255)
                                                gttGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGa                                                                                                       >  2:1670880/3‑150 (MQ=255)
                                                    gAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCa                                                                                                   <  1:4946674/150‑1 (MQ=255)
                                                    gAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCa                                                                                                   <  2:2390738/150‑1 (MQ=255)
                                                     aGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAg                                                                                                  <  2:4856244/150‑1 (MQ=255)
                                                       ttGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCAt                                                                                                         >  1:2145045/1‑141 (MQ=255)
                                                       ttGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCAt                                                                                                         <  2:2145045/141‑1 (MQ=255)
                                                          cccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgc                                                                                                                                  >  1:929735/1‑113 (MQ=255)
                                                          cccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgc                                                                                                                                  <  2:929735/113‑1 (MQ=255)
                                                          cccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCTAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCaa                                                                                             <  1:1118106/150‑1 (MQ=255)
                                                          cccTTGGCGGGAGCGGCTGGCTGTTGATGGGGTGCAATGGCGCCCGGCATCTTCGCGGGTAAACCCGCTCCTAATGGGGTTGTGTTTGCTTTTGTCGGCGCGGGTTAACCCGGGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCaa                                                                                             <  1:1118479/150‑1 (MQ=255)
                                                           ccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAg                                                                                            >  1:3357952/1‑150 (MQ=255)
                                                           ccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAg                                                                                            >  2:3869375/1‑150 (MQ=255)
                                                                gCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                                       >  1:3771538/1‑150 (MQ=255)
                                                                gCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                                       >  1:3771535/1‑150 (MQ=255)
                                                                  tGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGt                                                                                     <  2:3662338/150‑1 (MQ=255)
                                                                   ggAGCGGCAGGCAGTTGATGGGGTGCAACGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGtt                                                                                    <  1:2668359/150‑1 (MQ=255)
                                                                     aGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTAc                                                                                  <  2:862695/150‑1 (MQ=255)
                                                                            ggcagTGGATGGGGGGCAAAGGCGCCAGGCCTCTTCGCGGGTAACCCCGCTCCTACCGGGGATGTGGCGGGTGTTGTCGGCGCGGGTTCACCCGCGACGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCAcccc                                                                           <  1:1389216/150‑1 (MQ=255)
                                                                              aagctctgatcctttttcctgtaccggcctcttcgcgggtaaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGc                                                                                                                                              >  2:4276383‑M2/56‑81 (MQ=255)
                                                                              aagctctgatcctttttcctgtaccggcctcttcgcgggtaaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGc                                                                                                                                              <  1:4276383‑M2/26‑1 (MQ=255)
                                                                                     tGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGTGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTaa                                                                  <  1:841583/150‑1 (MQ=255)
                                                                                     tGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTaa                                                                  <  1:833619/150‑1 (MQ=255)
                                                                                        ggTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGAATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGa                                                               <  1:3796950/150‑1 (MQ=255)
                                                                                             aaaGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCa                                                          <  2:2724272/150‑1 (MQ=255)
                                                                                             aaaGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCa                                                          <  1:5318329/150‑1 (MQ=255)
                                                                                             aaaGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCa                                                          <  2:2724202/150‑1 (MQ=255)
                                                                                             aaaGGCGCCAGGCCTCTTCGCGGGTAAAACCCCTCCGACAGGGGATGTGGATGCTGTTGTTGGAGCGGGTTCACCCGCGAAGAGGCCGGAAGCCATGATGCAAGAATCAGCCCAACGCCAGGTTACAAAACCCACCAACTAAAAAAAtaa                                                          >  1:4184292/1‑148 (MQ=25)
                                                                                                 gcgcCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCg                                                      >  2:3308496/1‑150 (MQ=255)
                                                                                                      aGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTcc                                                 <  1:2194678/150‑1 (MQ=255)
                                                                                                                    gtaaacccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGt                                   >  1:678264‑M2/18‑150 (MQ=255)
                                                                                                                     taaacccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTc                                  >  1:696522‑M2/17‑150 (MQ=255)
                                                                                                                         cccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCcgcg                              <  2:835382‑M2/138‑1 (MQ=255)
                                                                                                                         cccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCAACCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCcgcg                              <  2:5294896‑M2/138‑1 (MQ=255)
                                                                                                                          ccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCgg                             >  1:4860165‑M2/12‑150 (MQ=255)
                                                                                                                          ccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCgg                             >  2:2751746‑M2/12‑150 (MQ=255)
                                                                                                                          ccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCgg                             >  2:2002519‑M2/12‑150 (MQ=255)
                                                                                                                            gctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                                       <  1:1369418‑M2/81‑1 (MQ=255)
                                                                                                                            gctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                                       >  2:1369418‑M2/10‑90 (MQ=255)
                                                                                                                            gctcctacaGGGGATCTGGATGCTGATGTAGGAGCGGTTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTGTCCAGCACCGGCGCCGTCCgcgggt                           >  2:5261811‑M2/10‑148 (MQ=255)
                                                                                                                              tcctaccGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCg                         <  2:1920856‑M2/143‑1 (MQ=255)
                                                                                                                                ctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGt                       >  1:5048189‑M2/6‑150 (MQ=255)
                                                                                                                                ctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGt                       <  2:1319004‑M2/145‑1 (MQ=255)
                                                                                                                                 tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGt                                                     >  1:3193987‑M2/5‑119 (MQ=255)
                                                                                                                                 tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGt                                                     <  2:3193987‑M2/115‑1 (MQ=255)
                                                                                                                                 tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTc                      >  1:1426941‑M2/5‑150 (MQ=255)
                                                                                                                                  acaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCa                     >  2:1020698‑M2/4‑150 (MQ=255)
                                                                                                                                    aGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGc                   <  2:3054518‑M2/149‑1 (MQ=255)
                                                                                                                                       ggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcggc                >  2:5159797/1‑150 (MQ=255)
                                                                                                                                       ggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcggc                >  2:5159889/1‑150 (MQ=255)
                                                                                                                                           gtgGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTAc            <  2:5130601/150‑1 (MQ=255)
                                                                                                                                              gATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGa         <  1:3598016/150‑1 (MQ=255)
                                                                                                                                                 gCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACAt      >  2:3265201/1‑150 (MQ=255)
                                                                                                                                                   tgttgtAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATAc    <  2:5048189/150‑1 (MQ=255)
                                                                                                                                                     ttgtAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAg  >  1:201315/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACATACAG  >  NC_002947/5121453‑5121751

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