Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,121,589 G→C 13.6% intergenic (+52/+79) 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,5890GC13.6% 362.6 / 45.0 132intergenic (+52/+79)PP_4505/PP_4506sensor histidine kinase/nitroreductase family protein
Reads supporting (aligned to +/- strand):  ref base G (50/64);  new base C (4/14);  total (54/78)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

Only 100 of 137 total aligned reads displayed.

CGGAGCGCCATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGC  >  NC_002947/5121440‑5121737
                                                                                                                                                     |                                                                                                                                                    
cGGAGCGCCATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAgggg                                                                                                                                                      >  1:5294896/1‑150 (MQ=255)
  gAGCGCCATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGAt                                                                                                                                                    >  2:2868211/1‑150 (MQ=255)
       ccATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGa                                                                                                                                               >  1:5130601/1‑150 (MQ=255)
       ccATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGa                                                                                                                                               >  1:5130580/1‑150 (MQ=255)
       ccATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGa                                                                                                                                               <  2:1906468/150‑1 (MQ=255)
       ccATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGa                                                                                                                                               >  2:4438618/1‑150 (MQ=255)
        cATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCGAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGAt                                                                                                                                              >  1:835382/1‑150 (MQ=255)
        cATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGAt                                                                                                                                              >  2:4879364/1‑150 (MQ=255)
            gCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtg                                                                                                                                          <  1:1434143/150‑1 (MQ=255)
            gCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtg                                                                                                                                          <  1:5363479/150‑1 (MQ=255)
             cGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgt                                                                                                                                         <  2:5182467/150‑1 (MQ=255)
              ggTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgtt                                                                                                                                        <  1:345111/150‑1 (MQ=255)
                tAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgttgt                                                                                                                                      <  1:3317826/150‑1 (MQ=255)
                tAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgttgt                                                                                                                                      <  1:2193008/150‑1 (MQ=255)
                tAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgttgt                                                                                                                                      <  2:853470/150‑1 (MQ=255)
                 aGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTa                                                                                                                                     <  2:2751930/150‑1 (MQ=255)
                    tgtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtgccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGa                                                                                                                                  <  1:2154978‑M2/24‑1 (MQ=255)
                    tgtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGa                                                                                                                                  <  1:49391‑M2/24‑1 (MQ=255)
                      tctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtgccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGc                                                                                                                                <  1:4209709‑M2/26‑1 (MQ=255)
                      tctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGc                                                                                                                                <  1:2238959‑M2/26‑1 (MQ=255)
                      gATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGc                                                                                                                                <  1:4491405/150‑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)
                           tGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGc                                                                                                                                >  1:3787805/1‑145 (MQ=255)
                           cggacggcgccggtgctggaaacggctgatttcgttagtagctggggtgagtaacctggccttgggctgattcatgcatcaaggcatccGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGtt                                                                                                                           <  2:1709973‑M1/61‑1 (MQ=255)
                            ggAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTc                                                                                                                          <  1:2230653/150‑1 (MQ=255)
                             gagaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCa                                                                                                                         >  2:26094/1‑150 (MQ=255)
                             gagaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCa                                                                                                                         >  2:3371884/1‑150 (MQ=255)
                                aaTCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCAccc                                                                                                                      <  1:2076510/150‑1 (MQ=255)
                                aaTCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGAGTTCAccc                                                                                                                      <  1:4353531/150‑1 (MQ=255)
                                  cgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtaaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCgg                                                                                                                              <  2:610606‑M2/28‑1 (MQ=255)
                                  cgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtaaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCgg                                                                                                                              >  1:610606‑M2/113‑140 (MQ=255)
                                  cgccggtgctggaaacggctgatttcgttagtagctggggtgagtaacctggccttgggctgattcatgcatcaaggcatccGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgc                                                                                                                    >  1:1709973‑M1/83‑150 (MQ=255)
                                   catcCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgcg                                                                                                                   <  2:2078978/150‑1 (MQ=255)
                                   catcCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgcg                                                                                                                   <  2:4930956/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                                                                                                                        <  2:95446/141‑1 (MQ=255)
                                       cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCAc                                                                                                                        >  1:95446/1‑141 (MQ=255)
                                       cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                               >  1:863442/1‑150 (MQ=255)
                                       cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                               >  1:862695/1‑150 (MQ=255)
                                          gggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCGGGTTTAcc                                                                                                                       >  2:1176531‑M2/105‑139 (MQ=255)
                                          gggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCGGGTTTAcc                                                                                                                       <  1:1176531‑M2/35‑1 (MQ=255)
                                                 gtTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGc                                                                                                     <  2:4366002/150‑1 (MQ=255)
                                                     gtgGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTg                                                                                                 <  2:1037230/150‑1 (MQ=255)
                                                       ggCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGAt                                                                                               >  2:3796950/1‑150 (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                                                                                     <  2:2390738/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)
                                                                       cccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCTAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCaa                                                                               <  1:1118106/150‑1 (MQ=255)
                                                                       cccTTGGCGGGAGCGGCTGGCTGTTGATGGGGTGCAATGGCGCCCGGCATCTTCGCGGGTAAACCCGCTCCTAATGGGGTTGTGTTTGCTTTTGTCGGCGCGGGTTAACCCGGGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCaa                                                                               <  1:1118479/150‑1 (MQ=255)
                                                                        ccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAg                                                                              >  2:3869375/1‑150 (MQ=255)
                                                                        ccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAg                                                                              >  1:3357952/1‑150 (MQ=255)
                                                                             gCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                         >  1:3771538/1‑150 (MQ=255)
                                                                               tGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGt                                                                       >  1:1319004/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)
                                                                                  aGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTAc                                                                    <  2:863442/150‑1 (MQ=255)
                                                                                         ggcagTGGATGGGGGGCAAAGGCGCCAGGCCTCTTCGCGGGTAACCCCGCTCCTACCGGGGATGTGGCGGGTGTTGTCGGCGCGGGTTCACCCGCGACGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCAcccc                                                             <  1:1389216/150‑1 (MQ=255)
                                                                                                  tGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTaa                                                    <  1:833619/150‑1 (MQ=255)
                                                                                                  tGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTaa                                                    <  1:841637/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)
                                                                                                                  cAGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGGGTTGTGGATGCTGTTGTAGTCGGGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTc                                    <  1:2433493/150‑1 (MQ=39)
                                                                                                                   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)
                                                                                                                                      cccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCAACCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCcgcg                <  2:5294896‑M2/138‑1 (MQ=255)
                                                                                                                                       ccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCgg               >  1:4860165‑M2/12‑150 (MQ=255)
                                                                                                                                         gctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                         <  1:1369418‑M2/81‑1 (MQ=255)
                                                                                                                                         gctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                         >  2:1369418‑M2/10‑90 (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)
                                                                                                                                             ctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGt         >  2:4537435‑M2/6‑150 (MQ=255)
                                                                                                                                              tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGt                                       >  1:3193987‑M2/5‑119 (MQ=255)
                                                                                                                                              tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGt                                       <  2:3193987‑M2/115‑1 (MQ=255)
                                                                                                                                              tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTc        >  1:360384‑M2/5‑150 (MQ=255)
                                                                                                                                              tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTc        >  2:785446‑M2/5‑150 (MQ=255)
                                                                                                                                              tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTc        >  2:3598016‑M2/5‑150 (MQ=255)
                                                                                                                                               acaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCa       >  2:1020698‑M2/4‑150 (MQ=255)
                                                                                                                                                  ggggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcg    >  1:4649263/1‑150 (MQ=255)
                                                                                                                                                   gggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcgg   >  1:4909114/1‑150 (MQ=255)
                                                                                                                                                    ggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcggc  >  2:5159889/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGAGCGCCATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGC  >  NC_002947/5121440‑5121737

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: