Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,806,022 G→T 7.0% T16T (ACG→ACT emrB → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,806,0220GT 7.0% 1205.4 / 42.5 312T16T (ACG→ACTemrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (145/145);  new base T (11/11);  total (156/156)
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.98e-01

Only 100 of 315 total aligned reads displayed.

CGGTAGCGGTAAGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATGCCATCTCGAT  >  CP009273/2805873‑2806171
                                                                                                                                                     |                                                                                                                                                     
 ggTAGCGGTAAGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGa                                                                                                                                                      >  1:4674094/1‑150 (MQ=255)
  gTAGCGGTAAGCCCCGCGCGTGAAATCAGCCTGGCACCTGTCAATAACCTGATCGACGAGATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGAt                                                                                                                                                     <  1:8445/150‑1 (MQ=255)
   tAGCGGTAAGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGAtt                                                                                                                                                    <  1:1327359/150‑1 (MQ=255)
    aGCGGTAAGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTg                                                                                                                                                   <  1:1551292/150‑1 (MQ=255)
    aGCGGTAAGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTg                                                                                                                                                   <  2:4317064/150‑1 (MQ=255)
     gCGGTAAGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTgc                                                                                                                                                  <  2:4951109/150‑1 (MQ=255)
     gCGGTAAGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTgc                                                                                                                                                  <  2:1279848/150‑1 (MQ=255)
        gTAAGCACCGCGCGTGAAATCCGCCTGGCACCTGTCAATCAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCt                                                                                                                                               <  2:3128068/150‑1 (MQ=255)
        gTAAGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCt                                                                                                                                               <  2:1907002/150‑1 (MQ=255)
        gTAAGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCt                                                                                                                                               <  1:187079/150‑1 (MQ=255)
          aaGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGt                                                                                                                                             >  2:2063464/1‑150 (MQ=255)
          aaGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGt                                                                                                                                             >  2:2691771/1‑150 (MQ=255)
          aaGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGt                                                                                                                                             <  2:2414217/150‑1 (MQ=255)
           aGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATACACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTc                                                                                                                                            <  2:3479018/150‑1 (MQ=255)
              aCCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACt                                                                                                                                         >  2:3950868/1‑150 (MQ=255)
               ccGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACTATTGCGCTGTCACTg                                                                                                                                        >  2:4809329/1‑150 (MQ=255)
                 gcgcgTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGc                                                                                                                                      <  1:2017407/150‑1 (MQ=255)
                 gcgcgTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGc                                                                                                                                      <  2:4156986/150‑1 (MQ=255)
                  cgcgTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACTATTGCGCTGTCACTGGCg                                                                                                                                     <  1:4502709/150‑1 (MQ=255)
                  cgcgTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCg                                                                                                                                     <  1:1217736/150‑1 (MQ=255)
                  cgcgTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCg                                                                                                                                     <  2:1767369/150‑1 (MQ=255)
                    cgTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGAc                                                                                                                                   >  2:624781/1‑150 (MQ=255)
                     gTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTCACGCTGGCTAATCCAGAGGTGCGTGTGCTGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACa                                                                                                                                  <  2:3411560/150‑1 (MQ=255)
                         aaTCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACTATTGCGCTGTCACTGGCGACATTCa                                                                                                                              >  1:4601372/1‑150 (MQ=255)
                         aaTCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCa                                                                                                                              >  2:3860688/1‑150 (MQ=255)
                             aGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCACCAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCa                                                                                                                          >  2:3709697/1‑150 (MQ=255)
                                cTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGt                                                                                                                       <  1:3866280/150‑1 (MQ=255)
                                 tGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTg                                                                                                                      >  1:3734534/1‑150 (MQ=255)
                                   gCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCt                                                                                                                    >  2:807707/1‑150 (MQ=255)
                                   gCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCt                                                                                                                    <  1:2647487/150‑1 (MQ=255)
                                     aCCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTgg                                                                                                                  <  1:2281045/150‑1 (MQ=255)
                                       cTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGAc                                                                                                                <  2:3644014/150‑1 (MQ=255)
                                        tGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACt                                                                                                               >  1:1744047/1‑150 (MQ=255)
                                        tGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACt                                                                                                               <  1:1143018/150‑1 (MQ=255)
                                         gTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTc                                                                                                              >  1:1097524/1‑150 (MQ=255)
                                             aTAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTccacc                                                                                                          <  2:3439987/150‑1 (MQ=255)
                                             aTAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTccacc                                                                                                          >  1:3194731/1‑150 (MQ=255)
                                               aaaCTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCAt                                                                                                        <  1:4363660/150‑1 (MQ=255)
                                                aaCTGATCGACGATATAGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCAtt                                                                                                       >  2:3838207/1‑150 (MQ=255)
                                                aaCGGATCGACGATATCGTAAAAGCTACCGCGGGCTACTCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCAtt                                                                                                       <  2:4278649/150‑1 (MQ=255)
                                                  cTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGc                                                                                                     <  1:5028962/150‑1 (MQ=255)
                                                      tCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAAc                                                                                                 <  1:1725223/150‑1 (MQ=255)
                                                      tCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAAc                                                                                                 >  1:3961958/1‑150 (MQ=255)
                                                       cgacgaTATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACg                                                                                                <  1:4684770/150‑1 (MQ=255)
                                                        gacgaTATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGt                                                                                               >  1:2918356/1‑150 (MQ=255)
                                                          cgaTATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTgg                                                                                             >  2:3777451/1‑150 (MQ=255)
                                                            atatCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCg                                                                                           <  2:3056332/150‑1 (MQ=255)
                                                                         tAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCg                                                                              <  2:5120082/150‑1 (MQ=255)
                                                                          aaCGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACTATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGc                                                                             >  1:2162479/1‑150 (MQ=255)
                                                                           aCGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGcc                                                                            <  2:2591543/150‑1 (MQ=255)
                                                                            cGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACTATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCg                                                                           <  1:601012/150‑1 (MQ=255)
                                                                            cGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCg                                                                           >  1:4966554/1‑150 (MQ=255)
                                                                            cGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCg                                                                           >  1:3786242/1‑150 (MQ=255)
                                                                             gCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCgg                                                                          >  1:3788095/1‑150 (MQ=255)
                                                                              cTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTAGCGCCggg                                                                         <  2:4736622/150‑1 (MQ=255)
                                                                                ggCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGaa                                                                       >  2:112341/1‑150 (MQ=255)
                                                                                 gCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAAt                                                                      <  1:3426890/150‑1 (MQ=255)
                                                                                  cTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATc                                                                     >  1:2481275/1‑150 (MQ=255)
                                                                                     aTCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTgg                                                                  <  1:3462197/150‑1 (MQ=255)
                                                                                      tCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTggg                                                                 >  1:2696987/1‑150 (MQ=255)
                                                                                       ccAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGc                                                                >  2:4910939/1‑150 (MQ=255)
                                                                                       ccAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGc                                                                <  2:4974926/150‑1 (MQ=255)
                                                                                         agagGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTc                                                              >  2:2810951/1‑150 (MQ=255)
                                                                                          gagGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCa                                                             >  1:27771/1‑150 (MQ=255)
                                                                                            ggTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATc                                                           >  1:3599215/1‑150 (MQ=255)
                                                                                              tGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGc                                                         >  1:1882978/1‑150 (MQ=255)
                                                                                                   gtgATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGc                                                    >  2:5065142/1‑150 (MQ=255)
                                                                                                    tgATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACTATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGcc                                                   >  1:2568629/1‑150 (MQ=255)
                                                                                                       tGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAgg                                                >  2:460069/1‑150 (MQ=255)
                                                                                                       tGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAgg                                                <  2:4573282/150‑1 (MQ=255)
                                                                                                         cAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGa                                              >  2:4084725/1‑150 (MQ=255)
                                                                                                         cAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGa                                              <  2:4617348/150‑1 (MQ=255)
                                                                                                          aaCAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGaa                                             >  1:578192/1‑150 (MQ=255)
                                                                                                          aaCAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGaa                                             >  2:1335517/1‑150 (MQ=255)
                                                                                                            cAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACg                                           <  2:4014155/150‑1 (MQ=255)
                                                                                                             aGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGt                                          <  2:2702319/150‑1 (MQ=255)
                                                                                                              gCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTg                                         <  2:3370840/150‑1 (MQ=255)
                                                                                                               cAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTgg                                        <  1:4108409/150‑1 (MQ=255)
                                                                                                                 aaaaCCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGt                                      >  1:2345058/1‑150 (MQ=255)
                                                                                                                  caaCCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTa                                     <  2:1209414/149‑1 (MQ=255)
                                                                                                                     ccGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATc                                  <  1:105560/150‑1 (MQ=255)
                                                                                                                      cGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCa                                 <  2:1573703/150‑1 (MQ=255)
                                                                                                                        cTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCAct                               <  2:1718565/150‑1 (MQ=255)
                                                                                                                        cTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCAct                               <  1:5073129/150‑1 (MQ=255)
                                                                                                                         tGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCActt                              >  2:3236633/1‑150 (MQ=255)
                                                                                                                          ggAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCActtc                             <  2:1427611/150‑1 (MQ=255)
                                                                                                                             aGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTGttt                          <  2:2114020/150‑1 (MQ=255)
                                                                                                                             aGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCttt                          <  1:2401692/150‑1 (MQ=255)
                                                                                                                                 gcgcAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCggg                      >  2:1661376/1‑150 (MQ=255)
                                                                                                                                    cAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTg                   >  2:4059684/1‑150 (MQ=255)
                                                                                                                                       cTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCg                <  2:2951964/150‑1 (MQ=255)
                                                                                                                                        tGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGa               <  1:687111/150‑1 (MQ=255)
                                                                                                                                         ggTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGaa              >  2:2284030/1‑150 (MQ=255)
                                                                                                                                           tCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATg            <  1:984058/150‑1 (MQ=255)
                                                                                                                                           tCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATg            >  2:4347709/1‑150 (MQ=255)
                                                                                                                                            cATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATGc           <  1:1166152/150‑1 (MQ=255)
                                                                                                                                            cATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATGc           >  1:2110309/1‑150 (MQ=255)
                                                                                                                                               tATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCTCTTCTTTCGGGGTGGCGAATGCCAt        >  2:729929/1‑150 (MQ=255)
                                                                                                                                               tATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATGCCAt        >  1:3070201/1‑150 (MQ=255)
                                                                                                                                                aTGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATGCCAtc       >  1:4001946/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGTAGCGGTAAGCACCGCGCGTGAAATCAGCCTGGCACCTGTCAATAAACTGATCGACGATATCGTAAAAGCTAACGCTGGCTAATCCAGAGGTGCGTGTGATGCAACAGCAAAAACCGCTGGAAGGCGCGCAACTGGTCATTATGACGATTGCGCTGTCACTGGCGACATTCATGCAGGTGCTGGACTCCACCATTGCTAACGTGGCGATCCCCACTATCGCCGGGAATCTGGGCTCATCGCTCAGCCAGGGAACGTGGGTAATCACTTCTTTCGGGGTGGCGAATGCCATCTCGAT  >  CP009273/2805873‑2806171

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: