Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,701,591 C→T 8.3% intergenic (+38/‑174) pyrC' → / → hslV dihydroorotase‑like protein/ATP‑dependent HslVU protease peptidase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,701,5910CT 8.3% 397.1 / 18.8 132intergenic (+38/‑174)pyrC'/hslVdihydroorotase‑like protein/ATP‑dependent HslVU protease peptidase subunit
Reads supporting (aligned to +/- strand):  ref base C (38/83);  new base T (4/7);  total (42/90)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 140 total aligned reads displayed.

TGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCAC  >  NC_002947/5701442‑5701735
                                                                                                                                                     |                                                                                                                                                 
tGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCGTAGCGGGCACGCCCGCTCCCCCAGGATc                                                                                                                                                   >  2:1613853/1‑150 (MQ=255)
  cTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATcac                                                                                                                                                 >  2:4103702/1‑150 (MQ=255)
  cTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATcac                                                                                                                                                 <  1:1582625/150‑1 (MQ=255)
   tGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACCCGATcacc                                                                                                                                                >  1:2940278/1‑150 (MQ=255)
    ggCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGGTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCAGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATcacca                                                                                                                                               <  2:5621075/150‑1 (MQ=255)
     gCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATcacccc                                                                                                                                              >  1:4585726/1‑148 (MQ=255)
     gCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATcaccac                                                                                                                                              <  2:3540123/150‑1 (MQ=255)
      cGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACcaca                                                                                                                                             >  2:157221/1‑150 (MQ=255)
         aCAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACaagt                                                                                                                                          >  1:1053813/1‑148 (MQ=255)
           aaTGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCCCAATTcc                                                                                                                                        >  2:5057723/1‑150 (MQ=255)
           aaTGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTcc                                                                                                                                        <  1:4439316/150‑1 (MQ=255)
           aaTGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTcc                                                                                                                                        <  2:3085308/150‑1 (MQ=255)
            aTGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCa                                                                                                                                       <  2:4091881/150‑1 (MQ=255)
            aTGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCa                                                                                                                                       <  1:3617199/150‑1 (MQ=255)
             tGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCCCAGGCTCACCCCAATTCCaa                                                                                                                                      >  1:2828271/1‑150 (MQ=255)
             tGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCaa                                                                                                                                      <  1:1644079/150‑1 (MQ=255)
             tGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCaa                                                                                                                                      <  1:1443139/150‑1 (MQ=255)
             tGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCaa                                                                                                                                      <  1:4918681/150‑1 (MQ=255)
             tGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCaa                                                                                                                                      >  2:2265736/1‑150 (MQ=255)
               gTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAgg                                                                                                                                    <  2:1331246/150‑1 (MQ=255)
                 tctcCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCa                                                                                                                                  >  1:699495/1‑150 (MQ=255)
                  ctcCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCAt                                                                                                                                <  2:1197460/150‑1 (MQ=255)
                    cccGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCAtgt                                                                                                                              <  1:1494247/150‑1 (MQ=255)
                     ccGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACCATTCCAAGGCAtgtg                                                                                                                             >  2:3382966/1‑150 (MQ=255)
                     ccGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCAtgtg                                                                                                                             <  1:2255376/150‑1 (MQ=255)
                     ccGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGGCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACACTTCCAAGGCAtgtg                                                                                                                             >  2:2531596/1‑150 (MQ=255)
                         ggCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTaa                                                                                                                         <  2:4801688/150‑1 (MQ=255)
                         ggCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTaa                                                                                                                         <  1:3765049/150‑1 (MQ=255)
                          gCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCCTGCTTAGCGGGCACGCCCGCTCCCACAGGATCCCCACAATTCCCAGGCCGGTGGTAAc                                                                                                                        >  1:2388890/1‑150 (MQ=255)
                          gCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAAc                                                                                                                        <  2:699495/150‑1 (MQ=255)
                             agaACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCt                                                                                                                     <  1:1198468/150‑1 (MQ=255)
                               aaCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgt                                                                                                                   <  1:1066890/150‑1 (MQ=255)
                               aaCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgt                                                                                                                   <  1:1066839/150‑1 (MQ=255)
                                aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgtg                                                                                                                  <  2:1011586/150‑1 (MQ=255)
                                aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgtg                                                                                                                  <  2:972204/150‑1 (MQ=255)
                                aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgtg                                                                                                                  <  2:386768/150‑1 (MQ=255)
                                 cTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTgg                                                                                                                 <  1:293556/150‑1 (MQ=255)
                                  tGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTggc                                                                                                                <  2:65026/150‑2 (MQ=255)
                                    cccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGAg                                                                                                              <  2:3537171/150‑1 (MQ=21)
                                    cccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGAg                                                                                                              <  2:3537199/150‑1 (MQ=21)
                                    cccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAg                                                                                                              >  2:4564554/1‑150 (MQ=255)
                                     ccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGc                                                                                                             <  2:3096767/150‑1 (MQ=255)
                                       gTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCgg                                                                                                           <  1:4396560/150‑1 (MQ=255)
                                        ttCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCggg                                                                                                          <  1:3850756/150‑1 (MQ=255)
                                         tCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGc                                                                                                         <  2:1000252/150‑1 (MQ=255)
                                           aTCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCAt                                                                                                       <  2:5253364/150‑1 (MQ=255)
                                           aTCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCAt                                                                                                       <  1:3306268/150‑1 (MQ=255)
                                           aTCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAAATCCAAGGCATGTGGTAACCCAGTGGGAGCGGGcct                                                                                                       >  1:3917348/1‑148 (MQ=255)
                                            tCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGAGCgatgact                                                                                                      >  1:1578427/1‑144 (MQ=21)
                                             cGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAAGTCCAAGGCATGTGGTAACCCAGTGGGAGCGGGccggc                                                                                                     >  1:2966236/1‑146 (MQ=255)
                                                gcACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCTACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCcg                                                                                                  >  2:5622587/2‑150 (MQ=255)
                                                ccACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAAATCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCAGGCCcg                                                                                                  >  2:757892/1‑150 (MQ=255)
                                                   ctgcctgcCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCcgctc                                                                                               <  1:3399542/150‑3 (MQ=255)
                                                    tgcctgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGaa                                                                                              <  1:2628999/150‑1 (MQ=255)
                                                    tgcctgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGaa                                                                                              <  1:2629010/150‑1 (MQ=255)
                                                    tgcctgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGaa                                                                                              <  1:4474012/150‑1 (MQ=255)
                                                        tgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                                               <  1:2847083/129‑1 (MQ=18)
                                                        tgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                                               >  2:2847083/1‑129 (MQ=18)
                                                         gccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                                               <  1:2365596/128‑1 (MQ=18)
                                                         gccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                                               <  1:2365386/128‑1 (MQ=18)
                                                         gccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                                               >  2:2365386/1‑128 (MQ=18)
                                                           cGGGCGCGGTGCGTTATCCCTGGGTCGCGGGCCACCTCTGCCACGAGGCCTGAGTACCATGCTGAGCGGGCACGCCCGCTCCCCCAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGGGGGAGCGGGCATGCCCGCGAATGCCTTc                                                                                       <  2:4548037/150‑1 (MQ=255)
                                                           cGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCCGGCCCGCGACTGGCTTc                                                                                       >  2:2410211/1‑150 (MQ=255)
                                                            gggCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACCATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCCTGCCCGCGCAGGCCGTcc                                                                                      >  1:1250347/1‑149 (MQ=255)
                                                            gagcgcGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACCGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTca                                                                                      <  2:958142/148‑1 (MQ=255)
                                                               cgcgGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACAt                                                                                   <  2:5378931/150‑1 (MQ=255)
                                                               cgcgGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACAt                                                                                   <  2:5378991/150‑1 (MQ=255)
                                                                gcgGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATa                                                                                  <  1:2612426/150‑1 (MQ=255)
                                                                  ggTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCCTGCCCGCGAATGGCTGCACCTAgg                                                                                >  2:3668215/1‑149 (MQ=255)
                                                                      cGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACACTTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCCTGCCCGCGAATGCCGGCGCATAgggcgc                                                                            >  2:2486333/1‑145 (MQ=38)
                                                                        ttATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCCCAGGATCACCACCATTCCAAGGCCTGTGGTACCCCTGTGGGAGCGGGCCGGCCCGCGAATGCCTGCCCCGACTcccccc                                                                          >  2:3775920/1‑150 (MQ=255)
                                                                        ttATACCGTGGTCGAGGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTcccccc                                                                          <  1:5343198/150‑1 (MQ=255)
                                                                         tataCCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCCTGCCCGCGACTGCCGTCACCAAGTCccgccg                                                                         >  2:4594188/1‑150 (MQ=39)
                                                                         tataCCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGCAGGCCTTCACCTAGtgcggcgg                                                                         >  2:1759573/1‑143 (MQ=39)
                                                                            aCCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACACTTCCCAGGCCTGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGACTGCCGTCCCCTCGGCCCCCCggcc                                                                      >  2:5029094/1‑147 (MQ=37)
                                                                            aCCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCCAGGCATGTGGTAACCCAGTGGGAGCGGGCATGCCCGCGAATGCCGTCACAGAGGCGCCCCGccg                                                                      >  1:4489266/1‑149 (MQ=38)
                                                                             ccTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCaa                                                                     <  1:611445/150‑1 (MQ=255)
                                                                                          ccACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCa                                                        <  1:4491210/150‑1 (MQ=255)
                                                                                              aTCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATccc                                                    <  1:1451265/150‑1 (MQ=255)
                                                                                                     aCGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAAtt                                                                                                                                          >  1:1323729/1‑58 (MQ=255)
                                                                                                     aCGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAAtt                                                                                                                                          <  2:1323729/58‑1 (MQ=255)
                                                                                                     aCGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGt                                             <  1:1415337/150‑1 (MQ=255)
                                                                                                     aCGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGt                                             <  1:4253805/150‑1 (MQ=255)
                                                                                                         ggCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCCCAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGaa                                         <  2:1636510/150‑1 (MQ=255)
                                                                                                               aGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCtt                                   <  2:3465623/150‑1 (MQ=255)
                                                                                                               aGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCtt                                   <  2:3465497/150‑1 (MQ=255)
                                                                                                                 tACCATGCTTAGCGGGCACGCCCGCTCCCCCAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTcc                                 <  1:2410211/150‑1 (MQ=255)
                                                                                                                   ccATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTcccc                               <  1:1859235/150‑1 (MQ=255)
                                                                                                                      tGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCccccca                            <  1:769747/150‑1 (MQ=255)
                                                                                                                           aGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccc                       <  1:5325347/150‑1 (MQ=255)
                                                                                                                           aGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccc                       <  1:5325295/150‑1 (MQ=255)
                                                                                                                            gCGGGCACGCCCGCTCCCACAGGATCACCCCAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGACTGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccca                      <  2:1964420/150‑1 (MQ=255)
                                                                                                                            gCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccca                      <  1:4436310/150‑1 (MQ=255)
                                                                                                                            gCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccca                      <  2:1964754/150‑1 (MQ=255)
                                                                                                                              cggCACGCCCGCTCCCCCAGGAACACCACACGTCCACGGCATGTGGCAACCCGGTGGGAGCGGGCATGCCCGCGAATGCCTTCACAGTGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCata                    <  2:2828271/149‑1 (MQ=255)
                                                                                                                                 cACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCTTGTGGTAACCCTGTGGGATCGGGCATGCCCGCGTATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCACCGGTTGAATTCCTTGCCCCCACCCCCATATgc                 <  2:4830569/150‑2 (MQ=255)
                                                                                                                                  aCGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTca                                                                                      <  1:4285363/80‑1 (MQ=255)
                                                                                                                                   cGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCAGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGt               <  1:2265736/150‑1 (MQ=255)
                                                                                                                                      ccGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTg            <  1:3334444/150‑1 (MQ=255)
                                                                                                                                                cggATCACCACCCGTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGCGTCTGCATCAGGCAc  <  2:1224110/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TGCTGGCGAACAATGGTTCTCCCGTGGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCAC  >  NC_002947/5701442‑5701735

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