Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,584,491 C→G 5.6% G175R (GGG→CGG)  PP_1388 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,584,4910CG 5.6% 226.8 / 3.6 71G175R (GGG→CGG) PP_1388EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base C (36/30);  new base G (2/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

CCTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGT  >  NC_002947/1584347‑1584637
                                                                                                                                                |                                                                                                                                                     
cctggcccACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGc                                                                                                                                                   >  3:59271/1‑149 (MQ=255)
           aaGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTTGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGAt                                                                                                                                        <  3:349870/149‑1 (MQ=255)
           aaGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCTGCCGCTGGGGCTGGTTgct                                                                                                                                        <  7:109344/149‑3 (MQ=255)
             ggTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATcc                                                                                                                                      <  4:59271/149‑1 (MQ=255)
                 aTGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGaa                                                                                                                                  <  2:101329/149‑1 (MQ=255)
                     ccAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGTTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGaa                                                                                                                                  <  8:252963/145‑1 (MQ=255)
                     ccAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGTTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGaa                                                                                                                                  >  7:252963/1‑145 (MQ=255)
                     ccAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACAcc                                                                                                                              >  8:123651/1‑149 (MQ=255)
                     ccAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACAcc                                                                                                                              >  3:82906/1‑149 (MQ=255)
                         cagcagGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAgc                                                                                                                          <  4:313357/149‑1 (MQ=255)
                             agGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATc                                                                                                                                       <  3:139283/132‑1 (MQ=255)
                             agGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATc                                                                                                                                       >  4:139283/1‑132 (MQ=255)
                                   ccAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAg                                                                                                                <  4:82906/149‑1 (MQ=255)
                                        gCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACt                                                                                                           >  6:21619/1‑149 (MQ=255)
                                          cGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCAATGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCg                                                                                                         >  4:158453/1‑149 (MQ=255)
                                           gagaATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATc                                                                                                                                       <  5:74379/118‑1 (MQ=255)
                                           gagaATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATc                                                                                                                                       >  6:74379/1‑118 (MQ=255)
                                           gagaATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCgg                                                                                                        <  7:215290/149‑1 (MQ=255)
                                           gagaATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGAGCCAGAACACCCAGCGCCACGACAGGTACTCgg                                                                                                        >  6:291206/1‑149 (MQ=255)
                                                 cAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCCCAAAGCAACCAGCCCCCGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAAcc                                                                                                  >  3:266050/1‑149 (MQ=255)
                                                          ggCGCCAAGTGAGTCCAGCTGTGCGTCACGTCGCTGCAGTGGGATGTCAGCGAGCGCTCGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCGAGTACACCCAGCGCCACGACAGGTACTCGGTCAACCAgcggccg                                                                                          <  4:266050/148‑2 (MQ=255)
                                                              ccAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCAc                                                                                     >  6:320100/1‑149 (MQ=255)
                                                                aaGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTg                                                                                   >  4:301015/1‑149 (MQ=255)
                                                                      gTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGAc                                                                                                                  >  6:234216/1‑112 (MQ=255)
                                                                      gTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGAc                                                                                                                  <  5:234216/112‑1 (MQ=255)
                                                                       tCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGc                                                                            <  5:320100/149‑1 (MQ=255)
                                                                          aCCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCAc                                                                         >  3:94775/1‑149 (MQ=255)
                                                                               tgCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCt                                                                      <  8:182053/147‑1 (MQ=255)
                                                                               tgCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCt                                                                      <  5:21619/147‑1 (MQ=255)
                                                                                gCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGc                                                                   >  5:24837/1‑149 (MQ=255)
                                                                                 ccTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGcc                                                                  >  6:83009/1‑149 (MQ=255)
                                                                                  cTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCcg                                                                 >  7:243198/1‑149 (MQ=255)
                                                                                   tCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACACGGACAACCAGCCGCCAATCACTCTCCCGGCCACGCAGGccac                                                                >  7:165745/1‑147 (MQ=255)
                                                                                          gCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTAc                                                      <  4:35293/149‑1 (MQ=255)
                                                                                              cACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACAtgc                                                  <  8:243709/149‑1 (MQ=255)
                                                                                              cACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACAtg                                                   <  3:301015/148‑1 (MQ=255)
                                                                                               aCTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCACAGCGGCAAGTTGATAAAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCACGGCCACGCTGGCCGAGGA‑‑‑GTTCAtgct                                                 >  1:99660/1‑149 (MQ=255)
                                                                                                   gCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATgctgct                                              >  8:141658/1‑148 (MQ=255)
                                                                                                   gCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATgctgct                                              >  7:211794/1‑148 (MQ=255)
                                                                                                     aTGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACACAGCGCCCCGAAAGGTAATCGGTCAACCAGCCGCCCCGAACTGGCCCGCCCCCGCTGCCCGCCGC‑‑‑GTAAAGgctcctc                                             >  8:305946/1‑146 (MQ=255)
                                                                                                      tGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAg                                          >  8:253423/1‑149 (MQ=255)
                                                                                                      tGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAAAACCCAGCGCCCCGACAGGTACTCGGTCAACCAGCCGCCCCGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTTCTGAAg                                          >  7:167105/1‑149 (MQ=255)
                                                                                                        tCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCa                                                                                      >  8:276558/1‑106 (MQ=255)
                                                                                                        tCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCa                                                                                      <  7:276558/106‑1 (MQ=255)
                                                                                                        tCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAg                                          >  4:178896/1‑147 (MQ=255)
                                                                                                           gcgagcgCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGcc                                     >  6:100944/1‑149 (MQ=255)
                                                                                                            cgagcgCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGccc                                    <  7:123651/149‑1 (MQ=255)
                                                                                                                cgcACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGt                                <  6:24837/149‑1 (MQ=255)
                                                                                                                cgcACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGt                                <  5:83009/149‑1 (MQ=255)
                                                                                                                 gcACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGaa                                           <  1:90045/137‑1 (MQ=255)
                                                                                                                 gcACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGaa                                           >  2:90045/1‑137 (MQ=255)
                                                                                                                 gcACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGc                                      >  5:88725/1‑142 (MQ=255)
                                                                                                                 gcACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGc                                      <  6:88725/142‑1 (MQ=255)
                                                                                                                 gcAAGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTa                               <  7:305946/149‑1 (MQ=255)
                                                                                                                              cccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcg                  >  5:209736/1‑149 (MQ=255)
                                                                                                                                cacaAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcgcg                >  3:340532/1‑149 (MQ=255)
                                                                                                                                 acaAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcgcgc               >  2:166495/1‑149 (MQ=255)
                                                                                                                                    gatcAACCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGaa                                           >  2:235680/4‑118 (MQ=255)
                                                                                                                                    gatcAACCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGaa                                           <  1:235680/115‑1 (MQ=255)
                                                                                                                                     aGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCgggg           <  7:141658/149‑1 (MQ=255)
                                                                                                                                      tcAACCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCggggg          >  7:167861/2‑149 (MQ=255)
                                                                                                                                       cAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGc         <  3:178896/149‑1 (MQ=255)
                                                                                                                                        aaCCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCa        >  6:319245/1‑149 (MQ=255)
                                                                                                                                         aCCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCTAGGGACATGATGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCggggg          >  1:119196/1‑149 (MQ=255)
                                                                                                                                         aCCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCAc       <  7:253423/149‑1 (MQ=255)
                                                                                                                                         aCCAGCCCCAGCGGCAGGTTGAGCCAGAACACCCAGCGCCACGACAGGTACTCGGTCACCCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCAc       <  5:291206/149‑1 (MQ=255)
                                                                                                                                          ccAGCCCCAGCTGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACg      <  8:167105/149‑1 (MQ=255)
                                                                                                                                          ccAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACg      <  5:319245/149‑1 (MQ=255)
                                                                                                                                           cAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACg      >  4:7866/1‑148 (MQ=255)
                                                                                                                                           cAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACg      >  2:121929/1‑148 (MQ=255)
                                                                                                                                            aGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACg      >  8:151532/1‑147 (MQ=255)
                                                                                                                                              ccccAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGt  >  4:305162/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CCTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGC‑‑‑GTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGCACGAAGT  >  NC_002947/1584347‑1584637

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