Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,121,598 G→C 13.7% intergenic (+61/+70) 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,5980GC13.7% 362.3 / 42.7 130intergenic (+61/+70)PP_4505/PP_4506sensor histidine kinase/nitroreductase family protein
Reads supporting (aligned to +/- strand):  ref base G (47/65);  new base C (4/14);  total (51/79)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 132 total aligned reads displayed.

ATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACAT  >  NC_002947/5121449‑5121747
                                                                                                                                                     |                                                                                                                                                     
aTGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATg                                                                                                                                                       <  1:5434590/150‑1 (MQ=255)
   gCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtg                                                                                                                                                    <  1:5363479/150‑1 (MQ=255)
     ggTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgtt                                                                                                                                                  <  1:345111/150‑1 (MQ=255)
       tAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgttgt                                                                                                                                                <  2:853470/150‑1 (MQ=255)
       tAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgttgt                                                                                                                                                <  1:3317826/150‑1 (MQ=255)
        aGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTa                                                                                                                                               <  2:2751930/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:2590247‑M2/28‑1 (MQ=255)
               tggcgttggtcgtgagatgggttcagcggccatcgagcagtacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCgg                                                                                                                                        <  2:4885730‑M2/28‑1 (MQ=255)
               tGCTGGAGCATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTCCAGGGGATGTGGATGCTGTTGTAGGAGCgg                                                                                                                                        <  2:4949432/150‑1 (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)
                  cggacggcgccggtgctggaaacggctgatttcgttagtagctggggtgagtaacctggccttgggctgattcatgcatcaaggcatccGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGtt                                                                                                                                     <  2:1709973‑M1/61‑1 (MQ=255)
                   ggAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTc                                                                                                                                    <  1:2230653/150‑1 (MQ=255)
                    gagaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCa                                                                                                                                   >  2:148566/1‑150 (MQ=255)
                    gagaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCa                                                                                                                                   >  2:26094/1‑150 (MQ=255)
                    gagaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCa                                                                                                                                   >  2:3371884/1‑150 (MQ=255)
                      gaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCAcc                                                                                                                                 <  1:2508027/150‑1 (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)
                          catcCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgcg                                                                                                                             <  2:2078978/150‑1 (MQ=255)
                              gtgctggaaacggctgatttcgttagtagctggggtgagtaacctggccttgggctgattcatgcatcaaggcatccGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                         >  1:4592209‑M1/78‑150 (MQ=255)
                              caccGGGTGGGTTTGTGGCCAGGCTGGAGTGGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                         <  1:4804084/150‑1 (MQ=255)
                              cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAATGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAagt                                                                                                                         >  2:3282114/1‑149 (MQ=255)
                              cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAATCCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                         >  2:3282126/1‑150 (MQ=255)
                              cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCAc                                                                                                                                  >  1:95446/1‑141 (MQ=255)
                              cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCAc                                                                                                                                  <  2:95446/141‑1 (MQ=255)
                              cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                         >  1:862695/1‑150 (MQ=255)
                              cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTGAACCcgctccta                                                                                                                         <  1:4879364/150‑6 (MQ=255)
                                 gggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCGGGTTTAcc                                                                                                                                 >  2:1176531‑M2/105‑139 (MQ=255)
                                 gggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCGGGTTTAcc                                                                                                                                 <  1:1176531‑M2/35‑1 (MQ=255)
                                          ttGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCt                                                                                                             <  2:4476929/150‑1 (MQ=255)
                                            gtgGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTg                                                                                                           <  2:1037230/150‑1 (MQ=255)
                                              ggCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGAt                                                                                                         <  2:990432/150‑1 (MQ=255)
                                                cctggCGGGAGGTGCCCGGGGCTGGAGCGGCCGGCAGTTGCGGGGGTGCAACGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGGGTTGTGGATGGTGTTGTAGGCGGGGGTCCACCCGCGAAGAGGCCGGATGCCTTGATGc                                                                                                       <  1:148566/147‑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)
                                                         aGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAg                                                                                              <  2:4856244/150‑1 (MQ=255)
                                                           ttGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCAt                                                                                                     <  2:2145045/141‑1 (MQ=255)
                                                           ttGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCAt                                                                                                     >  1:2145045/1‑141 (MQ=255)
                                                              cccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgc                                                                                                                              <  2:929735/113‑1 (MQ=255)
                                                              cccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgc                                                                                                                              >  1:929735/1‑113 (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:841637/150‑1 (MQ=255)
                                                                                            ggTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGAATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGa                                                           <  1:3796950/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)
                                                                                                          aGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTcc                                             <  1:2194678/150‑1 (MQ=255)
                                                                                                           ggcctcttcgcgggtaaacccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCa                                            <  2:1447162‑M2/120‑1 (MQ=255)
                                                                                                                      gggtaaacccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCCCCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAGATCAGCCGTTTCCAGCACCggcgct                                 <  2:3942899‑M2/131‑2 (MQ=255)
                                                                                                                        gtaaacccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGt                               >  1:678264‑M2/18‑150 (MQ=255)
                                                                                                                         taaacccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTc                              >  1:696522‑M2/17‑150 (MQ=255)
                                                                                                                         taaacccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTc                              >  1:696518‑M2/17‑150 (MQ=255)
                                                                                                                             cccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCcgcg                          <  2:835382‑M2/138‑1 (MQ=255)
                                                                                                                             cccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCAACCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCcgcg                          <  2:5294896‑M2/138‑1 (MQ=255)
                                                                                                                              ccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCgg                         >  2:2002519‑M2/12‑150 (MQ=255)
                                                                                                                              ccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCgg                         >  2:2751746‑M2/12‑150 (MQ=255)
                                                                                                                              ccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCgg                         >  1:4860165‑M2/12‑150 (MQ=255)
                                                                                                                                gctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                                   >  2:1369418‑M2/10‑90 (MQ=255)
                                                                                                                                gctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                                   <  1:1369418‑M2/81‑1 (MQ=255)
                                                                                                                                gctcctacaGGGGATCTGGATGCTGATGTAGGAGCGGTTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTGTCCAGCACCGGCGCCGTCCgcgggt                       >  2:5261811‑M2/10‑148 (MQ=255)
                                                                                                                                    ctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGt                   >  2:4537435‑M2/6‑150 (MQ=255)
                                                                                                                                    ctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGt                   >  1:5048189‑M2/6‑150 (MQ=255)
                                                                                                                                     tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGt                                                 >  1:3193987‑M2/5‑119 (MQ=255)
                                                                                                                                     tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTc                  >  2:3598016‑M2/5‑150 (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)
                                                                                                                                                     |                                                                                                                                                     
ATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGACAT  >  NC_002947/5121449‑5121747

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: