Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,243,136 A→T 5.2% M147K (ATG→AAG) ‡ PP_1975 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,243,1360AT 5.2% 803.2 / 21.9 232M147K (ATG→AAG) ‡PP_1975EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base A (96/123);  new base T (9/3);  total (105/127)
Fisher's exact test for biased strand distribution p-value = 4.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 234 total aligned reads displayed.

TCTTCTCAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGCACAGCAGCGA  >  NC_002947/2242987‑2243284
                                                                                                                                                     |                                                                                                                                                    
tcttctCAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACa                                                                                                                                                      >  2:4894157/1‑150 (MQ=255)
 cttctCAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACAt                                                                                                                                                     <  1:736502/150‑1 (MQ=255)
     tcAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGcc                                                                                                                                                 <  2:3426461/150‑1 (MQ=255)
     tcAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGcc                                                                                                                                                 >  2:4789448/1‑150 (MQ=255)
      cAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCa                                                                                                                                                <  1:430558/150‑1 (MQ=255)
       aaGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCAt                                                                                                                                               <  1:1558470/150‑1 (MQ=255)
       aaGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGCCATGGCGATGTTCGCCAt                                                                                                                                               <  1:4988550/150‑1 (MQ=255)
         gCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGc                                                                                                                                             <  2:5390440/150‑1 (MQ=255)
         gCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGc                                                                                                                                             <  2:730425/150‑1 (MQ=255)
           ccGTATAACAAGCCCGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCt                                                                                                                                           >  1:4351136/1‑150 (MQ=255)
            cGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCtt                                                                                                                                          <  1:3708019/150‑1 (MQ=255)
              tataACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTg                                                                                                                                        <  1:5417096/150‑1 (MQ=255)
              tataACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTg                                                                                                                                        <  2:1940011/150‑1 (MQ=255)
              tataACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTg                                                                                                                                        <  1:3980896/150‑1 (MQ=255)
              gataACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTg                                                                                                                                        <  1:5252808/149‑1 (MQ=255)
                  aCAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCGGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGgcgc                                                                                                                                    >  1:4016441/1‑150 (MQ=255)
                     aggcGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCggt                                                                                                                                 <  1:176798/147‑1 (MQ=255)
                      cccGGCGATCATTAGCAGGCCCGGTGGTCTATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCTGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAATAGTGGCTGTCATGGCGAACATCGCCATGCATTTGTCGCggtg                                                                                                                                <  2:4970807/149‑1 (MQ=255)
                        cGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCggtggt                                                                                                                              <  2:1800099/150‑1 (MQ=255)
                            gATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTc                                                                                                                          >  1:4690519/1‑150 (MQ=255)
                            gATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTc                                                                                                                          <  1:455849/150‑1 (MQ=255)
                              tCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCgg                                                                                                                        >  1:1636519/1‑150 (MQ=255)
                              tCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCgg                                                                                                                        <  1:5191680/150‑1 (MQ=255)
                               cATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCggg                                                                                                                       <  2:4661677/150‑1 (MQ=255)
                               cATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCggg                                                                                                                       <  2:628529/150‑1 (MQ=255)
                                   aGCCGGCCCGGTGGGAGATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGc                                                                                                                   <  2:4882678/150‑1 (MQ=255)
                                   aGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGc                                                                                                                   <  1:4589610/150‑1 (MQ=255)
                                    gCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCt                                                                                                                  >  1:7499/1‑150 (MQ=255)
                                      aGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTg                                                                                                                <  2:1266093/150‑1 (MQ=255)
                                       ggCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGa                                                                                                               <  2:203647/150‑1 (MQ=255)
                                        gCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGAt                                                                                                              <  1:1185518/150‑1 (MQ=255)
                                            ggTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGc                                                                                                          <  1:3049387/150‑1 (MQ=255)
                                                    atTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGt                                                                                                  <  1:1984169/150‑1 (MQ=255)
                                                       gATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCggt                                                                                                                                 >  1:2688971/1‑116 (MQ=255)
                                                       gATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCggt                                                                                                                                 <  2:2688971/116‑1 (MQ=255)
                                                        aTGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAg                                                                                              <  1:689056/150‑1 (MQ=255)
                                                          gTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGc                                                                                            >  2:3869079/1‑150 (MQ=255)
                                                          gTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGc                                                                                            >  2:1886373/1‑150 (MQ=255)
                                                             gaagaTGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAg                                                                                         >  1:1032229/1‑150 (MQ=255)
                                                             gaagaTGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAg                                                                                         >  2:723555/1‑150 (MQ=255)
                                                             gaagaTGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAg                                                                                         >  2:5284939/1‑150 (MQ=255)
                                                                gaTGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCgg                                                                                      <  1:1612647/150‑1 (MQ=255)
                                                                gaTGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCgg                                                                                      >  1:2862060/1‑150 (MQ=255)
                                                                 aTGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGACAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCggg                                                                                     <  1:5078403/150‑1 (MQ=255)
                                                                   gTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGAt                                                                                   <  1:995902/150‑1 (MQ=255)
                                                                   gTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGAt                                                                                   >  2:2761681/1‑150 (MQ=255)
                                                                     aTTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCa                                                                                 <  1:3287865/150‑1 (MQ=255)
                                                                     aTTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCa                                                                                 <  1:5042018/150‑1 (MQ=255)
                                                                      ttCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGa                                                                                    >  1:1871852/1‑146 (MQ=255)
                                                                      ttCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAg                                                                                >  2:4539910/1‑150 (MQ=255)
                                                                       tCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAgc                                                                               >  1:4233402/1‑150 (MQ=255)
                                                                         ccAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAGGGCCAg                                                                                         >  2:1091631/1‑138 (MQ=255)
                                                                         ccAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAGGGCCAg                                                                                         <  1:1091631/138‑1 (MQ=255)
                                                                         ccAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAgcgc                                                                             >  1:2541947/1‑150 (MQ=255)
                                                                         ccAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAgcgc                                                                             >  2:2659263/1‑150 (MQ=255)
                                                                         ccAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAgcgc                                                                             >  2:2517325/1‑150 (MQ=255)
                                                                              cccAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGGACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCAc                                                                        >  2:4412522/1‑150 (MQ=255)
                                                                               ccAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCAcc                                                                       >  1:2887056/1‑150 (MQ=255)
                                                                               ccAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCAcc                                                                       >  1:1949150/1‑150 (MQ=255)
                                                                                cAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCg                                                                      <  1:1488456/150‑1 (MQ=255)
                                                                                  cttCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGt                                                                    <  2:2867026/149‑1 (MQ=255)
                                                                                          tCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACCTCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCt                                                            <  2:4751973/150‑1 (MQ=255)
                                                                                          tCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCt                                                            <  2:315324/150‑1 (MQ=255)
                                                                                               cagccGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAg                                                       <  2:1454648/150‑1 (MQ=255)
                                                                                                agccGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAgg                                                      >  2:504102/1‑150 (MQ=255)
                                                                                                 gccgccCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGc                                                     >  2:3334333/1‑150 (MQ=255)
                                                                                                  ccgccCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGc                                                                                        >  2:2990639/1‑114 (MQ=255)
                                                                                                   cgccCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCAc                                                   <  2:2380565/150‑1 (MQ=255)
                                                                                                   cgccCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCAc                                                   <  1:1302471/150‑1 (MQ=255)
                                                                                                    gccCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACg                                                  >  2:4254897/1‑150 (MQ=255)
                                                                                                     cccAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACgg                                                 >  2:82413/1‑150 (MQ=255)
                                                                                                      ccAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGa                                                <  2:3960070/150‑1 (MQ=255)
                                                                                                      ccAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGa                                                <  2:2861508/150‑1 (MQ=255)
                                                                                                         gggtgggtCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACa                                             <  1:4182222/150‑1 (MQ=255)
                                                                                                                 ccGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATg                                     <  1:1041490/150‑1 (MQ=255)
                                                                                                                 ccGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATg                                     <  1:2039290/150‑1 (MQ=255)
                                                                                                                   gATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGACCAGGATCATGct                                   >  2:5530667/1‑150 (MQ=255)
                                                                                                                   gATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGct                                   >  1:1722097/1‑150 (MQ=255)
                                                                                                                   gATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGct                                   >  2:4109511/1‑150 (MQ=255)
                                                                                                                       gAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTcc                               <  2:889664/150‑1 (MQ=255)
                                                                                                                        aGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCa                              <  1:849722/150‑1 (MQ=255)
                                                                                                                        aGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCa                              <  1:438310/150‑1 (MQ=255)
                                                                                                                        aGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCa                              <  1:2105990/150‑1 (MQ=255)
                                                                                                                          ggggCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAAAAGGATCATGCTCTCCAgg                            >  1:4814195/1‑150 (MQ=255)
                                                                                                                             gCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTc                         >  2:1748194/1‑150 (MQ=255)
                                                                                                                               aaaaGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCa                       >  2:380711/1‑150 (MQ=255)
                                                                                                                                aaaGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAg                      <  1:2499045/150‑1 (MQ=255)
                                                                                                                                     ggCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAg                      >  1:3921442/1‑145 (MQ=255)
                                                                                                                                      gCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATc                <  2:3496999/150‑1 (MQ=255)
                                                                                                                                           cATGGCGATGTTCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGca           >  2:2966494/1‑150 (MQ=255)
                                                                                                                                            aTGGCGATGTTCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGcac          >  1:3176822/1‑150 (MQ=255)
                                                                                                                                            aTGGCGATGTTCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGcac          >  1:4111029/1‑150 (MQ=255)
                                                                                                                                            aTGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGcac          <  2:5080990/150‑1 (MQ=255)
                                                                                                                                             tGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCCGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGACCAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGcaca         >  1:4495792/1‑150 (MQ=255)
                                                                                                                                             tGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGcaca         <  1:2831910/150‑1 (MQ=255)
                                                                                                                                             tGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGcaca         >  1:2786250/1‑150 (MQ=255)
                                                                                                                                             tGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGcaca         <  1:1158078/150‑1 (MQ=255)
                                                                                                                                              ggCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGCAcag        >  1:2987430/1‑150 (MQ=255)
                                                                                                                                              ggCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGCAcag        >  1:2170319/1‑150 (MQ=255)
                                                                                                                                                  aaCATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGCACagcagc    <  2:1475388/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTTCTCAAGCCCGTATAACAAGCCGGCGATCATTAGCAGGCCCGGTGGGATATTGATGTAGAAGATGTATTCCCAGCCCCAGTTCTCGGTCAGCCAGCCGCCCAGGGTGGGTCCGATGGAGGGGGCAAAAGTGGCGGTCATGGCGAACATCGCCATGCCTTTGGCGCGGTGGTGCTCGGGCAGCTTGATCAGCGTCAGGGTGAAGGCCAGCGGGATCAGCGCGCCACCGGTGAAACCCTGCAAGGCACGGAACAGGATCATGCTCTCCAGGTTCCAGGCCATCGAGCACAGCAGCGA  >  NC_002947/2242987‑2243284

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: