Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,121,597 T→A 13.8% intergenic (+60/+71) 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,5970TA13.8% 346.8 / 40.8 131intergenic (+60/+71)PP_4505/PP_4506sensor histidine kinase/nitroreductase family protein
Reads supporting (aligned to +/- strand):  ref base T (47/65);  new base A (4/14);  total (51/80)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

CATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGA  >  NC_002947/5121448‑5121744
                                                                                                                                                     |                                                                                                                                                   
cATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCGAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGAt                                                                                                                                                     >  1:835382/1‑150 (MQ=255)
cATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGAt                                                                                                                                                     >  2:4879364/1‑150 (MQ=255)
    gCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtg                                                                                                                                                 <  1:5363479/150‑1 (MQ=255)
      ggTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgtt                                                                                                                                               <  1:345111/150‑1 (MQ=255)
        tAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgttgt                                                                                                                                             <  1:3317826/150‑1 (MQ=255)
        tAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCtgttgt                                                                                                                                             <  2:853470/150‑1 (MQ=255)
         aGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTa                                                                                                                                            <  2:1738830/150‑1 (MQ=255)
         aGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTa                                                                                                                                            <  2:2751930/150‑1 (MQ=255)
            tgtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtgccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGa                                                                                                                                         <  1:2154937‑M2/24‑1 (MQ=255)
            tgtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtgccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGa                                                                                                                                         <  1:2154978‑M2/24‑1 (MQ=255)
            tgtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGa                                                                                                                                         <  1:49391‑M2/24‑1 (MQ=255)
            tgtctggcgttggtcgtgagatgggttcagcggccatcgagcattacacggagctgaaatcggtgttgatcaagctctgatcctttttcctgtaccggcctcttcgcgggtgaacccgctcctacaGGGCATGTGATACCTGGTGTAGGa                                                                                                                                         <  1:3988938‑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)
                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                                                                                                                                       <  2:3787805/145‑1 (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: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)
                        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)
                               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)
                                         gtTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGc                                                                                                            <  2:4366002/150‑1 (MQ=255)
                                            tgtgGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCTGTGGATGCTGTTGGAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCtt                                                                                                         >  1:523341/1‑150 (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)
                                                   cggCTGGTGGTGCCCTTGGCTGGCGCGGCAGGCAGTTGAGGGGGTGCAACGGCGCCCGGCCTCTGCGCGGGTAACCCCGCTCCTACCGGGGTTGTGGTTGCTGTTGTAGGAGCGGGTACACCCGCGAAGAGGCCGGATGCCTTGATGCAt                                                                                                  <  1:2543809/149‑1 (MQ=255)
                                                         gAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCa                                                                                            <  2:2390738/150‑1 (MQ=255)
                                                         gAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCa                                                                                            <  1:4946674/150‑1 (MQ=255)
                                                          aGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAg                                                                                           <  2:4856244/150‑1 (MQ=255)
                                                            ttGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCAt                                                                                                  >  1:2145045/1‑141 (MQ=255)
                                                            ttGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCAt                                                                                                  <  2:2145045/141‑1 (MQ=255)
                                                               cccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgc                                                                                                                           >  1:929735/1‑113 (MQ=255)
                                                               cccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCcgc                                                                                                                           <  2:929735/113‑1 (MQ=255)
                                                               cccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCTAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCaa                                                                                      <  1:1118106/150‑1 (MQ=255)
                                                               cccTTGGCGGGAGCGGCTGGCTGTTGATGGGGTGCAATGGCGCCCGGCATCTTCGCGGGTAAACCCGCTCCTAATGGGGTTGTGTTTGCTTTTGTCGGCGCGGGTTAACCCGGGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCaa                                                                                      <  1:1118479/150‑1 (MQ=255)
                                                                ccTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAg                                                                                     >  1:3357952/1‑150 (MQ=255)
                                                                     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)
                                                                                          tGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGTGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTaa                                                           <  1:841583/150‑1 (MQ=255)
                                                                                          tGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTaa                                                           <  1:833619/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)
                                                                                                           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)
                                                                                                                              cccgctcctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCcgcg                       <  2:835382‑M2/138‑1 (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)
                                                                                                                                 gctcctacaGGGGATCTGGATGCTGATGTAGGAGCGGTTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTGTCCAGCACCGGCGCCGTCCgcgggt                    >  2:5261811‑M2/10‑148 (MQ=255)
                                                                                                                                   tcctaccGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCg                  <  2:1920856‑M2/143‑1 (MQ=255)
                                                                                                                                     ctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGt                >  2:4537435‑M2/6‑150 (MQ=255)
                                                                                                                                     ctacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGt                <  2:1319004‑M2/145‑1 (MQ=255)
                                                                                                                                      tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGt                                              >  1:3193987‑M2/5‑119 (MQ=255)
                                                                                                                                      tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGt                                              <  2:3193987‑M2/115‑1 (MQ=255)
                                                                                                                                      tacaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTc               >  1:1426941‑M2/5‑150 (MQ=255)
                                                                                                                                       acaGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCa              >  2:1020698‑M2/4‑150 (MQ=255)
                                                                                                                                         aGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGc            <  2:3054518‑M2/149‑1 (MQ=255)
                                                                                                                                          ggggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcg           >  1:4649263/1‑150 (MQ=255)
                                                                                                                                           gggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcgg          >  1:4909114/1‑150 (MQ=255)
                                                                                                                                            ggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcggc         >  2:5159889/1‑150 (MQ=255)
                                                                                                                                            ggATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAgcggc         >  2:5159797/1‑150 (MQ=255)
                                                                                                                                                gtgGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTAc     <  2:5130601/150‑1 (MQ=255)
                                                                                                                                                   gATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGa  <  1:3598016/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CATGGCGGTAGCCTGATGCTGGAGAATCATCCACAGGGTGGGTTTGTGGCCAGGCTGGAGTTGCCCTTGGCTGGAGCGGCAGGCAGTTGATGGGGTGCAAAGGCGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGATGTGGATGCTGTTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGATGCCTTGATGCATGAATCAGCCCAAGGCCAGGTTACTCACCCCAGCTACTAACGAAATCAGCCGTTTCCAGCACCGGCGCCGTCCGCGGTTCGGTCAGCGGCGTACCGA  >  NC_002947/5121448‑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: