Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,701,615 A→G 9.8% intergenic (+62/‑150) 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,6150AG 9.8% 326.9 / 19.3 112intergenic (+62/‑150)pyrC'/hslVdihydroorotase‑like protein/ATP‑dependent HslVU protease peptidase subunit
Reads supporting (aligned to +/- strand):  ref base A (23/78);  new base G (4/7);  total (27/85)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

Only 100 of 114 total aligned reads displayed.

GGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCACTTTCGCCCCGCTGCGTGGAGACTCTCCC  >  NC_002947/5701467‑5701763
                                                                                                                                                     |                                                                                                                                                    
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)
   gagaACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAAccc                                                                                                                                                  >  1:3096767/1‑150 (MQ=255)
    agaACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCt                                                                                                                                                 <  1:1119740/150‑1 (MQ=255)
    agaACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCt                                                                                                                                                 <  1:1198468/150‑1 (MQ=255)
    agaACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCt                                                                                                                                                 <  2:2710551/150‑1 (MQ=255)
      aaCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgt                                                                                                                                               <  1:1066890/150‑1 (MQ=255)
      aaCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgt                                                                                                                                               <  1:1066839/150‑1 (MQ=255)
       aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgtg                                                                                                                                              <  2:972204/150‑1 (MQ=255)
       aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgtg                                                                                                                                              <  2:386768/150‑1 (MQ=255)
       aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgtg                                                                                                                                              <  2:1011586/150‑1 (MQ=255)
       aCTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCtgtg                                                                                                                                              <  1:2263639/150‑1 (MQ=255)
        cTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTgg                                                                                                                                             <  1:293556/150‑1 (MQ=255)
         tGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTggg                                                                                                                                            >  2:4916036/1‑150 (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)
            ccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCCCAAATCCAAGGCATGTGGTAACCCGGTGGgcgc                                                                                                                                         >  1:4414036/1‑150 (MQ=255)
            ccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGc                                                                                                                                         <  2:3096767/150‑1 (MQ=255)
            ccGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGAGGGCCACCTCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGc                                                                                                                                         <  2:122828/150‑1 (MQ=255)
              gTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCgg                                                                                                                                       <  1:4396560/150‑1 (MQ=255)
              gTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCgg                                                                                                                                       <  1:2928664/150‑1 (MQ=255)
              gTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCgg                                                                                                                                       <  2:101902/150‑1 (MQ=255)
              gTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTAGACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCgg                                                                                                                                       <  1:3927121/150‑1 (MQ=255)
               ttCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCggg                                                                                                                                      <  1:3850756/150‑1 (MQ=255)
                tCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGGTCAACACACTTCCAAGGCAAGTGGTAACCCTGTGGGAGCGGGc                                                                                                                                     >  1:1576328/1‑150 (MQ=255)
                tCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGc                                                                                                                                     <  2:1000252/150‑1 (MQ=255)
                  aTCGGCCTCTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCGTCGAGGGCTGAGTACCATGCTTCGCGGGCATGCCCGCTCCCCCAGGGTTACCACATGCCCTGGAATTGTGGTGATCCTGTGGGAGCGGGcgt                                                                                                                                   <  2:3400827/150‑3 (MQ=12)
                  aTCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCAt                                                                                                                                   <  1:4572357/150‑1 (MQ=255)
                  aTCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCAt                                                                                                                                   <  1:3306268/150‑1 (MQ=255)
                  aTCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCAt                                                                                                                                   <  2:5253364/150‑1 (MQ=255)
                   tCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGAGCgatgact                                                                                                                                  >  1:1578427/1‑144 (MQ=21)
                       gcACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCTACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCcg                                                                                                                              >  2:5622587/2‑150 (MQ=255)
                       ccACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAAATCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCAGGCCcg                                                                                                                              >  2:757892/1‑150 (MQ=255)
                          ctgcctgcCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCcgctc                                                                                                                           <  1:3399542/150‑3 (MQ=255)
                           tgcctgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGaa                                                                                                                          <  1:2629010/150‑1 (MQ=255)
                           tgcctgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGaa                                                                                                                          <  1:4474012/150‑1 (MQ=255)
                            gcctgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAAt                                                                                                                         <  2:3843089/150‑1 (MQ=255)
                               tgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                                                                           >  2:2847083/1‑129 (MQ=18)
                               tgccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCATGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                                                                           <  1:2847083/129‑1 (MQ=18)
                                gccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                                                                           >  2:2365596/1‑128 (MQ=18)
                                gccGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGGTTACCACATGCCTTGGAATTGTGGTGATCCTGTGGGa                                                                                                                                           <  1:2365596/128‑1 (MQ=18)
                                  cGGGCGCGGTGCGTTATCCCTGGGTCGCGGGCCACCTCTGCCACGAGGCCTGAGTACCATGCTGAGCGGGCACGCCCGCTCCCCCAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGGGGGAGCGGGCATGCCCGCGAATGCCTTc                                                                                                                   <  2:4548037/150‑1 (MQ=255)
                                  cGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCCGGCCCGCGACTGGCTTc                                                                                                                   >  2:2410211/1‑150 (MQ=255)
                                   gggCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACCATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCCTGCCCGCGCAGGCCGTcc                                                                                                                  >  1:1250347/1‑149 (MQ=255)
                                   gggCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCACGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTATCCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTca                                                                                                                  <  2:1597294/150‑1 (MQ=255)
                                   gagcgcGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACCGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTca                                                                                                                  <  2:958142/148‑1 (MQ=255)
                                      cgcgGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACAt                                                                                                               <  2:5378931/150‑1 (MQ=255)
                                       gcgGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATa                                                                                                              <  1:2612426/150‑1 (MQ=255)
                                       gcgGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATa                                                                                                              <  1:245271/150‑1 (MQ=255)
                                         ggTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACCATTCCCAGGCCTGTGGTAACCCTGTGGGAGCGGGCCTGCCCGCGCCTGCCGTcgcctcgg                                                                                                            >  2:838794/1‑143 (MQ=38)
                                         ggTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCCTGCCCGCGAATGGCTGCACCTAgg                                                                                                            >  2:3668215/1‑149 (MQ=255)
                                            gCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTccc                                                                                                         <  1:2385730/150‑1 (MQ=255)
                                             cGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACACTTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCCTGCCCGCGAATGCCGGCGCATAgggcgc                                                                                                        >  2:2486333/1‑145 (MQ=38)
                                               ttATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCCCAGGATCACCACCATTCCAAGGCCTGTGGTACCCCTGTGGGAGCGGGCCGGCCCGCGAATGCCTGCCCCGACTcccccc                                                                                                      >  2:3775920/1‑150 (MQ=255)
                                               ttATACCGTGGTCGAGGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTcccccc                                                                                                      <  1:5343198/150‑1 (MQ=255)
                                                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)
                                                        ggTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAgcgc                                                                                             <  2:3535762/150‑1 (MQ=255)
                                                                 ccACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCa                                                                                    <  1:4491210/150‑1 (MQ=255)
                                                                     aTCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATccc                                                                                <  1:1451215/150‑1 (MQ=255)
                                                                     aTCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATccc                                                                                <  1:1451265/150‑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:3465497/150‑1 (MQ=255)
                                                                                      aGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCtt                                                               <  2:3465623/150‑1 (MQ=255)
                                                                                       gTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTc                                                              <  1:745167/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:5325295/150‑1 (MQ=255)
                                                                                                  aGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccc                                                   <  1:5325347/150‑1 (MQ=255)
                                                                                                   gCGGGCACGCCCGCTCCCACAGGATCACCCCAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGACTGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccca                                                  <  2:1964420/150‑1 (MQ=255)
                                                                                                   gCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCggg                                                                                                                                      <  1:1725226/66‑1 (MQ=255)
                                                                                                   gCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCggg                                                                                                                                      >  2:1725226/1‑66 (MQ=255)
                                                                                                   gCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccca                                                  <  1:4436310/150‑1 (MQ=255)
                                                                                                   gCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCcccccaccccca                                                  <  2:1964754/150‑1 (MQ=255)
                                                                                                     gggCACGCCCGCGCCCCCCGGCTCCCCCCACTGCCCCGGGAGGTGGGAACCCGGTGGGAGCGGGCAGGCCCGCGAAGGCCGTCACCTAGGCCCCCCGCCAAGCGCGACCATTCAGCCCCCCGGGTTGAAATCCTTCCCCCCACCCCCata                                                <  2:5223480/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)
                                                                                                         aCGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTca                                                                                                                  >  2:4285363/1‑80 (MQ=255)
                                                                                                          cGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCAGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGt                                           <  1:2265736/150‑1 (MQ=255)
                                                                                                             ccGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTg                                        <  1:3334444/150‑1 (MQ=255)
                                                                                                                               cacaATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCACTTTCGccc                      <  1:2933623/150‑1 (MQ=255)
                                                                                                                               cacaATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCACTTTCGccc                      <  1:583582/150‑1 (MQ=255)
                                                                                                                               cacaATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCAGATGAGTCTGCATCAGGCACTTTCGccc                      <  1:564951/150‑1 (MQ=255)
                                                                                                                                          ggCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCACTTTCGCCCCGCTGCGTGga           <  1:2531596/150‑1 (MQ=255)
                                                                                                                                           gaaTGCGGTATCCCAGGGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCACTTTCGCCCCGCTGCGTGgag          <  1:3668215/148‑1 (MQ=255)
                                                                                                                                           gCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCGTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCACTTTCGCCCCGCTGCGTGgag          <  1:4259385/150‑1 (MQ=255)
                                                                                                                                            cATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCACTTTCGCCCCGCTGCGTGgaga         <  1:2877950/150‑1 (MQ=255)
                                                                                                                                               tgtgGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCACTTTCGCCCCGCTGCGTGGAGAct       <  1:2962766/150‑1 (MQ=255)
                                                                                                                                                  ggTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCACTTTCGCCCCGCTGCGTGGAGACtctc    <  1:1885753/150‑1 (MQ=255)
                                                                                                                                                    tAACCCTGTGGGAGCGGGCATGCCCGCGAAGGCCTTCACATAGGCCCCCCGCCAAGCGCGACCAGCCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATCCGAGTCTGCATCAGGCACTGTCGCCCCGCTGCGTGGAGACTCGccc  <  2:2388890/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCGAGAACTGCCCGTTCATCGGCCACTGCCTGCCGGGCGCGGTGCGTTATACCTTGGTCGATGGCCACATCTGCCACGAGGCCTGAGTACCATGCTTAGCGGGCACGCCCGCTCCCACAGGATCACCACAATTCCAAGGCATGTGGTAACCCTGTGGGAGCGGGCATGCCCGCGAATGCCTTCACATAGTCCCCCCGCCAAGCGCGACCATTCATCCCCCCCGGTTGAAATCCTTCCCCCCACCCCCATATGAGTCTGCATCAGGCACTTTCGCCCCGCTGCGTGGAGACTCTCCC  >  NC_002947/5701467‑5701763

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: