Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,556,967 A→T 11.1% intergenic (+68/+113) 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,556,9670AT11.1% 427.4 / 37.2 144intergenic (+68/+113)pydA/rutRNADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family
Reads supporting (aligned to +/- strand):  ref base A (30/97);  new base T (9/7);  total (40/104)
Fisher's exact test for biased strand distribution p-value = 1.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 150 total aligned reads displayed.

AAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTG  >  NC_002947/4556820‑4557116
                                                                                                                                                   |                                                                                                                                                     
aaaTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGcaca                                                                                                                                                     <  2:4617125/150‑1 (MQ=255)
aaaTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGcaca                                                                                                                                                     <  1:1342886/150‑1 (MQ=255)
     gAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGc                                                                                                                                                <  1:7197541/150‑1 (MQ=255)
      aGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCGACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTGTGCGGcc                                                                                                                                               <  1:5381853/150‑1 (MQ=255)
        gCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGcccc                                                                                                                                             <  1:4939902/150‑1 (MQ=255)
        gCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGcccc                                                                                                                                             <  1:7410752/150‑1 (MQ=255)
         cacaGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTGTGCGGcccga                                                                                                                                            <  2:3812546/150‑3 (MQ=255)
         cacaGGAGACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGTACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTGTGCGGcccga                                                                                                                                            <  2:2725785/150‑3 (MQ=255)
           caGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGTGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcagc                                                                                                                                          >  2:3954663/1‑150 (MQ=255)
           caGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCTATCGGGCCGCACAGCGGCCCcagc                                                                                                                                          <  1:1533563/150‑1 (MQ=255)
            aGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTGTGCGGcccgatcg                                                                                                                                         <  2:1693706/150‑6 (MQ=255)
            aGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcagca                                                                                                                                         <  2:5738033/150‑1 (MQ=255)
               aTACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTc                                                                                                                                      <  1:2497728/150‑1 (MQ=255)
               aTACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTc                                                                                                                                      <  2:5988146/150‑1 (MQ=255)
               aTACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTc                                                                                                                                      <  1:4212482/150‑1 (MQ=255)
               aTACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTc                                                                                                                                      <  2:4619185/150‑1 (MQ=255)
                 aCCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCtt                                                                                                                                    <  1:368352/150‑1 (MQ=255)
                   cGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTg                                                                                                                                  >  1:5250296/1‑150 (MQ=255)
                   cGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTg                                                                                                                                  <  1:1509129/150‑1 (MQ=255)
                    ggCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGccc                                                                                                                                              <  2:4947249/137‑1 (MQ=255)
                      cAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGTGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCcg                                                                                                                               <  1:3954663/150‑1 (MQ=255)
                       aaGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCcgc                                                                                                                              <  2:6089068/150‑1 (MQ=255)
                         gCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGa                                                                                                                            <  2:5250296/150‑1 (MQ=255)
                         gCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGa                                                                                                                            <  2:1649510/150‑1 (MQ=255)
                         gCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGa                                                                                                                            <  1:5078929/150‑1 (MQ=255)
                         gCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGa                                                                                                                            <  1:1899468/150‑1 (MQ=255)
                          ccGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGAt                                                                                                                           <  2:4453475/150‑1 (MQ=255)
                             tACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACa                                                                                                                        <  2:5140172/150‑1 (MQ=255)
                             tACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACa                                                                                                                        <  2:5140189/150‑1 (MQ=255)
                             tACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACa                                                                                                                        <  1:4384804/150‑1 (MQ=255)
                               ccTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGa                                                                                                                      <  2:7427696/150‑1 (MQ=255)
                               ccTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGa                                                                                                                      <  2:7300478/150‑1 (MQ=255)
                               ccTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGa                                                                                                                      <  2:4468351/150‑1 (MQ=255)
                                  gAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACAt                                                                                                                   <  2:1610864/150‑1 (MQ=255)
                                  gAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACAt                                                                                                                   <  2:1610785/150‑1 (MQ=255)
                                  gAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACAt                                                                                                                   <  1:6865658/150‑1 (MQ=255)
                                      tGGACGCAGGATCCGCGTAACCCCTCCCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAg                                                                                                               <  2:6358835/150‑1 (MQ=255)
                                      tGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAg                                                                                                               >  1:3679294/1‑150 (MQ=255)
                                      tGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAg                                                                                                               <  2:4810865/150‑1 (MQ=255)
                                       ggACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAgg                                                                                                              <  2:122179/150‑1 (MQ=255)
                                        gACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTTAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACACCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAggg                                                                                                             <  1:6072982/150‑1 (MQ=255)
                                         aCGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAgggg                                                                                                            <  2:1265518/150‑1 (MQ=255)
                                          cGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGc                                                                                                           <  2:2436367/150‑1 (MQ=255)
                                            cAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCg                                                                                                         <  1:6552512/150‑1 (MQ=255)
                                             aGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCTGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGa                                                                                                                            >  1:5457726/1‑130 (MQ=255)
                                             aGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCTGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGa                                                                                                                            <  2:5457726/130‑1 (MQ=255)
                                                 tCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcg                                                                                                    <  2:4263339/150‑1 (MQ=255)
                                                 tCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcg                                                                                                    <  2:5358386/150‑1 (MQ=255)
                                                  ccGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgc                                                                                                   <  1:1246559/150‑1 (MQ=255)
                                                  ccGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgc                                                                                                   <  2:3679294/150‑1 (MQ=255)
                                                   cgcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgcg                                                                                                  <  1:4668078/150‑1 (MQ=255)
                                                   cgcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgcg                                                                                                  <  2:6296647/150‑1 (MQ=255)
                                                   cgcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgcg                                                                                                  <  2:7297204/150‑1 (MQ=255)
                                                    gcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCgg                                                                                                 <  1:1141003/150‑1 (MQ=255)
                                                     cgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGc                                                                                                <  2:2776902/150‑1 (MQ=255)
                                                     cgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGc                                                                                                <  2:2776986/150‑1 (MQ=255)
                                                     cgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGc                                                                                                <  2:6491067/150‑1 (MQ=255)
                                                      gTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAGCGGcc                                                                                               <  1:6193404/150‑1 (MQ=255)
                                                       tAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGccc                                                                                              <  1:4429677/150‑1 (MQ=255)
                                                       tAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGccc                                                                                              <  1:3140533/150‑1 (MQ=255)
                                                       tAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGccc                                                                                              <  1:1207610/150‑1 (MQ=255)
                                                       tAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGccc                                                                                              <  1:2252801/150‑1 (MQ=255)
                                                       tAACCCCGACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGccc                                                                                              <  1:5218379/150‑1 (MQ=255)
                                                       gacccccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGGACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAGCGGccc                                                                                              <  2:427030/147‑1 (MQ=38)
                                                       cgtccccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGccc                                                                                              <  2:5702587/147‑1 (MQ=255)
                                                         aCCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAGCGGcccct                                                                                            <  2:471242/150‑3 (MQ=255)
                                                           cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATc                                                                                          <  2:6822846/150‑1 (MQ=255)
                                                           cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATc                                                                                          <  1:1353342/150‑1 (MQ=255)
                                                           cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATc                                                                                          <  2:3154391/150‑1 (MQ=255)
                                                             cgaCCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATcgc                                                                                        <  1:7138349/148‑1 (MQ=255)
                                                             cTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATcgc                                                                                        <  2:369817/150‑1 (MQ=255)
                                                             cTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATcgc                                                                                        <  2:3055302/150‑1 (MQ=255)
                                                                ccGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGac                                                                                     <  2:2692305/150‑1 (MQ=255)
                                                                    gcggCCAGTTTAAGCGGGTCACCGCGTGAACAGCTGCCGCGGCCCGTGTGGGAGCGGCCTTGTGTCGCGATGGGGCCGCACAGCGGCCCCAGCTGTCTTTGCCGCGATACAGACATCGTGGGGCCGCTTCGCGGCCCGATCGCGACACaa                                                                                 <  2:6146631/148‑1 (MQ=37)
                                                                                       cACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGc                                                                                                                                                   <  1:4962088/65‑1 (MQ=255)
                                                                                       cACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGc                                                                                                                                                   >  2:4962088/1‑65 (MQ=255)
                                                                                                    gCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAggg                                                                                                             <  1:2569812/90‑1 (MQ=255)
                                                                                                                        aGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgc                                                                                                   >  1:6450152/1‑80 (MQ=255)
                                                                                                                        aGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgc                                                                                                   <  2:6450152/80‑1 (MQ=255)
                                                                                                                               ttGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCCGCCCCGTCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAAAGCCGGATCAAGGCTCCAGCCCGa                      >  2:3538144/1‑150 (MQ=38)
                                                                                                                                 gtgtCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGCTc                    >  2:6073334/1‑150 (MQ=255)
                                                                                                                                    tCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCAGGAATGCCGGATCACGGCTCCAGCCCGATcccc                 >  1:1002413/1‑150 (MQ=255)
                                                                                                                                      gcgATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCg               >  1:6243240/1‑150 (MQ=255)
                                                                                                                                       tgctggggCCGCTGTGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGc                                           <  2:6895737/116‑1 (MQ=38)
                                                                                                                                       tgctggggCCGCTGTGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGc              >  1:5147388/6‑150 (MQ=255)
                                                                                                                                         ctggggCCGCTGTGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCaa            >  2:396162/4‑150 (MQ=255)
                                                                                                                                          tggggCCGCTGTGCGGCCCGAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATcg                                                                                         <  1:5718255/70‑1 (MQ=38)
                                                                                                                                          tggggCCGCTGTGCGGCCCGAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATcg                                                                                         >  2:5718255/3‑72 (MQ=38)
                                                                                                                                          tggggCCGCTGTGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAg           >  2:3406060/3‑150 (MQ=255)
                                                                                                                                          tggggCCGCTGTGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAg           >  1:3191300/3‑150 (MQ=255)
                                                                                                                                             ggCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATc        >  1:6764376/1‑150 (MQ=255)
                                                                                                                                              gCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCa       >  2:768536/1‑150 (MQ=255)
                                                                                                                                               ccGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCAc      >  2:2717728/1‑150 (MQ=255)
                                                                                                                                                cGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACg     >  1:868287/1‑150 (MQ=255)
                                                                                                                                                cGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGACCGATCCCCCGCAAGATCACg     >  2:3316107/1‑150 (MQ=255)
                                                                                                                                                 gCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGc    >  2:2632544/1‑150 (MQ=255)
                                                                                                                                                 gCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGc    >  1:531854/1‑150 (MQ=255)
                                                                                                                                                 gCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGc    >  1:5587419/1‑150 (MQ=255)
                                                                                                                                                 gCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGc    >  1:5587427/1‑150 (MQ=255)
                                                                                                                                                   acaGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTg  >  2:2474432/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTG  >  NC_002947/4556820‑4557116

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