Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,575,429 A→C 8.1% intergenic (+46/‑142) 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,4290AC 8.1% 443.7 / 4.8 158intergenic (+46/‑142)pcaF‑I/pcaTbeta‑ketoadipyl‑CoA thiolase/alpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base A (56/84);  new base C (9/8);  total (65/93)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

TGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACCTC  >  NC_002947/1575280‑1575578
                                                                                                                                                       |                                                                                                                                                     
tgactgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCa                                                                                                                                                       >  1:430530/1‑150 (MQ=255)
 gactgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCAc                                                                                                                                                      <  1:508505/150‑1 (MQ=255)
 gactgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAAGGATacac                                                                                                                                                      >  1:6446395/1‑150 (MQ=255)
    tgCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCa                                                                                                                                                   >  2:4604190/1‑150 (MQ=255)
      cGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCCCGCCAACCGGGATCCACCCaaa                                                                                                                                                 >  1:3647837/1‑149 (MQ=255)
       gTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCACTc                                                                                                                                                >  2:5550894/1‑150 (MQ=255)
        ttGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCa                                                                                                                                               <  2:6534798/150‑1 (MQ=255)
         tGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCATCAGGGTTCCACGCAATCAc                                                                                                                                              <  1:5908846/150‑1 (MQ=255)
           cACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCt                                                                                                                                            <  1:1556501/150‑1 (MQ=255)
             ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                                          <  1:3563171/150‑1 (MQ=255)
             ccAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgt                                                                                                                                          <  1:3562959/150‑1 (MQ=255)
              cAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCtgtg                                                                                                                                         <  1:5710567/150‑1 (MQ=255)
               aGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTgg                                                                                                                                        <  1:1026489/150‑1 (MQ=255)
               aGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTgg                                                                                                                                        <  2:6340629/150‑1 (MQ=255)
                 cTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGa                                                                                                                                      >  2:4814488/1‑150 (MQ=255)
                 cTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGAGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGa                                                                                                                                      <  2:2646565/150‑1 (MQ=255)
                    gagaAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCGCGCCAAAAGGGGTCCACGCACACACCCGGGGGGGCg                                                                                                                                   >  1:6563806/1‑150 (MQ=255)
                    gagaAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCg                                                                                                                                   >  2:1584316/1‑150 (MQ=255)
                     agaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCgg                                                                                                                                  <  2:4873794/150‑1 (MQ=255)
                      gacgcgTGGCGGTCGCAAGGGCCTGGCGACCCGGTGTGTGGGGGTCGGCCAAGGTCTGGCGTTGGCCATCGCGCGGGTTTGACCTGGACCGGCCTCTTCGCGGGGACGCCCGCTCCAACCGGGCTCCCCGCAATCACCTGTGGGAGCggg                                                                                                                                 <  2:368646/145‑1 (MQ=255)
                      gaagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGTAGCCCGCGCCCACAGGGGTCCACGCAATCACCgggggggggggg                                                                                                                                 >  1:4256400/1‑150 (MQ=255)
                       aagaGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                                                <  1:3089340/150‑1 (MQ=255)
                        agagTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCt                                                                                                                               <  1:2842390/150‑1 (MQ=255)
                        agagTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCACCAACAGGGATCCCCCCAAACACCAGTGGGGGCggggt                                                                                                                               >  1:5451674/1‑148 (MQ=255)
                         gagTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGCAGCCCGCGCCAACCCGGATCCCCGCAATCACCGGTGGGAGCGGGcgt                                                                                                                              >  2:5155764/1‑148 (MQ=255)
                         gagTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCtt                                                                                                                              <  2:3782730/150‑1 (MQ=255)
                              gCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCcg                                                                                                                         >  2:1285509/1‑150 (MQ=255)
                              atggTCGCCAGGCCCTTGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAAAAGGGGTCCACGCAATCACCAGTGGGAGCGGGCGTCCCcg                                                                                                                         >  1:1417171/3‑150 (MQ=255)
                                 gTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGa                                                                                                                      <  1:5511724/150‑1 (MQ=255)
                                  tCGCACGGGCCTGGCGACCTTGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTTTACCGGCCTCTTCGCGGGGAAGCCCGCTCCCACAGGTGATTGCGTGGATCCCTGTTGGAGCGGGCTTCCCCGCGaa                                                                                                                     <  1:2685597/150‑1 (MQ=17)
                                  tCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCCCGCCCACCGGGGTCCCCCCCAAAACCAGGGGGGGCGGGCGTGCCcgccga                                                                                                                     >  1:2986633/1‑147 (MQ=255)
                                    gCACGGGCCTGGCGACCTTGTGTGTGGGTGTCGGCCACGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTCCCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTCCCCCGCGAaga                                                                                                                   <  1:5499905/150‑1 (MQ=255)
                                    gCACGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGGGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGTATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAaga                                                                                                                   <  2:4163401/150‑1 (MQ=255)
                                    gCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACCGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAaga                                                                                                                   <  2:3844744/150‑1 (MQ=255)
                                      aaGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGTCCCCCGCGCCAACAGGGGAACCCCCCCACAACAGGGGGCGCGGGCTTCCCCCCGCagacg                                                                                                                 >  1:5312682/1‑148 (MQ=255)
                                      aaGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCGCTCCAACAGGGCTCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAgg                                                                                                                 >  1:4489459/1‑150 (MQ=255)
                                      aaGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCGCCAACACGGATCCACGCCATCACCGGTGGGAGCGGGCGTCCCCGCGACGAgg                                                                                                                 >  1:5551376/1‑150 (MQ=255)
                                                 gACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCtgtg                                                                                                      >  2:5525068/1‑150 (MQ=255)
                                                 gACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCAGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCtgtg                                                                                                      >  2:500212/1‑150 (MQ=255)
                                                  aCCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTga                                                                                                     <  2:430530/150‑1 (MQ=255)
                                                     aTGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAAACACCTGTGGGAGCGGGCTTCCCCCCGAAGAGGCCTCTCCTGTGAGAc                                                                                                  >  1:4469543/1‑150 (MQ=255)
                                                       gtgtgtgGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACAt                                                                                                >  2:4483870/1‑150 (MQ=255)
                                                         gtgtgGGTGGCGGCCACGGTCGGGCGTGGGCCATCGTGGGGGTTTGACCTGTACGGGCCTCTTCGGGGGGACGCCCGCTCCAACCGGGATCCACGCAATCCCCGGTGGGAGGGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATaa                                                                                              <  2:375707/150‑1 (MQ=255)
                                                          tgtgGGTGTCGGCCAAGGTCTGGGGTTGGCCTCCGAGCGGGTTTGACCTGTACCGGCCTATTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACGTGTGGGAGCGGGATTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAg                                                                                             <  2:1949147/150‑1 (MQ=255)
                                                            tgGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCCCACCAAAAAGGCTACACCCAATAACCAGTGGGGTCGGGCTTACCCCCGAAGAGGCCACTCCAGTGAGAAAAAAGGc                                                                                           >  1:6600483/1‑150 (MQ=255)
                                                             gggTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCCCGCCAAAACGGATAACCGCAAACACCGGGGGGCGCGGGCGTCCCCCCGAAGAGGCCGCGCCGGTGGGACATAAGGcc                                                                                          >  2:2847414/1‑150 (MQ=255)
                                                               gtgtCGGCCATGGTCTGGCGTTGTCCATCGAGGGGTTTTGACCTGTACCGGCCTCTTGGCGGGTAAGCCCGCTCCAACTGTGATCCACGAAAGCACGTGTGGGTGGGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCa                                                                                        <  2:213611/150‑1 (MQ=255)
                                                                 gtCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCAtt                                                                                      <  2:2794608/150‑1 (MQ=255)
                                                                     gCCAAGGTCTGGCGTTGGCCATCGCGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGTAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCa                                                                                  <  1:6217944/150‑1 (MQ=255)
                                                                      ccAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCGCGCCAACCGGGGTGCCCCCCAACACCTGGGGGGGCGGGGGTGCCCGCGCAGGGGCCTCGCCGGGGGGGCCGAAGGCCCCTTggccg                                                                                 >  1:4904174/1‑146 (MQ=255)
                                                                       cAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGCAGCCCCCGCCAACAGGGCTCCCCCCAATAAACAGTGGGTGCGGGCATCCCCCCGAAGAGGCCTCGCCTGTGAGACCTAAGGCCCAGTGCCacc                                                                                >  2:844848/1‑148 (MQ=255)
                                                                        atggTGTGGCGTTGGCCATCGCGGGGGTTTGCCCTGTACGGGCGTCTTGGCGGGGAAGCCCGCTCCAACTGGGATCCACGCAATCACGTGTGGGTGCGGGCTTCCCGGCGTAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGcc                                                                               <  2:5885925/148‑1 (MQ=255)
                                                                        aaGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGTTCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGcc                                                                               <  2:1793050/150‑1 (MQ=255)
                                                                             cTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATa                                                                          <  2:6127122/150‑1 (MQ=255)
                                                                              tGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAg                                                                         >  2:5606617/1‑150 (MQ=255)
                                                                               ggCGTTGGCCATGGCGGGGGTTTGACCGGTACCGGCCTCTTCGCGGGGACGCCCGCTCCAACCGGGATCCCCGCACTCACCGGGGGGCGCGGGCTGCCCCGCGAAGAGGCCTCTCCTGTGAGACCTAAGGCCCATTGCCAGCCCGCTAGc                                                                        <  2:1882089/150‑1 (MQ=255)
                                                                               ggCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGTAAGCCCGCTCCAACAGGGATCCACGAAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGc                                                                        <  1:6804198/150‑1 (MQ=255)
                                                                               ggCGTGGGCCAGCGCGGGGGTTGGACCGGTCCCGGCCTGTGCGGGGGGACGCCCGCGCCACCCGGGCGCCCCGCAACCCCCGGGGGGCGCGGGCTGCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGc                                                                        <  2:1882141/150‑1 (MQ=255)
                                                                                     ggCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGAGCCCGCGCCCACAGGGGTCCAAGCCAAAAACGGGGGGGGCGGGCGTCCCCCCCGAGAGGCCACTCCGGTGGGACAAAAGGGCCAATGGCCGCCCCATAGCGGaaag                                                                  >  2:2869252/1‑148 (MQ=255)
                                                                                          tCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGGACCCCCCGCCCACACGGATCCAAGCCAACACCGGGTGGAGCGGGCATCCCCCCGAAGAGGCCACTCCAGTGAGTCAAAAGGCCCAATTCCAGCCCCATAGCGGAACGAggtt                                                             >  2:2289027/1‑147 (MQ=255)
                                                                                           cgcgCGGGTTTGCCCTGTCCCGGCCTCTTCGCGGGGACGCCCGCTCCAACCGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAgtgtg                                                            <  2:6303215/147‑1 (MQ=255)
                                                                                            gAGGGGGTTTGTGCTGTACCGGCCTTTTCGCGGGGTTGCCGGCTCAAAGCGGGATCCTCGAATTCAGGTGTGGGTGCGGGCTTCCGCTCGAATTGGACTCTGCTTTGAGACATTAGGCCCATTGCCAGCCCGATAGCGGAACgtgtgtgt                                                           <  2:3159907/150‑1 (MQ=255)
                                                                                                gggTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACt                                                       <  1:2808440/150‑1 (MQ=255)
                                                                                                 ggTTTGACCTGTACGGGCCTCTTCGGGGGGAAGCCCGCTCAAACTGGGATCCACGAAATCACCTTTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTa                                                      <  2:1701599/150‑1 (MQ=255)
                                                                                                 ggTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTa                                                      <  2:6726798/150‑1 (MQ=255)
                                                                                                 ggTTTGACCGGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTa                                                      <  2:6543326/150‑1 (MQ=255)
                                                                                                  gTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAg                                                     <  2:6646603/150‑1 (MQ=255)
                                                                                                  gTTTGACCGGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAg                                                     <  2:1417171/150‑1 (MQ=255)
                                                                                                      gcccTGTCCCGGCCTCTTCGCGGGGACGCCCGCTCCAACAGGGATCCACGCAATCCCCTGTGGGTGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTAt                                                 <  1:6280789/148‑1 (MQ=255)
                                                                                                        cgggTACGGGCCAGTTCGCGGGGACGCCCGCCACTACAGTTACCCACGCAATCACCTGTGGGAGGGGGCTTCCCCGGGATGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAAAGAGTGTGTTACTAGGTATGt                                               <  1:3401493/147‑1 (MQ=255)
                                                                                                            tcccGGCCTCTGCGGGGGCACCCCCGCACAAACCGGTACCCACCAAATCACGTGGGGGCGCGGGCTCCCCCGCTAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAg                                           <  1:1466670/148‑1 (MQ=255)
                                                                                                              ccGGCCTCTTCGGGGGGACCCCCGCACAAACCGGTAACCACGAAACCCCGTGGGGGCGGGGGTTCCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGaa                                         <  1:1842613/150‑1 (MQ=255)
                                                                                                              ccGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCGGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGaa                                         >  2:5015029/1‑150 (MQ=255)
                                                                                                                    tctTCGCGGGGAAGCCCGCTCCCACAGGTGATTGCGTGGATCCCTGTTGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGGGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTGTGTCTAGAACTACCt                                   >  2:662144/1‑150 (MQ=17)
                                                                                                                    tctTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGCTAGCGGCACGAGTGTGTTACTAGGTATGTCTAGAACTACCt                                   >  2:6684518/1‑150 (MQ=255)
                                                                                                                    tctTCGCGGGGAAGCCCGCTCCAACAGGCATCCCCGCAATCACCTGTGGGAGCGGGCTCCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCt                                   <  2:3037642/150‑1 (MQ=255)
                                                                                                                       tCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAgg                                >  2:5238015/1‑150 (MQ=255)
                                                                                                                        cgcgGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGa                               <  1:3866893/150‑1 (MQ=255)
                                                                                                                         gcgGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGc                                                                                                                                >  1:2337901/1‑52 (MQ=255)
                                                                                                                               aaGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGGATGTCTAGAACTACCTAGGACCCTCCa                        >  1:4665897/1‑150 (MQ=255)
                                                                                                                                aGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAg                       <  2:3875964/150‑1 (MQ=255)
                                                                                                                                aGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAg                       <  2:4172421/150‑1 (MQ=255)
                                                                                                                                  cccGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                     >  1:511170/1‑150 (MQ=255)
                                                                                                                                  cccGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                     <  2:2842575/150‑1 (MQ=255)
                                                                                                                                  cccGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                     <  2:4579831/150‑1 (MQ=255)
                                                                                                                                   ccGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAggct                    >  2:6759755/1‑148 (MQ=255)
                                                                                                                                   ccGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAg                    >  1:4887133/1‑150 (MQ=255)
                                                                                                                                    cGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGt                   >  2:1027508/1‑150 (MQ=255)
                                                                                                                                         cAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAAcc              >  2:2778431/1‑150 (MQ=255)
                                                                                                                                          aaCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAAccc             >  1:2117020/1‑150 (MQ=255)
                                                                                                                                          aaCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAAccc             >  2:4886359/1‑150 (MQ=255)
                                                                                                                                          aaCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAAccc             >  1:1237794/1‑150 (MQ=255)
                                                                                                                                          aaCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTAGTTCTAGACATACCTa                                  >  1:6432871/1‑129 (MQ=255)
                                                                                                                                          aaCAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTAGTTCTAGACATACCTa                                  <  2:6432871/129‑1 (MQ=255)
                                                                                                                                                ggATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                     >  2:6185292/1‑137 (MQ=255)
                                                                                                                                                ggATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGa                     <  1:6185292/137‑1 (MQ=255)
                                                                                                                                                ggATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCAt        >  1:2384193/1‑150 (MQ=255)
                                                                                                                                                    accACGCAAACCCGGGTGGGCGGGGGCTGCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGAc     <  1:1749373/149‑1 (MQ=255)
                                                                                                                                                       aCGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACCTc  >  2:5441066/1‑150 (MQ=255)
                                                                                                                                                       aCGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACCTc  >  2:77194/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TGACTGCGTTGCACCAGCTGGAGAAGAGTGGCGGTCGCAAGGGCCTGGCGACCATGTGTGTGGGTGTCGGCCAAGGTCTGGCGTTGGCCATCGAGCGGGTTTGACCTGTACCGGCCTCTTCGCGGGGAAGCCCGCTCCAACAGGGATCCACGCAATCACCTGTGGGAGCGGGCTTCCCCGCGAAGAGGCCTCTCCTGTGAGACATAAGGCCCATTGCCAGCCCGATAGCGGAACGAGTGTGTTACTAGGTATGTCTAGAACTACCTAGGACCCTCCAGGAGTAAACCCGTCATGACCTC  >  NC_002947/1575280‑1575578

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: