Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,553,940 A→T 9.0% 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,553,9400AT 9.0% 269.1 / 12.3 86intergenic (+68/+113)pydA/rutRNADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family
Reads supporting (aligned to +/- strand):  ref base A (26/52);  new base T (5/3);  total (31/55)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCT  >  NC_002947/4553793‑4554088
                                                                                                                                                   |                                                                                                                                                    
aaaTGGAGGCACAGGATACCGGCAAGCCGGACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGcaca                                                                                                                                                    <  2:5008954/150‑1 (MQ=255)
 aaTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAg                                                                                                                                                   <  2:4546797/150‑1 (MQ=255)
 aaTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAg                                                                                                                                                   >  1:346354/1‑150 (MQ=255)
    ggAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGACCCCTGTGCGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCgg                                                                                                                                                >  2:4022284/1‑150 (MQ=255)
       ggCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGccc                                                                                                                                             <  1:833489/150‑1 (MQ=255)
       ggCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGccc                                                                                                                                             <  2:4687475/150‑1 (MQ=255)
       ggCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGccc                                                                                                                                             <  2:4687520/150‑1 (MQ=255)
        tgtctgacatggccttcgagaatgcggtgcagacggtgaccagcgtgatcttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTTCGCGgcccg                                                                                                                                            >  2:5374629‑M1/107‑149 (MQ=255)
        gCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGcccc                                                                                                                                            <  1:2775250/150‑1 (MQ=255)
         cacaGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCca                                                                                                                                           <  2:3592693/150‑1 (MQ=255)
          acaGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcag                                                                                                                                          <  1:1288952/150‑1 (MQ=255)
           ctgacatggccttcgagaatgcggtgcagacggtgaccagcgtgatcttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTTCgcggcccgatc                                                                                                                                         >  1:3092925‑M1/104‑146 (MQ=255)
           ctgacatggccttcgagaatgcggtgcagacggtgaccagcgtgatcttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTTCgcggcccgatc                                                                                                                                         >  1:3592218‑M1/104‑146 (MQ=255)
           ctgacatggccttcgagaatgcggtgcagacggtgaccagcgtgatcttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTTCgcggcccgatc                                                                                                                                         >  2:3554057‑M1/104‑146 (MQ=255)
              gATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGt                                                                                                                                      <  1:3733277/150‑1 (MQ=255)
                tACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCt                                                                                                                                    <  2:3500584/150‑1 (MQ=255)
                   cGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTg                                                                                                                                 <  2:2294308/150‑1 (MQ=255)
                    ggCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGc                                                                                                                                <  1:3102313/150‑1 (MQ=255)
                         gCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGa                                                                                                                           >  1:1069427/1‑150 (MQ=255)
                          ccGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGAt                                                                                                                          <  1:3456449/150‑1 (MQ=255)
                          ccGTACCTGAACTGGACGCAGGATCCGCGTAACCCCGACCGTGAGGCCAGTTGAAGCGGGTCACCACCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGAt                                                                                                                          <  1:993765/150‑1 (MQ=255)
                             tACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACa                                                                                                                       <  2:3894488/150‑1 (MQ=255)
                             tACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACa                                                                                                                       <  1:3765464/150‑1 (MQ=255)
                             tACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACa                                                                                                                       <  2:2642804/150‑1 (MQ=255)
                             tACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACa                                                                                                                       <  2:2642729/150‑1 (MQ=255)
                                    aCTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGACCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCg                                                                                                                <  1:4022284/150‑1 (MQ=255)
                                       ggACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAgg                                                                                                             <  2:1465645/150‑1 (MQ=255)
                                        gACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAggg                                                                                                            <  2:4982978/150‑1 (MQ=255)
                                          ggtgaccagcgtgatcttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGc                                                                                                          >  2:3462034‑M1/73‑150 (MQ=255)
                                          cGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGc                                                                                                          <  2:1069427/150‑1 (MQ=255)
                                           gCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGcc                                                                                                         <  2:1727961/150‑1 (MQ=255)
                                             aGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGc                                                                                                       <  1:3532857/150‑1 (MQ=255)
                                              ggATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGa                                                                                                               >  2:2531481/1‑141 (MQ=255)
                                              ggATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGa                                                                                                               <  1:2531481/141‑1 (MQ=255)
                                               gATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCtt                                                                                                     <  1:3625085/150‑1 (MQ=255)
                                               gATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCtt                                                                                                     <  1:2096255/150‑1 (MQ=255)
                                                aTCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGcgaa                                                                                                    <  2:1635875/150‑4 (MQ=255)
                                                  gcgtgatcttgcgggggatcgggctggagccgtgttccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgc                                                                                                  >  1:2504935‑M1/65‑150 (MQ=255)
                                                  gcgtgatcttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTTCgcggcccgatc                                                                                                                                         <  2:1537972‑M1/47‑5 (MQ=255)
                                                  gcgtgatcttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCTTCgcggcccgatc                                                                                                                                         >  1:1537972‑M1/65‑107 (MQ=255)
                                                   cgtgatcttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTcgcg                                                                                                 >  1:1109144‑M1/64‑150 (MQ=255)
                                                   cgcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcg                                                                                                 >  2:3391248‑M2/1‑147 (MQ=255)
                                                   cgcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcg                                                                                                 >  2:346354‑M2/1‑147 (MQ=255)
                                                    gcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcgg                                                                                                >  2:4986638‑M2/1‑146 (MQ=255)
                                                    gcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcgg                                                                                                >  1:767965‑M2/1‑146 (MQ=255)
                                                      gTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcc                                                                                              >  1:5236537‑M2/1‑144 (MQ=255)
                                                      gTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcc                                                                                              >  2:2094882‑M2/1‑144 (MQ=255)
                                                         cttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGa                                                                                           >  1:3222274‑M1/58‑150 (MQ=255)
                                                          ttgcgggggatcgggctggagccgtgatccggcattccaggtggacgtggtcttttTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAGcggcccctc                                                                                          >  1:3986652‑M1/57‑147 (MQ=255)
                                                          ccccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAgcggcccctc                                                                                          >  1:2502783‑M2/1‑140 (MQ=255)
                                                          ccccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAgcggcccctc                                                                                          >  1:2502787‑M2/1‑140 (MQ=255)
                                                           cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatc                                                                                         >  2:4858568‑M2/1‑139 (MQ=255)
                                                           cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatc                                                                                         >  1:720746‑M2/1‑139 (MQ=255)
                                                               aCCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcga                                                                                     >  1:3537332‑M2/1‑135 (MQ=255)
                                                                       gCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCg                                                                                                                >  2:4276898/1‑115 (MQ=255)
                                                                       gCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCg                                                                                                                <  1:4276898/115‑1 (MQ=255)
                                                                                 gCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcg                                                                                                 <  2:4035707‑M2/120‑4 (MQ=255)
                                                                                 gCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcg                                                                                                 >  1:4035707‑M2/1‑117 (MQ=255)
                                                                                        aCCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGa                                                                                                                           <  1:4977586/87‑1 (MQ=255)
                                                                                        aCCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGa                                                                                                                           >  2:4977586/1‑87 (MQ=255)
                                                                                           gCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcga                                                                                     >  1:4492935‑M2/1‑107 (MQ=255)
                                                                                           gCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACAACGAGGGGCCGCTTCgcggcccgatcgcga                                                                                     <  2:4492935‑M2/122‑16 (MQ=255)
                                                                                                     cTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggccc                                                                                             <  1:66952‑M2/104‑8 (MQ=255)
                                                                                                     cTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggccc                                                                                             >  2:66952‑M2/1‑97 (MQ=255)
                                                                                                                    tgGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAgcggcccctc                                                                                          >  1:3637275‑M2/1‑82 (MQ=255)
                                                                                                                    tgGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCGAAgcggcccctc                                                                                          <  2:3637275‑M2/92‑11 (MQ=255)
                                                                                                                                    TCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATcccc                <  2:1295271‑M1/150‑1 (MQ=255)
                                                                                                                                    TCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATcccc                <  1:4354438‑M1/150‑1 (MQ=255)
                                                                                                                                     actgctgGGGCCGCTGTGCGGCCCGAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATccccc               <  1:1252848‑M1/143‑1 (MQ=255)
                                                                                                                                       tgctgGGGCCGCTGTGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGc             <  1:2791670‑M1/145‑1 (MQ=255)
                                                                                                                                       tgcTCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGc             <  1:262103‑M1/147‑1 (MQ=255)
                                                                                                                                       tgcTCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGc             <  1:1135857‑M1/147‑1 (MQ=255)
                                                                                                                                        gATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggg            <  2:721187‑M2/150‑89 (MQ=255)
                                                                                                                                        gATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggg            <  1:4658833‑M2/150‑89 (MQ=255)
                                                                                                                                        gATCGGGCCGCACAGCGGCCCCAGCAGTCGTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacgggcggg            <  2:3648997‑M2/150‑89 (MQ=255)
                                                                                                                                        GATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCa            <  1:5216457‑M1/150‑1 (MQ=255)
                                                                                                                                         aTCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtagggg           <  2:2414742‑M2/150‑90 (MQ=255)
                                                                                                                                            GGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGAt        <  2:695289‑M1/150‑1 (MQ=255)
                                                                                                                                             GGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATc       <  1:1796327‑M1/150‑1 (MQ=255)
                                                                                                                                             GGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATc       <  1:1163057‑M1/150‑1 (MQ=255)
                                                                                                                                              GCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCa      <  2:53784‑M1/150‑1 (MQ=255)
                                                                                                                                               CCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCAc     <  1:3927719‑M1/150‑1 (MQ=255)
                                                                                                                                                cGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccgcaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggttacgcg    <  2:1803120‑M2/150‑97 (MQ=255)
                                                                                                                                                cGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggttacgcg    <  2:1825420‑M2/150‑97 (MQ=255)
                                                                                                                                                 gCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggttacgcgg   <  2:60657‑M2/150‑98 (MQ=255)
                                                                                                                                                  cacaGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggc                        >  2:3135414‑M2/1‑52 (MQ=255)
                                                                                                                                                  cacaGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggc                        <  1:3135414‑M2/128‑77 (MQ=255)
                                                                                                                                                  CACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCt  <  1:877598‑M1/150‑1 (MQ=255)
                                                                                                                                                  CACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCt  <  1:877172‑M1/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCT  >  NC_002947/4553793‑4554088

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: