Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,575,407 A→G 6.0% intergenic (+24/‑164) pcaF‑I → / → pcaT beta‑ketoadipyl‑CoA thiolase/alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,575,4070AG 6.0% 612.2 / 5.6 193intergenic (+24/‑164)pcaF‑I/pcaTbeta‑ketoadipyl‑CoA thiolase/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base A (67/110);  new base G (10/2);  total (79/114)
Fisher's exact test for biased strand distribution p-value = 3.96e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 200 total aligned reads displayed.

ATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCA  >  NC_002947/1575258‑1575556
                                                                                                                                                     |                                                                                                                                                     
aTGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGa                                                                                                                                                       >  1:4557495/1‑150 (MQ=255)
 tGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGaa                                                                                                                                                      <  1:2622174/150‑1 (MQ=255)
 tGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGaa                                                                                                                                                      <  2:4665251/150‑1 (MQ=255)
   ttcttcgccggtgtaataggttgaggtcatgacgggtttactcctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCACTTCGCGGGGGAGc                                                                                                                                                    <  1:393593‑M1/19‑1 (MQ=255)
   aGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGc                                                                                                                                                    >  1:5071556/1‑150 (MQ=255)
   aGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGc                                                                                                                                                    >  1:7160842/1‑150 (MQ=255)
    tcttcgccggtgtaataggttgaggtcatgacgggtttactcctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGcc                                                                                                                                                   >  2:1719972‑M1/131‑150 (MQ=255)
    tcttcgccggtgtaataggttgaggtcatgacgggtttactcctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGcc                                                                                                                                                   <  2:6772264‑M1/20‑1 (MQ=255)
    gCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGcc                                                                                                                                                   <  2:332108/150‑1 (MQ=255)
    gCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGcc                                                                                                                                                   <  2:4084299/150‑1 (MQ=255)
     cttcgccggtgttataggttgaggtcatgacggctttactcctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccctatgtctcacaggagaGGCCTCTTCGCGGGGAAGccc                                                                                                                                                  >  1:2250912‑M1/130‑150 (MQ=255)
      ttcgccggtgtaataggttgaggtcatgacgggtttactcctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCg                                                                                                                                                 <  2:3375603‑M1/22‑1 (MQ=255)
      ggTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCg                                                                                                                                                 <  1:7108972/150‑1 (MQ=255)
         gccggtgtaataggttgaggtcatgacgggtttactcctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTc                                                                                                                                              >  2:7434877‑M1/126‑150 (MQ=255)
         gCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTc                                                                                                                                              <  1:1705545/150‑1 (MQ=255)
           aCGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCa                                                                                                                                            >  2:5686167/1‑150 (MQ=255)
           aCGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCa                                                                                                                                            <  1:5187608/150‑1 (MQ=255)
            cGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTccca                                                                                                                                           <  1:2543601/150‑3 (MQ=255)
             gCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCCCACCAAc                                                                                                                                          >  1:4159583/1‑150 (MQ=255)
              ccTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGTAGCCCGCGCCaaaa                                                                                                                                         >  1:5911729/1‑148 (MQ=255)
              ccTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACa                                                                                                                                         <  1:4932664/150‑1 (MQ=255)
                tGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAgg                                                                                                                                       <  1:845200/150‑1 (MQ=255)
                   tACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGACGCCCGCGCCAACAgggct                                                                                                                                    >  1:5666159/1‑148 (MQ=255)
                     ctgactgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGc                                                                                                                                                    >  2:620170/1‑132 (MQ=255)
                       gactgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCAc                                                                                                                                <  1:121160/150‑1 (MQ=255)
                       gactgCGATGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCAc                                                                                                                                <  2:5210017/150‑1 (MQ=255)
                          tgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCa                                                                                                                             >  1:6422607/1‑150 (MQ=255)
                          tgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACACGGATACACCCa                                                                                                                             >  1:2954626/1‑150 (MQ=255)
                              ttGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCa                                                                                                                         <  1:2102354/150‑1 (MQ=255)
                               tGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCAc                                                                                                                        <  1:4735139/150‑1 (MQ=255)
                                 cACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCt                                                                                                                      <  2:6422607/150‑1 (MQ=255)
                                  acgggtttactcctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtg                                                                                                                     >  1:3301742‑M1/101‑150 (MQ=255)
                                   cgggtttactcctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCCACAGGGATCCACGCAATCACCtgt                                                                                                                    >  1:6203916‑M1/100‑150 (MQ=255)
                                   cgggtttactcctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCCACAGGGATCCACGCAATCACCtgt                                                                                                                    >  1:7217801‑M1/100‑150 (MQ=255)
                                   cgggtttactcctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                    >  1:2336244‑M1/100‑150 (MQ=255)
                                   ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATc                                                                                                                          >  2:1058813/1‑144 (MQ=255)
                                   ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATc                                                                                                                          <  1:1058813/144‑1 (MQ=255)
                                   ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                    <  1:5325223/150‑1 (MQ=255)
                                   ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                    <  1:4726572/150‑1 (MQ=255)
                                   ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                    <  1:4115500/150‑1 (MQ=255)
                                   ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                    >  2:5452268/1‑150 (MQ=255)
                                   ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                    >  2:2414157/1‑150 (MQ=255)
                                   ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                    <  1:1928953/150‑1 (MQ=255)
                                    cAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgtg                                                                                                                   <  1:3206048/150‑1 (MQ=255)
                                    cAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgtg                                                                                                                   <  1:2061608/150‑1 (MQ=255)
                                         ggAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGc                                                                                                              <  2:7209336/150‑1 (MQ=255)
                                         ggAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCGCCAACAGGGATCCACGCAATCACCTGTGGGAGc                                                                                                              <  2:3587040/150‑1 (MQ=255)
                                          gagaAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCCCGCCAAAAAGGATCCACCCAATAACCCGGGGGGGCg                                                                                                             >  1:1905248/1‑150 (MQ=255)
                                          gagaAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                             >  2:4843511/1‑150 (MQ=255)
                                          gagaAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                             <  2:859730/150‑1 (MQ=255)
                                          gagaAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                             <  1:1167138/150‑1 (MQ=255)
                                           agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCgg                                                                                                            <  2:5973095/150‑1 (MQ=255)
                                           agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                             >  1:205774/1‑149 (MQ=255)
                                            gaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCggg                                                                                                           <  1:1255483/150‑1 (MQ=255)
                                             aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                          <  2:4785310/150‑1 (MQ=255)
                                             aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                          <  2:1751370/150‑1 (MQ=255)
                                             aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                          <  1:5753634/150‑1 (MQ=255)
                                              ctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGc                                                                                                                                                    <  2:20846‑M1/19‑1 (MQ=255)
                                              ctggagggtcctaggtagttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCgg                                                                                                            <  1:764956‑M1/59‑1 (MQ=255)
                                               gagTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGAACCACGCAAACACCAGTGGGAGCGGGCtt                                                                                                        >  1:6772264/1‑150 (MQ=255)
                                                agTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTc                                                                                                       <  1:4697426/150‑1 (MQ=255)
                                                  tGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCGCCTGTGGGAGCGGGCTTccc                                                                                                     >  2:996875/1‑150 (MQ=255)
                                                  tGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTccc                                                                                                     <  1:2394444/150‑1 (MQ=255)
                                                    gCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                             >  1:1302337/1‑140 (MQ=255)
                                                    gCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                             <  2:1302528/140‑1 (MQ=255)
                                                          ccAAGGCCCAGCCCACCATGTGTGTGGGTGACGGCCAAGGACTGGCGTTGGCCATCGAGCGGCATTGACCTGTACCGGCCTCTTCGCGGCGAAGCCCGATCCAACACGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAaga                                                                                             <  2:7176005/149‑1 (MQ=255)
                                                            aaGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCGTCGCGGGGGCGCCCCCGCCAACAGGGATCCACCCCAACACCAGTGGGCGCGGGCGGCCCCCCCCAGAgg                                                                                           >  1:3824485/1‑150 (MQ=255)
                                                            aaGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCGGGACCGGCCTCGTCGCGGGGAAGCCCGCTCCACCAGGGACCCACGCAATCACCGGTGGGAGCGGGCTTCCCCGCGAAGAgg                                                                                           <  2:6594158/150‑1 (MQ=255)
                                                             aGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGc                                                                                          <  2:2954714/150‑1 (MQ=255)
                                                                ttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTgg                                                                                                                  >  2:1506599‑M1/71‑123 (MQ=255)
                                                                ttctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTgg                                                                                                                  <  1:1506599‑M1/53‑1 (MQ=255)
                                                                  ctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGGAGCCCGCTCCaa                                                                                                                                           >  2:2247963‑M1/69‑96 (MQ=255)
                                                                  ctagacatacctagtaacacactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCaa                                                                                                                                           <  1:5998295‑M1/28‑1 (MQ=255)
                                                                  cTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACCCAAACACCAGTGGGAGCCGGCGTACCCCCGAAGAGGCCtctc                                                                                     >  1:4833113/1‑150 (MQ=255)
                                                                      cGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCAGCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCtgt                                                                                 <  2:5800130/150‑1 (MQ=255)
                                                                        aCCTTGTGTGTGGGTGTCGGCAAAGTTCTGGCGTTGGAAATCGAGCGGTTTTGACCTGTACCGGCCTATTCGCGGGGAACCCCGCTCCAACAGGCATCCACGAAATCACCTGTGGTAGCGGGATTCCCCGCGAAGAGGCCTCTCCTGTga                                                                               <  1:6074659/150‑1 (MQ=255)
                                                                        aCCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTCGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCCCACCAACCGCGATCCACCCAAACACCACTGGGAGCGGGCTTCCCCGCGAAGACGCCTCTccactga                                                                               >  2:4130242/1‑145 (MQ=255)
                                                                              tgtgtgGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCCCCCCAACCGGGTTCACCGCAATCACCAGTGGGAGCGGGCGACCCCGCGAACAGGCCACTCCAgagagaaaaa                                                                         >  2:393593/1‑146 (MQ=255)
                                                                                    ggTGTCGGCGAAGGTCGGGGGTTGGCAATCGCGGGGGTTTGACCTGTACCGGCCTCTTCGGGGGGTAGCCCGCTCCAACAGGGATCCACGCAATCACGTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGccc                                                                   <  2:7001355/150‑1 (MQ=255)
                                                                                    ggTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCGTCGCGGGGGAGCCCCCGCCAACACGGGGACCCCCAAAAACCAGTGGGGGCGGGCTTACCCCCCCACAGGCCACAACAGGGAGACAAAAAGccc                                                                   >  1:6886232/1‑150 (MQ=255)
                                                                                    acactcgttccgctatcgggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCgg                                                                                                            <  2:1189003‑M1/59‑1 (MQ=255)
                                                                                          ggCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGGGAGACATAAGGCCCATTGcc                                                             >  2:7306471/1‑150 (MQ=255)
                                                                                           gcaaaGGTCTGGCGTGGGCCATCGCGGGGTTTTGACCTGTACCGGCCACTGCGCGGGGAAGCCCGCTCCAACAGGCATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCa                                                            <  2:1404466/147‑1 (MQ=255)
                                                                                           gCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCGCGCCAACAGGGATCCACGCAATCACCGGTGGGAGCGGGCTTCCCCCCGCAGAGGCCTCTCCGGGGAGACATAAGGCCCCATGCCa                                                            >  2:1991119/1‑150 (MQ=255)
                                                                                                     gggctggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGc                                                                                                                                                <  2:3183352‑M1/23‑1 (MQ=255)
                                                                                                         tggcaatgggccttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGc                                                                                                              <  2:4019673‑M1/57‑1 (MQ=255)
                                                                                                          tGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCGCCACCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTGCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAAc                                             <  2:6277319/150‑1 (MQ=255)
                                                                                                              cATCGAGCGGGTTGGACCTGTACCGGCCTCTTCGCGGGGATGCCCGCTCCAACCGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgt                                         <  2:3718721/150‑1 (MQ=255)
                                                                                                              atgggccttatgtctcataggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgtg                                                                                                                   <  2:1299799‑M1/52‑1 (MQ=255)
                                                                                                                tCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgtgt                                       <  2:4557495/150‑1 (MQ=255)
                                                                                                                 cGAGGGGGTTTGACCTGTACCGGCGTATTCGCGGGGTAGCCCGCTCCAACAGGGATCATCGCAATCACCTGTGGTAGCGGGCTTCCCCGCGAAGAGGGCTCTCCTTTGATACTTAAGGCCCATTGCCAGCCCGATAGCGGAACGAgtgtg                                      <  2:517257/150‑1 (MQ=255)
                                                                                                                     ttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                             >  1:3626178‑M1/18‑75 (MQ=255)
                                                                                                                     ttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                             <  2:3626078‑M1/58‑1 (MQ=255)
                                                                                                                     ttatgtctcacaggagaGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                             >  2:265816‑M1/18‑75 (MQ=255)
                                                                                                                                  tACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAg                     <  2:5239289/150‑1 (MQ=255)
                                                                                                                                           CTTCGCGGCGAACCCCGCTCCAACAGCCATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGAGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTa            >  2:5911729‑M1/1‑150 (MQ=255)
                                                                                                                                            TTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAg           <  2:3659670‑M1/150‑1 (MQ=255)
                                                                                                                                              CGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCCCCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGa         >  1:3058258‑M1/1‑150 (MQ=255)
                                                                                                                                                   GGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGgaccccc    <  1:4629508‑M1/150‑3 (MQ=255)
                                                                                                                                                     AAGCACGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCa  <  2:6001158‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGAGCGGTGCACGCCTGGTACTGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCA  >  NC_002947/1575258‑1575556

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: