Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,121,166 G→C 5.2% E324Q (GAG→CAG)  PP_4505 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,121,1660GC 5.2% 715.8 / 20.5 213E324Q (GAG→CAG) PP_4505sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (107/95);  new base C (6/5);  total (113/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 216 total aligned reads displayed.

CCGCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCAGC  >  NC_002947/5121018‑5121315
                                                                                                                                                    |                                                                                                                                                     
ccGCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGTGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGa                                                                                                                                                      <  1:5132148/150‑1 (MQ=255)
ccGCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGa                                                                                                                                                      <  2:433510/150‑1 (MQ=255)
ccGCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGa                                                                                                                                                      <  1:428201/150‑1 (MQ=255)
  gCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAgg                                                                                                                                                    <  2:1475194/150‑1 (MQ=255)
  gCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAgg                                                                                                                                                    >  1:1586647/1‑150 (MQ=255)
    cTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGgcg                                                                                                                                                  <  2:4949471/150‑1 (MQ=255)
    cTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGgcg                                                                                                                                                  >  1:1929891/1‑150 (MQ=255)
    atGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGACCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGgcg                                                                                                                                                  >  2:4042457/2‑150 (MQ=255)
     tGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGgcgc                                                                                                                                                 >  1:3193730/1‑150 (MQ=255)
      ggAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCa                                                                                                                                                >  2:2419141/1‑150 (MQ=255)
       gAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAg                                                                                                                                               >  1:587106/1‑150 (MQ=255)
          gACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCACTGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTCCAGGCGCacgtc                                                                                                                                            <  1:2224517/150‑5 (MQ=255)
           aCCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAg                                                                                                                                           <  1:2426558/150‑1 (MQ=255)
           aCCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAg                                                                                                                                           <  1:4406551/150‑1 (MQ=255)
           aCCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAg                                                                                                                                           <  1:2799035/150‑1 (MQ=255)
             cTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCACGGCCCTCCTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTCCAGGCGCacgtcctg                                                                                                                                         <  1:3262124/150‑8 (MQ=255)
             cTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGt                                                                                                                                         >  1:4331348/1‑150 (MQ=255)
             cTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGt                                                                                                                                         >  2:123133/1‑150 (MQ=255)
             cTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCAGGCCGAACCGGTCGACCTCCATGCCCTCCTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGt                                                                                                                                         <  2:4381536/150‑1 (MQ=255)
               gATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGttt                                                                                                                                       >  1:340742/1‑150 (MQ=255)
                 tCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTCAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGa                                                                                                                                     >  1:5109956/1‑150 (MQ=255)
                   aGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGc                                                                                                                                   >  2:3484433/1‑150 (MQ=255)
                     cGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCtt                                                                                                                                 <  1:394361/150‑1 (MQ=255)
                      gAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTg                                                                                                                                >  1:4082357/1‑150 (MQ=255)
                      gAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTg                                                                                                                                >  2:1386351/1‑150 (MQ=255)
                       aaaTCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGc                                                                                                                               >  1:3559951/1‑150 (MQ=255)
                        aaTCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCa                                                                                                                              <  2:2173468/150‑1 (MQ=255)
                        aaTCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCa                                                                                                                              <  2:293048/150‑1 (MQ=255)
                        aaTCCTCGCCCTGGCCCGGGTGGAGGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCa                                                                                                                              <  2:4179882/150‑1 (MQ=255)
                         aTCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAg                                                                                                                             <  1:2206878/150‑1 (MQ=255)
                             tCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCACTGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGTGCTTGCCGGGtt                                                                                                                         >  1:4580133/1‑150 (MQ=255)
                             tCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGtt                                                                                                                         >  1:988520/1‑150 (MQ=255)
                             tCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGtt                                                                                                                         >  1:4715729/1‑150 (MQ=255)
                             tCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTCCAGGCGCacgtcc                                                                                                                                           <  2:1784958/132‑6 (MQ=255)
                              cGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTg                                                                                                                        <  1:123133/150‑1 (MQ=255)
                                       ccGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCt                                                                                                               >  1:3737348/1‑150 (MQ=255)
                                           gTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATc                                                                                                           <  2:3773347/150‑1 (MQ=255)
                                               aTGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGc                                                                                                       <  2:587106/150‑1 (MQ=255)
                                               aTGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAAGGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGc                                                                                                       <  1:2141816/150‑1 (MQ=255)
                                                    gAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGcc                                                                                                  >  2:1669801/1‑150 (MQ=255)
                                                       cAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGtt                                                                                               <  1:5029550/150‑1 (MQ=255)
                                                        aaGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTg                                                                                              <  1:5392985/150‑1 (MQ=255)
                                                        aaGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTg                                                                                              >  2:2017392/1‑150 (MQ=255)
                                                        aaGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTg                                                                                              >  2:4755217/1‑150 (MQ=255)
                                                          gCCCATGCCGAACCGGTCTACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGAc                                                                                            >  1:282111/1‑150 (MQ=255)
                                                               tGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACct                                                                                       <  2:3711889/150‑1 (MQ=255)
                                                                gCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctg                                                                                      <  2:2877559/150‑1 (MQ=255)
                                                                  cGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgct                                                                                    >  2:1066444/1‑150 (MQ=255)
                                                                  cGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgct                                                                                    <  1:1386351/150‑1 (MQ=255)
                                                                    aaCCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCtgctgc                                                                                  <  2:1966037/150‑1 (MQ=255)
                                                                       cGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCa                                                                               >  2:3790141/1‑150 (MQ=255)
                                                                        ggTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCaa                                                                              <  2:4667011/150‑1 (MQ=255)
                                                                         gTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGa                                                                                                         >  1:1158406/1‑122 (MQ=255)
                                                                            gACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGcc                                                                          >  2:463878/1‑150 (MQ=255)
                                                                            gACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGcc                                                                          <  1:84933/150‑1 (MQ=255)
                                                                            gACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGcc                                                                          >  2:3175597/1‑150 (MQ=255)
                                                                              ccTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCt                                                                        >  2:4067432/1‑150 (MQ=255)
                                                                              ccTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCt                                                                        <  1:2850594/150‑1 (MQ=255)
                                                                              ccTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCt                                                                        <  2:2418714/150‑1 (MQ=255)
                                                                               cTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGa                                                                                                                                     <  1:3175166/88‑1 (MQ=255)
                                                                               cTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTg                                                                       <  2:2821187/150‑1 (MQ=255)
                                                                                tCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgc                                                                      <  1:2017392/150‑1 (MQ=255)
                                                                                tCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgc                                                                      >  2:4829182/1‑150 (MQ=255)
                                                                                tCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTGTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTgc                                                                      >  2:4829095/1‑150 (MQ=255)
                                                                                       ccTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCaa                                                               >  1:1015610/1‑150 (MQ=255)
                                                                                       ccTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCaa                                                               >  2:2230871/1‑150 (MQ=255)
                                                                                         tACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAcc                                                             >  1:540486/1‑150 (MQ=255)
                                                                                          aCTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAccc                                                            <  2:802697/150‑1 (MQ=255)
                                                                                             cGGCAGCGTGCGCACGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCGGCTGCGCAATGCCCTGCGCTTCAACCcggc                                                         <  2:5122901/150‑1 (MQ=255)
                                                                                              ggCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGtt                                                                                               >  1:4857705/1‑111 (MQ=255)
                                                                                                 aGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGc                                                     <  2:1571115/150‑1 (MQ=255)
                                                                                                      gcgcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGcc                                                <  2:3910320/150‑1 (MQ=255)
                                                                                                       cgcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCg                                               >  1:5302457/1‑150 (MQ=255)
                                                                                                        gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCgg                                              >  1:1551425/1‑150 (MQ=255)
                                                                                                         cAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGt                                             >  2:3271193/1‑150 (MQ=255)
                                                                                                          aaGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGTCCTGCGCTTCAACCCGGCGGGCCAGCCGGTc                                            >  2:2814389/1‑150 (MQ=255)
                                                                                                             gACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAg                                         >  1:468990/1‑150 (MQ=255)
                                                                                                             gACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAg                                         <  2:383078/150‑1 (MQ=255)
                                                                                                              aCGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAgg                                        >  2:4154472/1‑150 (MQ=255)
                                                                                                               cGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGt                                       <  2:247102/150‑1 (MQ=255)
                                                                                                                gCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTc                                      >  1:789012/1‑150 (MQ=255)
                                                                                                                  ccTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAg                                    >  2:3118214/1‑150 (MQ=255)
                                                                                                                      ctgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGcgcc                                <  2:540209/150‑1 (MQ=255)
                                                                                                                         aGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCtgctgc                                                                                  >  1:382622/1‑97 (MQ=255)
                                                                                                                          gcgcgcCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGTAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCcgcgc                            <  2:1531018/150‑1 (MQ=255)
                                                                                                                              gcCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGa                        <  1:903118/150‑1 (MQ=255)
                                                                                                                              gcCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGa                        >  2:721987/1‑150 (MQ=255)
                                                                                                                                  gATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAg                    <  1:2583466/150‑1 (MQ=255)
                                                                                                                                    tCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGa                  >  2:3736608/1‑150 (MQ=255)
                                                                                                                                       ggACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGc                                                                                                                      >  2:491257/1‑47 (MQ=255)
                                                                                                                                       ggACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGc                                                                                                                      <  1:491342/47‑1 (MQ=255)
                                                                                                                                       ggACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCg               >  1:2619518/1‑150 (MQ=255)
                                                                                                                                       ggACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCg               >  2:1351448/1‑150 (MQ=255)
                                                                                                                                       ctggctGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCg               >  2:2098333/6‑150 (MQ=255)
                                                                                                                                        tggctGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGc              >  1:2373391/5‑150 (MQ=255)
                                                                                                                                         ggctGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCa             >  2:5016783/4‑150 (MQ=255)
                                                                                                                                             gcgcCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGt         >  1:1386286/1‑150 (MQ=255)
                                                                                                                                               gcCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGa       >  1:3648060/1‑150 (MQ=255)
                                                                                                                                                  tGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCa    >  1:389376/1‑150 (MQ=255)
                                                                                                                                                    gAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCAgc  <  2:5302457/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGCCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCAGC  >  NC_002947/5121018‑5121315

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