Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,121,596 A→T 14.8% intergenic (+59/+72) 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,5960AT14.8% 318.7 / 43.2 133intergenic (+59/+72)PP_4505/PP_4506sensor histidine kinase/nitroreductase family protein
Reads supporting (aligned to +/- strand):  ref base A (49/63);  new base T (4/16);  total (53/80)
Fisher's exact test for biased strand distribution p-value = 5.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

Only 100 of 137 total aligned reads displayed.

CCATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGA  >  NC_002947/5121447‑5121744
                                                                                                                                                     |                                                                                                                                                    
ccATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGa                                                                                                                                                      >  1:5130601/1‑150 (MQ=255)
ccATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGa                                                                                                                                                      >  1:5130580/1‑150 (MQ=255)
 cATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCGAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGAt                                                                                                                                                     >  1:835382/1‑150 (MQ=255)
 cATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGAt                                                                                                                                                     >  2:4879364/1‑150 (MQ=255)
  aTGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATg                                                                                                                                                    <  1:5434590/150‑1 (MQ=255)
     gCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtg                                                                                                                                                 <  1:5363479/150‑1 (MQ=255)
     gCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtg                                                                                                                                                 <  1:1434143/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:1738830/150‑1 (MQ=255)
             tgtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtgccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGa                                                                                                                                         <  1:2154978‑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)
               gATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGc                                                                                                                                       <  1:4491405/150‑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)
                 ccgcggacggcgccggtgctggaaacggctgatttcgttagtagctggggtgagtaacctggccttgggctgattcatgcatcaaggcatccGGCCTCTTCGCGGGTGAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCgg                                                                                                                                     >  1:3243862‑M1/93‑150 (MQ=255)
                  gCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCggg                                                                                                                                    <  1:332855/150‑1 (MQ=255)
                    tGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGc                                                                                                                                       <  2:3787805/145‑1 (MQ=255)
                    tGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGc                                                                                                                                       >  1:3787805/1‑145 (MQ=255)
                    cggacggcgccggtgctggaaacggctgatttcgttagtagctggggtgagtaacctggccttgggctgattcatgcatcaaggcatccGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGtt                                                                                                                                  <  2:1709973‑M1/61‑1 (MQ=255)
                      gagaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCa                                                                                                                                >  2:3371884/1‑150 (MQ=255)
                      gagaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCa                                                                                                                                >  2:148566/1‑150 (MQ=255)
                        gaATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCAcc                                                                                                                              <  1:2508027/150‑1 (MQ=255)
                         aaTCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCAccc                                                                                                                             <  1:2076510/150‑1 (MQ=255)
                         aaTCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGAGTTCAccc                                                                                                                             <  1:4353531/150‑1 (MQ=255)
                           cgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtaaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCgg                                                                                                                                     <  2:610606‑M2/28‑1 (MQ=255)
                           cgccggtgctggaaacggctgatttcgttagtagctggggtgagtaacctggccttgggctgattcatgcatcaaggcatccGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgc                                                                                                                           >  1:1709973‑M1/83‑150 (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:863442/1‑150 (MQ=255)
                                cacaGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAaga                                                                                                                      >  1:862695/1‑150 (MQ=255)
                                   gggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCGGGTTTAcc                                                                                                                              <  1:1176531‑M2/35‑1 (MQ=255)
                                   gggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGAGCGGGTTTAcc                                                                                                                              >  2:1176531‑M2/105‑139 (MQ=255)
                                       tgggtTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGa                                                                                                               >  1:3054518/1‑150 (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)
                                                  cctggCGGGAGGTGCCCGGGGCTGGAGCGGCCGGCAGTTGCGGGGGTGCAACGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGGGTTGTGGATGGTGTTGTAGGCGGGGGTCCACCCGCGAAGAGGCCGGATGCCTTGATGc                                                                                                    <  1:148566/147‑1 (MQ=255)
                                                      gttGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGa                                                                                                >  2:1670880/3‑150 (MQ=255)
                                                          gAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCa                                                                                            <  2:2390738/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)
                                                                cccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCTAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCaa                                                                                      <  1:1118106/150‑1 (MQ=255)
                                                                cccTTGGCGGGAGCGGCTGGCTGTTGATGGGGTGCAATGGCGCCCGGCATCTTCGCGGGTAAACCCGCTCCTAATGGGGTTGTGTTTGCTTTTGTCGGCGCGGGTTAACCCGGGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCaa                                                                                      <  1:1118479/150‑1 (MQ=255)
                                                                      gCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                                >  1:3771538/1‑150 (MQ=255)
                                                                        tGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGt                                                                              <  2:3662338/150‑1 (MQ=255)
                                                                        tGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGt                                                                              >  1:1319004/1‑150 (MQ=255)
                                                                         ggAGCGGCAGGCAGTTGATGGGGTGCAACGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGtt                                                                             <  1:2668359/150‑1 (MQ=255)
                                                                           aGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTAc                                                                           <  2:863442/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)
                                                                                           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                                                   <  2:2724202/150‑1 (MQ=255)
                                                                                                   aaaGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCa                                                   <  1:5318329/150‑1 (MQ=255)
                                                                                                   aaaGGCGCCAGGCCTCTTCGCGGGTAAAACCCCTCCGACAGGGGATGTGGATGCTGTTGTTGGAGCGGGTTCACCCGCGAAGAGGCCGGAAGCCATGATGCAAGAATCAGCCCAACGCCAGGTTACAAAACCCACCAACTAAAAAAAtaa                                                   >  1:4184292/1‑148 (MQ=25)
                                                                                                       gcgcCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCg                                               >  2:3308496/1‑150 (MQ=255)
                                                                                                           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:696518‑M2/17‑150 (MQ=255)
                                                                                                                               cccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCAACCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCcgcg                       <  2:5294896‑M2/138‑1 (MQ=255)
                                                                                                                                ccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCgg                      >  2:2751746‑M2/12‑150 (MQ=255)
                                                                                                                                ccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCgg                      >  1:4860165‑M2/12‑150 (MQ=255)
                                                                                                                                ccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCgg                      >  2:2002519‑M2/12‑150 (MQ=255)
                                                                                                                                  gctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                                >  2:1369418‑M2/10‑90 (MQ=255)
                                                                                                                                  gctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAg                                                                                <  1:1369418‑M2/81‑1 (MQ=255)
                                                                                                                                    tcctaccGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCg                  <  2:1920856‑M2/143‑1 (MQ=255)
                                                                                                                                      ctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGt                >  1:5048189‑M2/6‑150 (MQ=255)
                                                                                                                                      ctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGt                >  2:4537435‑M2/6‑150 (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)
                                                                                                                                        acaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCa              >  2:1020698‑M2/4‑150 (MQ=255)
                                                                                                                                          aGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGc            <  2:3054518‑M2/149‑1 (MQ=255)
                                                                                                                                           ggggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcg           >  1:4649263/1‑150 (MQ=255)
                                                                                                                                            gggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcgg          >  1:4909114/1‑150 (MQ=255)
                                                                                                                                             ggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcggc         >  2:5159889/1‑150 (MQ=255)
                                                                                                                                                 gtgGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTAc     <  2:5130601/150‑1 (MQ=255)
                                                                                                                                                 gtgGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTAc     <  2:5130580/150‑1 (MQ=255)
                                                                                                                                                    gATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGa  <  1:3598016/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGA  >  NC_002947/5121447‑5121744

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: