Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,818,821 C→T 10.5% intergenic (+29/+100) PP_5092 → / ← pilT NLP/P60 protein/twitching motility protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,818,8210CT10.5% 529.0 / 40.3 191intergenic (+29/+100)PP_5092/pilTNLP/P60 protein/twitching motility protein
Reads supporting (aligned to +/- strand):  ref base C (69/102);  new base T (12/8);  total (81/110)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 191 total aligned reads displayed.

CCACGCCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCC  >  NC_002947/5818673‑5818968
                                                                                                                                                    |                                                                                                                                                   
ccacgcCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCt                                                                                                                                                    <  1:68713/150‑1 (MQ=255)
 cacgcCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCtt                                                                                                                                                   >  1:1063356/1‑150 (MQ=255)
 cacgcCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCtt                                                                                                                                                   >  1:5509092/1‑150 (MQ=255)
  acgcCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTc                                                                                                                                                  <  2:6950068/150‑1 (MQ=255)
  acgcCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTc                                                                                                                                                  <  2:2332340/150‑1 (MQ=255)
    gcCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTcgc                                                                                                                                                <  2:1089905/150‑1 (MQ=255)
      cGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCgg                                                                                                                                              <  2:6768345/150‑1 (MQ=255)
          gCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCAC                                                                                                                                          >  2:8918324‑M1/1‑150 (MQ=255)
          gCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCAC                                                                                                                                          >  2:1319546‑M1/1‑150 (MQ=255)
          gCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCAC                                                                                                                                          <  1:6229859‑M1/150‑1 (MQ=255)
          gCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCAC                                                                                                                                          <  1:6229605‑M1/150‑1 (MQ=255)
            aaGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGC                                                                                                                                        >  1:2299759‑M1/1‑150 (MQ=255)
              gCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGGACGCCC                                                                                                                                      >  1:7174055‑M1/1‑150 (MQ=255)
               cGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCG                                                                                                                                     <  2:5771268‑M1/150‑1 (MQ=255)
                gtgtgAGTATCGACACCTTGATCAAGCCTTTCTGTCAGAAGCACTACGTGGTGGCCAAGCTGGTGCTGCCGAAAGAGCAGCATCTTAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGC                                                                                                                                    <  2:2311949‑M1/148‑1 (MQ=255)
                 ggTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCT                                                                                                                                   <  2:191685‑M1/150‑1 (MQ=255)
                 ggTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCT                                                                                                                                   <  1:330527‑M1/150‑1 (MQ=255)
                  gTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTATCGGCCCCTTCGCGGGCACGCCCGCTC                                                                                                                                  <  1:4150863‑M1/150‑1 (MQ=255)
                  gTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTC                                                                                                                                  <  1:4979564‑M1/150‑1 (MQ=255)
                  gTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTC                                                                                                                                  <  1:2706767‑M1/150‑1 (MQ=255)
                   tGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCGTGCCCGCgaa                                                                                                                                 <  2:6353399‑M1/150‑4 (MQ=255)
                   tGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCC                                                                                                                                 <  2:2025767‑M1/150‑1 (MQ=255)
                     aGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCGTCGCGGGCACGCCCGCTCCCA                                                                                                                               <  2:300332‑M1/150‑1 (MQ=255)
                      gTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCAC                                                                                                                              <  2:2000827‑M1/150‑1 (MQ=255)
                       tATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACA                                                                                                                             <  1:3750920‑M1/150‑1 (MQ=255)
                         tCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGG                                                                                                                           <  1:4648556‑M1/150‑1 (MQ=255)
                           gACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTA                                                                                                                         <  1:8027002‑M1/150‑1 (MQ=255)
                            acacCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAg                                                                                                                        >  2:5364201/1‑150 (MQ=255)
                              acCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgt                                                                                                                      <  2:4087385‑M1/150‑4 (MQ=255)
                               ccTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAgcgc                                                                                                                     <  1:4998957/150‑1 (MQ=255)
                               ccTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAgcgc                                                                                                                     <  1:4592117/150‑1 (MQ=255)
                                 ttGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgag                                                                                                                   >  2:6676688‑M1/1‑144 (MQ=255)
                                     gCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGcg                                                                                                               >  2:7244500/1‑149 (MQ=255)
                                       aaGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgc                                                                                                             >  1:5946168‑M1/1‑138 (MQ=255)
                                            ttACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTcgcgc                                                                                                        <  2:8370773/150‑1 (MQ=255)
                                             tACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgaga                                                                                                       <  2:2593705‑M1/150‑19 (MQ=255)
                                                  gCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGTAGCGCGGCGCGCACGCGCGCCGCt                                                                                                  >  2:5009528/1‑150 (MQ=255)
                                                  gCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCt                                                                                                  >  2:8190332/1‑150 (MQ=255)
                                                  gCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCt                                                                                                  <  1:4510406/150‑1 (MQ=255)
                                                   cAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGTCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTc                                                                                                 <  1:5239493/150‑1 (MQ=255)
                                                      aaGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGGAGCGCGTCGCTCTCGCGCGCAGCTcaca                                                                                              >  2:5912528/1‑150 (MQ=255)
                                                         ggcggcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATAc                                                                                           >  2:1014205‑M1/87‑150 (MQ=255)
                                                          aCTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCGCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAaacc                                                                                          <  1:6035620‑M2/150‑5 (MQ=255)
                                                           cTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAtacct                                                                                         >  1:5504946‑M2/1‑145 (MQ=255)
                                                            tACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAtaccgg                                                                                        <  2:8460688‑M2/150‑7 (MQ=255)
                                                            tACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCAAGCCCGCAACCACAGGGAGAGCGGCGCTCTCGCGCGCAGCGCACAtaccgg                                                                                        <  1:7543964‑M2/150‑7 (MQ=255)
                                                             gcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgt                                                                                       >  1:3377515‑M1/83‑150 (MQ=255)
                                                             gcggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgt                                                                                       >  1:1597449‑M1/83‑150 (MQ=255)
                                                             aCGTGGTGGCCAAGGGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgt                                                                                       <  2:4661788‑M1/150‑35 (MQ=255)
                                                             aCGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgt                                                                                       <  2:5946168‑M1/150‑35 (MQ=255)
                                                             aCGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAtacctgt                                                                                       >  2:4826181‑M2/1‑143 (MQ=255)
                                                             aCGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCGTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAtacctgt                                                                                       >  2:4272492‑M2/1‑143 (MQ=255)
                                                              cggatgtgatgctgccgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTTCGCGGGCACGCGCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgtg                                                                                      >  2:6286411‑M1/82‑150 (MQ=255)
                                                              cGTTGTGGCCAAGCGGGTGCTGCCGTAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAtacctgtg                                                                                      >  2:329117‑M2/1‑142 (MQ=255)
                                                              cGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAtacctgtg                                                                                      >  2:1735526‑M2/1‑142 (MQ=255)
                                                              cGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAtacctgtg                                                                                      >  1:3555433‑M2/1‑142 (MQ=255)
                                                               gtggtgGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCTTTCGCGGGCGCGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgtgg                                                                                     <  2:9066490‑M1/150‑37 (MQ=255)
                                                               gtggtgGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCTTTCGCGGGCGCGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgtgg                                                                                     <  1:3738565‑M1/150‑37 (MQ=255)
                                                               gtggtgGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgtgg                                                                                     <  1:6676688‑M1/150‑37 (MQ=255)
                                                               gtggtgGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgtgg                                                                                     <  1:9156691‑M1/150‑37 (MQ=255)
                                                               gtggtgGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAtacctgtgg                                                                                     >  2:2928708‑M2/1‑141 (MQ=255)
                                                               ggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                     >  1:6531549‑M1/81‑150 (MQ=255)
                                                               ggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                     >  2:7285601‑M1/81‑150 (MQ=255)
                                                               ggaggggatgctgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTCCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                     >  1:1751867‑M1/81‑150 (MQ=255)
                                                               ggaggggatgccgacgatggagggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTgg                                                                                     >  2:4385542‑M1/81‑150 (MQ=255)
                                                                tggtggCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgtggg                                                                                    <  2:7938181‑M1/150‑38 (MQ=255)
                                                                tggtggCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgtggg                                                                                    <  2:5507466‑M1/150‑38 (MQ=255)
                                                                   gggcCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAtacctgtgggagc                                                                                 >  2:3561937‑M2/2‑137 (MQ=255)
                                                                      ccAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgtgggagcggg                                                                              <  1:8373333‑M1/150‑44 (MQ=255)
                                                                        aaGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCGGGCCCGCgaa                                                                                                                                 >  1:8825792‑M1/1‑94 (MQ=255)
                                                                        aaGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCCCCCCCCCACCCACACGTAGCCCGGAGCACACGCCCGCCCCGCACAAACCCCTGGGGGcgcgcg                                                                            >  1:5771122/1‑150 (MQ=255)
                                                                        aaGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgtgggag                                                                                  >  2:5386430‑M1/1‑105 (MQ=255)
                                                                        aaGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgtgggag                                                                                  <  1:5386430‑M1/144‑40 (MQ=255)
                                                                         aGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGgcgc                                                                           <  1:3843055/150‑1 (MQ=255)
                                                                           gacgagggtgggggcgtggcggaggctgtgggagcgagggttggacaaaggcttgtagtacctgcaaaGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTGGCGCGCAGCTCACATACCTGTGGGAGCGGgcgcgc                                                                         >  1:3809325‑M1/69‑150 (MQ=255)
                                                                                 ctgcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCGCCCACCGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCGGTGGGAGCGGGCGCGCCCGCGa                                                                   >  2:8763112/1‑150 (MQ=255)
                                                                                  tgcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACCTACCCGTGGGAGCGGGCGCGCCCGCGaa                                                                  >  2:7117531/1‑150 (MQ=255)
                                                                                    gggggcgttgcggaggctgagggagcgagggttgatcaaaggcttgtagtacctgcaaaGGCCCTTTCGCGGGCGCGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCtgt                                                                                       <  1:29544‑M1/68‑1 (MQ=255)
                                                                                    ccGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGGCGCTCGCGCGCGCAGCTCACATACCTGTGGGGGCGGGCGCGCCCGCGAAAg                                                                >  2:2148407/1‑150 (MQ=255)
                                                                                      gAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAggg                                                              >  2:7437950/1‑150 (MQ=255)
                                                                                      gAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACAAACCTGTGGGAGCGGGCGCGCCCGCGAAAagg                                                              >  2:5277205/1‑150 (MQ=255)
                                                                                       aaaGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCTCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGc                                                             <  2:2115464/150‑1 (MQ=255)
                                                                                       aaaGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCCCGCCCACCCGGAGGGCCTCGCGGTCGCGCGCAGCGCCCCGACCAGGGGGGGCGGGGGCGGCCCCGGAAgggg                                                             >  1:549696/1‑149 (MQ=255)
                                                                                           agcagcagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCttt                                                         <  2:4754893/150‑1 (MQ=255)
                                                                                           agcagcagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACCCCCGCTCCCCCCGGTAGCGCGGCGCTCTCGCGCGCAGCGCACCTACCGGGGGGCGCGGGGGCGCCCGCGGAAGGGccggg                                                         >  2:2628345/1‑147 (MQ=255)
                                                                                           agcagcagcTGAGTCTGGCAAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCttt                                                         <  2:1486303/150‑1 (MQ=255)
                                                                                              agcagcTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCGCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACCTACCGGTGGGAGCGGGCGCGCCCCCGAAACGGCCGTTGCa                                                      >  2:3207226/1‑150 (MQ=38)
                                                                                               gcagcTGAGTCTGGCCATGCGTTATGCTGATCGGCGCTATCCGGTACCGGCCCCTGCGCGGGCACGCCCGCTCCCACCGGTCGCGCGACGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAg                                                     <  1:8659297/150‑1 (MQ=38)
                                                                                                        tCTGGCCAAGCTTTCAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACCGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGTCTTTGCAGGTACTACaa                                            <  1:2568340/150‑1 (MQ=255)
                                                                                                            gCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCAGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCt                                        <  1:5656044/150‑1 (MQ=255)
                                                                                                             ccAAGCTTTATGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCCCGCCCGCTCCCACCGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCtt                                       <  1:6545451/150‑1 (MQ=255)
                                                                                                              cAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCttt                                      <  1:6089816/150‑1 (MQ=255)
                                                                                                                aGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGa                                    <  2:7555152/150‑1 (MQ=255)
                                                                                                                 gCGTAAAGCTGATCGGCGCTAGCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGAt                                   <  1:7904773/150‑1 (MQ=255)
                                                                                                                   gTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACCGGGAGCGCGTCGCGCTCGCGCGCAGCTCACCTACCGGGGGGCGCGGGCGCGCCCGCCCAAGGGCCGTGGCGGGGACTCCCAGCCCTTGATcc                                 >  2:6562038/1‑149 (MQ=255)
                                                                                                                     tacgCTTATCGGCGCTATCCGGTCCCGGCCCCTTCGCGGGCACGCCCGCTCCCACCGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAAc                               <  1:3195099/147‑1 (MQ=255)
                                                                                                                      acgCTGATCGGCGCTAGCCGGTACCGGCCCCGTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCGCGCGCGCAGCTCACATACCTGGGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAAcc                              <  2:8633250/148‑1 (MQ=255)
                                                                                                                               ggCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAtgtgagctgcgcgcgagagcgacgcgctacctgt                                                                                       <  1:2012545‑M1/84‑35 (MQ=255)
                                                                                                                                    tATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGGCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGc                <  1:1236906/150‑1 (MQ=255)
                                                                                                                                                  CCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGcccc  >  2:7174039‑M1/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCACGCCGTGGCAAGCGGGTGAGTATCGACACCTTGAGCAAGCCTTACTGGCAGAAGCACTACGTGGTGGCCAAGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCC  >  NC_002947/5818673‑5818968

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: