Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,673,046:1 +T 5.5% intergenic (+72/‑95) prpF → / → prpD aconitate isomerase/2‑methylcitrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,673,0461.T 5.5% 799.4 / 36.7 183intergenic (+72/‑95)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base . (74/99);  new base T (2/8);  total (76/107)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 185 total aligned reads displayed.

TAACCAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACA  >  NC_002947/2672898‑2673193
                                                                                                                                                     |                                                                                                                                                   
tAACCAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCGCCCCCACGTaa                                                                                                                                                    >  1:5593810/1‑150 (MQ=255)
  aCCAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGGTTGCCCCCCCCCACACGGAAta                                                                                                                                                  >  1:2173141/1‑149 (MQ=255)
    cAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATcac                                                                                                                                                >  2:4897078/1‑150 (MQ=255)
     aaGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATcaca                                                                                                                                               <  2:8470916/150‑1 (MQ=255)
      aGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAg                                                                                                                                              <  2:1371876/150‑1 (MQ=255)
      aGGAAATCATGAGCCGTTGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCTCTGGCGCTCGGCTTCAAGAGGGGCCGGTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAg                                                                                                                                              <  1:6941701/150‑1 (MQ=255)
       ggCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAgg                                                                                                                                             <  1:117384/150‑1 (MQ=255)
       ggCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAgg                                                                                                                                             <  1:5496490/150‑1 (MQ=255)
       ggCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAgg                                                                                                                                             <  2:2392687/150‑1 (MQ=255)
       ggCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAgg                                                                                                                                             <  1:6669420/150‑1 (MQ=255)
       ggCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAgg                                                                                                                                             <  2:5413132/150‑1 (MQ=255)
       ggCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTATCACAgca                                                                                                                                            >  2:3815700/1‑148 (MQ=255)
        gCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAggg                                                                                                                                            <  2:2530228/150‑1 (MQ=255)
            tCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCTCCCACAGGTAATCACAGGGTTTg                                                                                                                                        >  1:2105376/1‑150 (MQ=255)
             cATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGa                                                                                                                                       <  1:7636823/150‑1 (MQ=255)
               tGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGatc                                                                                                                                     >  1:6803614/1‑148 (MQ=255)
               tGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAAt                                                                                                                                     <  2:4760862/150‑1 (MQ=255)
                gAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATg                                                                                                                                    >  2:5308032/1‑150 (MQ=255)
                  gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                  <  2:4070351/150‑1 (MQ=255)
                    cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCTCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                <  2:2337260/150‑1 (MQ=255)
                    cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                <  2:371205/150‑1 (MQ=255)
                    cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                <  2:2134032/150‑1 (MQ=255)
                    cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                <  2:2098807/150‑1 (MQ=255)
                    cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGGTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGTTAATCACAGGGTTTGAATGCtgt                                                                                                                                <  2:5791545/150‑1 (MQ=255)
                     gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                               <  1:4375248/150‑1 (MQ=255)
                       aGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTTCCCCCGCCCACAGGTAATAACAGGGTTTGAATGCTGTGat                                                                                                                             >  1:4117228/1‑150 (MQ=255)
                            tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATccc                                                                                                                        <  1:6339694/150‑1 (MQ=255)
                             cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                       <  1:2522371/150‑1 (MQ=255)
                             cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCACCCACACGGAAACAAACCGGTTTAAAGCTGTGCTATcacg                                                                                                                       >  1:6614064/1‑147 (MQ=255)
                              gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                      <  2:8668073/150‑1 (MQ=255)
                              gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCACCCACACGGAATAACACGGTTTTAATGCTGTTattaaccgg                                                                                                                      >  1:8487092/1‑143 (MQ=255)
                               cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                     <  2:4411118/150‑1 (MQ=255)
                               cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                     <  2:6337549/150‑1 (MQ=255)
                                gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                   <  1:4916224/150‑1 (MQ=25)
                                gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCTCCCACAGGGAATAACAGGGTTTGAATTCTGTGATATACCggtg                                                                                                                    >  1:1385922/1‑150 (MQ=255)
                                 tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCTGCTGTGCAGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                   <  2:6346370/150‑1 (MQ=17)
                                   cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                 <  2:8031394/150‑1 (MQ=18)
                                    tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGGTTGAATGCTGTGATATCCCGGTGGGAg                                                                                                                >  2:1478553/1‑150 (MQ=255)
                                    tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGAATCACAGGGGTTGAATGCGGTGATATCCCTGTGGGAg                                                                                                                >  1:7240420/1‑150 (MQ=255)
                                     gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGc                                                                                                               <  2:4811268/150‑1 (MQ=18)
                                     gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTGCCCCCGCCCCCCCGGAAAAAAAGGGGTTGTATTGTGTGATATCCCGGTGGGGGc                                                                                                               >  1:5942817/1‑150 (MQ=255)
                                       tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                             <  1:3084381/150‑1 (MQ=255)
                                       tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                             <  2:7377570/150‑1 (MQ=255)
                                       tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                             <  2:7818233/150‑1 (MQ=18)
                                       tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                             <  2:1350570/150‑1 (MQ=18)
                                           gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAACGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCaa                                                                                                         <  1:8449119/150‑1 (MQ=17)
                                             gCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGc                                                                                                       <  1:6458850/150‑1 (MQ=255)
                                              cTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGcc                                                                                                      >  2:2437845/1‑150 (MQ=255)
                                                  gTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATTCTGTGATATACCGGTGGGAGCGGGCAAGCCcgcg                                                                                                  >  2:1061464/1‑150 (MQ=255)
                                                      gcgTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCTACCACAGGTAAAAACAGGGTTTGAAAGCGGTGATAAACCTGTGGGAGCGGGCAACCCCGCGAAgg                                                                                              >  2:3945009/1‑150 (MQ=255)
                                                        gTGCCTGGCGATCGCTACTAAAGCACCACGGGCGCTCGGCTCCAAGCGGGGCCGCTTTGCGGCCCTTCGCGGGCTGGCCCGCTCCCCCAGGTAATCCCTGGGTTTGACTGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGc                                                                                            <  2:3560846/150‑1 (MQ=37)
                                                          gCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGAAAGCCCGCGAAGGGCTg                                                                                          <  2:2040176/150‑1 (MQ=255)
                                                           ccTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGc                                                                                         >  2:8250300/1‑150 (MQ=255)
                                                               gCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                    >  1:2432908/1‑119 (MQ=255)
                                                                cGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATCCCAGGGTTTGAATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCGGCAcagc                                                                                    >  2:2198684/1‑150 (MQ=255)
                                                                    aGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCACAGGTAATCACAGGGGTTGTATGCTGTGATATCCCTGGGGGAGCGGGCAAGCCCGCGAAGGGCTGCACCGCAGcc                                                                                >  2:3370791/1‑150 (MQ=38)
                                                                      cTTCTAACGCACCACTGGCGCTCGGGTTCACGCGGGGCCGCTGTGCGGCCCTGCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGcccc                                                                              <  2:6953653/150‑1 (MQ=25)
                                                                      cTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCACAGGTAATAACAGGGTTTTTATGCTGTGATATACCGGTGGGAGCGGGCCAGCCCCCGAAGGGCTGCACAGCAGcccc                                                                              >  2:2397193/1‑150 (MQ=37)
                                                                          tAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCCCAGGTAATCACAGGGGTTGAATGCTGTGATATCCCGGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCa                                                                          >  2:1582688/1‑150 (MQ=255)
                                                                           aaaGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGGCCGCGCCCCCAGGTAATCACAGGGGTTGTATTGGGTGGTATACCAGTGGGTGCGGGCAACCCCCCCAAGGGCGGCCCAGCAGCCCCCGCaa                                                                         >  1:1870784/1‑150 (MQ=255)
                                                                                  cccTGGCGCTGGGATGCAAGGGGGGCGGATTTGCGGCCATTCGCGGGCTTGCCCGCTCCCACAGTTAATCACTGGTTTTGATTGGTGTGATATCCGTTTGGGAGGGTCCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                  <  1:8076591/148‑1 (MQ=255)
                                                                                  caggggCGCTCGGCTACACGGGGGGCGGGTTGGCGGCCCTGGGGGGGCTGGCCGGCGCCCACCGTTAATCACCGGGTTTTAAGGGTGTGATATCCGTGTGGTAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                  <  1:7119329/146‑1 (MQ=255)
                                                                                  cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGTATTGTGTGATATACCGGTGGGTGCGGGCAAGCCCGCGAATGGCTGCACAGCAGCCCCGGCAAttattaa                                                                  >  2:1907279/1‑145 (MQ=37)
                                                                                  cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                             >  2:5813098/1‑107 (MQ=255)
                                                                                  cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                             <  1:5812365/107‑1 (MQ=255)
                                                                                        cgcTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGc                                                            <  2:3541998/150‑1 (MQ=255)
                                                                                            cGGCTGCAAGAGGGGCCGGTTTGCGGCCCTTGGCGGGCTTGCCCGCTCCCACAGTTAATCACAGGGTTTGAATGCTGTTATATCCCTGTGGGAGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGa                                                        <  2:6614064/150‑1 (MQ=38)
                                                                                                ttCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATCACAGGGTTTGAATGCTGTGATATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGGGCTGGATATg                                                    >  2:6215556/1‑150 (MQ=255)
                                                                                                 tctagcggggCGGCTTTGCGGCCCTGCGGGGGGTTGCCCGCTCCCACCGGTATCCACAGGGTTTTAATGCGGTGATAGCCCGGGGGTGGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGa                                                   <  2:5152089/144‑1 (MQ=255)
                                                                                                 tcatgcggggCCGCTTTGCGGCCCTGCGGGGGCTTGCCCGCTCCCCCTGTTAATCCCTGGGTTTGAATGGTGTGATTTCCCTGTGGGCGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGa                                                   <  1:3617110/144‑1 (MQ=255)
                                                                                                   atgtggggCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGGTATCACAGGGTTTTAAGGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgc                                                 <  2:7761105/146‑1 (MQ=25)
                                                                                                          gCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACg                                          >  2:5119573/1‑150 (MQ=255)
                                                                                                           ccGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                               <  1:7867240/64‑1 (MQ=255)
                                                                                                           ccGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                               >  2:7867240/1‑64 (MQ=255)
                                                                                                           ccGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCGGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGCAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGt                                         >  2:534031/1‑150 (MQ=255)
                                                                                                             gCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCCCTTCTGAAGGGAGCTGGATATGAGCGCCAACGTAg                                       >  2:7880593/1‑150 (MQ=255)
                                                                                                             gCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGGGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGCGCTGGATATGAGCGCCAACGTAg                                       >  2:7880655/1‑150 (MQ=255)
                                                                                                               tttGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGTTAATCCCTGGTTTTTATTGCTGTGATATCCCTGTGGGAGGGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                     <  1:3503866/150‑1 (MQ=38)
                                                                                                               tgggCGGCCCTTCGCGGGCTTGCCCGCTCCCCCCGGTAATCCCAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                     <  2:7023825/147‑1 (MQ=255)
                                                                                                                     gCCCTTCGGGGGCTTGCCCGCTCCCACAGGTAAACACTGGTTTTGAATGCTGTGATAGCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAAc                               <  2:113603/150‑1 (MQ=38)
                                                                                                                         ttCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTGAACGATc                           <  2:2026748/150‑1 (MQ=255)
                                                                                                                                  ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                  >  2:2882792/1‑150 (MQ=255)
                                                                                                                                  ggccccGCCCCCCCCGGAAACCACTGGTTTTGTTTGGGGTGATAGCCGTGTGGGCGGGGGCAACCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                  <  1:5221786/147‑1 (MQ=255)
                                                                                                                                   tGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                 <  2:6046271/150‑1 (MQ=255)
                                                                                                                                   tGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                 >  2:12231/1‑150 (MQ=255)
                                                                                                                                   ggCCCGCGCCCCCAGGTAAACACCGGTTTTGAATGCTGGGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                 <  2:6045785/149‑1 (MQ=39)
                                                                                                                                    gCCCGCTCCCACCGTTAATCACAGGGTTTGAATGCTGTGATAGCCCTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                <  2:449603/150‑1 (MQ=255)
                                                                                                                                    gCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                <  2:4286743/150‑1 (MQ=255)
                                                                                                                                      ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa              >  1:333526/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa              >  2:8211635/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGACCGGATTACGa              >  2:2238621/1‑150 (MQ=255)
                                                                                                                                        gCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc            >  1:2621857/1‑150 (MQ=255)
                                                                                                                                        gCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc            >  1:5132962/1‑150 (MQ=255)
                                                                                                                                            ccACAGGTACTCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCGGAAGGGAGCTGGATATGAGCGCCAACGTCGACCTCAACGATCGCCCGGATTACGACCGGGt        >  2:304781/1‑150 (MQ=255)
                                                                                                                                            ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt        >  1:4497434/1‑150 (MQ=255)
                                                                                                                                            ccACAGGTAATCACAGGGTTTGAATGCTATGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt        >  1:7692562/1‑150 (MQ=255)
                                                                                                                                             cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg       >  2:6976076/1‑150 (MQ=255)
                                                                                                                                             cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg       >  1:3192688/1‑150 (MQ=255)
                                                                                                                                               caGGTAATCACAGGGTTTGAATGCTGTGTTATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCt     <  2:4637556/150‑1 (MQ=255)
                                                                                                                                                  gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACa  >  1:5586221/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TAACCAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACA  >  NC_002947/2672898‑2673193

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: