Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,553,942 A→T 8.2% intergenic (+70/+111) pydA → / ← rutR NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,553,9420AT 8.2% 257.7 / 9.0 82intergenic (+70/+111)pydA/rutRNADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family
Reads supporting (aligned to +/- strand):  ref base A (18/55);  new base T (3/4);  total (22/60)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTG  >  NC_002947/4553795‑4554089
                                                                                                                                                   |                                                                                                                                                   
aTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGc                                                                                                                                                   <  1:3119495/150‑1 (MQ=255)
aTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGc                                                                                                                                                   <  2:1921958/150‑1 (MQ=255)
aTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGAACGGGCCGCACAGc                                                                                                                                                   <  2:1921930/150‑1 (MQ=255)
   gAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACCGCGGc                                                                                                                                                <  1:4467175/150‑1 (MQ=255)
   gAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGc                                                                                                                                                <  1:3465010/150‑1 (MQ=255)
    aGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGcc                                                                                                                                               <  1:1852620/150‑1 (MQ=255)
     ggCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGccc                                                                                                                                              <  2:1831978/150‑1 (MQ=255)
      gCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGTGGcccc                                                                                                                                             <  2:3833647/150‑1 (MQ=255)
        acaGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcag                                                                                                                                           <  1:517636/150‑1 (MQ=255)
        acaGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcag                                                                                                                                           <  2:1836080/150‑1 (MQ=255)
         caGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTGTGCGGcccgatc                                                                                                                                          <  1:507480/150‑5 (MQ=255)
          aGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcagca                                                                                                                                         <  1:2830718/150‑1 (MQ=255)
          aGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcagca                                                                                                                                         <  1:2304029/150‑1 (MQ=255)
          aGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcagca                                                                                                                                         <  1:2297507/150‑1 (MQ=255)
          aGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcagca                                                                                                                                         <  1:4138950/150‑1 (MQ=255)
           ggATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTGTGCGGcccgatcgc                                                                                                                                        <  2:3389107/150‑7 (MQ=255)
           ggATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTGTGCGGcccgatcgc                                                                                                                                        <  2:1859584/150‑7 (MQ=255)
            gATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGt                                                                                                                                       <  2:3270210/150‑1 (MQ=255)
              tACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCt                                                                                                                                     <  1:4084551/150‑1 (MQ=255)
               aCCGGCAAGCCGTACCTGATCTGGACGCAGGATCCGCGTAACCCCTACCGTGTGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCtt                                                                                                                                    <  1:3883726/150‑1 (MQ=255)
               aCCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCtt                                                                                                                                    <  1:3883812/150‑1 (MQ=255)
               aCCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCtt                                                                                                                                    <  1:3883727/150‑1 (MQ=255)
               aCCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCtt                                                                                                                                    <  1:2366165/150‑1 (MQ=255)
                ccGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCttt                                                                                                                                   <  2:3213010/150‑1 (MQ=255)
                ccGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCttt                                                                                                                                   <  2:1889310/150‑1 (MQ=255)
                 cGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTg                                                                                                                                  <  2:2390712/150‑1 (MQ=255)
                  ggCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGc                                                                                                                                 >  2:2110853/1‑150 (MQ=255)
                          gTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATAc                                                                                                                         <  1:2640584/150‑1 (MQ=255)
                          gTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATAc                                                                                                                         <  1:2640574/150‑1 (MQ=255)
                             ccTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGa                                                                                                                      <  1:3073688/150‑1 (MQ=255)
                                gAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACAt                                                                                                                   <  1:3252121/150‑1 (MQ=255)
                                gAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACAt                                                                                                                   <  1:3252101/150‑1 (MQ=255)
                                   cTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGa                                                                                                                <  1:3258295/150‑1 (MQ=255)
                                   cTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGa                                                                                                                <  1:3258133/150‑1 (MQ=255)
                                    tGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAg                                                                                                               <  2:3720044/150‑1 (MQ=255)
                                    tGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAg                                                                                                               <  2:2177149/150‑1 (MQ=255)
                                     ggACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAgg                                                                                                              <  1:2727229/150‑1 (MQ=255)
                                     ggACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAgg                                                                                                              <  2:2251500/150‑1 (MQ=255)
                                        cGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGc                                                                                                           <  2:3477593/150‑1 (MQ=255)
                                        cGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGc                                                                                                           <  1:3858812/150‑1 (MQ=255)
                                        cGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGc                                                                                                           <  1:1458127/150‑1 (MQ=255)
                                            ggATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCt                                                                                                       >  2:4040933/1‑150 (MQ=255)
                                               tCCGCGGAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcg                                                                                                    <  2:1440347/150‑1 (MQ=255)
                                                ccGCGTAACGCCTACCGTGAGGCCAGTTGTAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgc                                                                                                   <  1:4199231/150‑1 (MQ=255)
                                                ccGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgc                                                                                                   <  2:1028573/150‑1 (MQ=255)
                                                ccGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgc                                                                                                   <  2:1111457/150‑1 (MQ=255)
                                                ccGCGTAACCCCGACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgc                                                                                                   <  2:1111409/150‑1 (MQ=255)
                                                 cgcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgcg                                                                                                  <  1:3642381/150‑1 (MQ=255)
                                                 cgcgTAACCCCGCCCGTGAGGCCAGCTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgcg                                                                                                  <  2:3180029/150‑1 (MQ=255)
                                                   cgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGc                                                                                                <  2:1818914/150‑1 (MQ=255)
                                                   cgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGc                                                                                                <  1:3330046/150‑1 (MQ=255)
                                                    gTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGcc                                                                                               <  2:1209139/150‑1 (MQ=255)
                                                    gTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGcc                                                                                               <  1:1821755/150‑1 (MQ=255)
                                                     tAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGccc                                                                                              <  1:3563466/150‑1 (MQ=255)
                                                       aCCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGa                                                                                            <  1:3011199/150‑1 (MQ=255)
                                                       aCCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGa                                                                                            <  1:2728501/150‑1 (MQ=255)
                                                       aCCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGa                                                                                            <  1:1614071/150‑1 (MQ=255)
                                                           cTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATcgc                                                                                        <  1:4040860/150‑1 (MQ=255)
                                                           cTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATcgc                                                                                        <  1:3845178/150‑1 (MQ=255)
                                                           cTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATcgc                                                                                        <  2:3548397/150‑1 (MQ=255)
                                                                  gtggCCTGTTTATGGGGGTCACCGCTTTAACTGCTGCCGGGGCCGGTGTGGGAGCGGCTTTGTTTGGCGATCGGGCCGCACAGCGGCCCCAGCATTCTTTGCCGCTATACAGACATCGGGGGGCCGCTTCGCGGCCCGATCGCGACACaa                                                                                 <  2:3934306/148‑1 (MQ=255)
                                                                      ccttttGAGGGGGGTCACCGCCTTAACAGCTGCCGGGGCGCGTGGGGGCGGGGCCTTTTGTGGGGTTGGGGCCGCACTGGGGCCCCTGCTTTGTTTGCCGGTATACAGACCTCGTGGGGCCGCTGGGCGGCCCGATCGCGACACAAGGcc                                                                             <  2:3302351/146‑1 (MQ=255)
                                                                                                                          gCCTTGTGTCGCGCTAGGGCCGCCCCGCGGCCCCAGCAGTCTTTGCCGCGATACACACATCGAGGGGCCGCGTCGCGGCCCGAGAGCGACACAAGGCCGCTCCCACCAAAAGAACACCTCCACCTCGACAGCCGGAACACCGCTCCAGcc                         >  2:3271296/1‑150 (MQ=255)
                                                                                                                             ttGTGTCGCGATCGGGCCGCACAGCGGCCCCCGCAGTCATTTCCGCGCGACAGAAAAAGAGGGGGCGCGTCGCGGCCCGATCGCGCCACAAGGCCGCGCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGa                      >  2:2281690/1‑150 (MQ=255)
                                                                                                                                  tCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATcccc                 >  2:999976/1‑150 (MQ=255)
                                                                                                                                    gcgATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCg               >  1:1241781/1‑150 (MQ=255)
                                                                                                                                     tgctggggCCGCTGTGCGGCCCGAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGc              >  1:3329030/6‑150 (MQ=39)
                                                                                                                                     tgctggggCCGCTGTGCGGCCCGAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGc              >  1:1040591/6‑150 (MQ=39)
                                                                                                                                      gctggggCCGCTGTGCGGCCCGAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACACGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCa             >  2:3780395/5‑150 (MQ=255)
                                                                                                                                      gATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCa             >  1:3089445/1‑150 (MQ=255)
                                                                                                                                      gATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCa             >  1:3036416/1‑150 (MQ=255)
                                                                                                                                          gggCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACCCAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGAt         >  1:3334376/1‑150 (MQ=255)
                                                                                                                                          gggCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGAt         >  2:1966560/1‑150 (MQ=255)
                                                                                                                                           ggCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATc        >  1:997157/1‑150 (MQ=255)
                                                                                                                                           ggCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATc        >  2:4212608/1‑150 (MQ=255)
                                                                                                                                             ccGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCAc      >  1:4058571/1‑150 (MQ=255)
                                                                                                                                             ccGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCAc      >  2:2796564/1‑150 (MQ=255)
                                                                                                                                             ccGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCAc      >  2:4059495/1‑150 (MQ=255)
                                                                                                                                             ccGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCAc      >  2:4134100/1‑150 (MQ=255)
                                                                                                                                               gCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGc    >  1:77150/1‑150 (MQ=255)
                                                                                                                                               gCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGc    >  1:77094/1‑150 (MQ=255)
                                                                                                                                               gCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGc    >  2:3348989/1‑150 (MQ=255)
                                                                                                                                                 acaGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTg  >  1:3583988/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTG  >  NC_002947/4553795‑4554089

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