Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,610,676 G→A E14K (GAA→AAA)  emrE → small multidrug resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,610,6760GA99.5% 1328.3 / ‑5.3 366E14K (GAA→AAA) emrEsmall multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (178/186);  total (180/186)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 370 total aligned reads displayed.

CAAGGTTCCACTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATC  >  NC_002947/5610529‑5610822
                                                                                                                                                   |                                                                                                                                                  
    gTTCCACTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATGACTTTACACCATCCAGGAAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCAAGAATGCCTATACCTACCTCGCCATTGCCATCAGCGCCAACGTCa                                                                                                                                              >  1:2144794/1‑150 (MQ=255)
     ttCCACTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCAt                                                                                                                                             >  2:5116447/1‑150 (MQ=255)
      tCCACTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATc                                                                                                                                            <  1:7120268/150‑1 (MQ=255)
      tCCACTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATc                                                                                                                                            <  2:864190/150‑1 (MQ=255)
      tCCACTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATc                                                                                                                                            <  1:3341406/150‑1 (MQ=255)
      gccACTGCTCGCCCATCAACCGATCTCCCATCTCGCTAGGTCATGACTTTACCCCAACCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTAGACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATc                                                                                                                                            <  2:5643793/149‑1 (MQ=255)
         aCTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGcc                                                                                                                                         >  2:2774617/1‑150 (MQ=255)
         aCTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGcc                                                                                                                                         >  1:6815055/1‑150 (MQ=255)
         aCTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCGAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGcc                                                                                                                                         >  1:2473483/1‑150 (MQ=255)
                    aTCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGAt                                                                                                                              <  2:2533694/150‑1 (MQ=255)
                        aCCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAg                                                                                                                          >  2:1530079/1‑150 (MQ=255)
                        aCCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAg                                                                                                                          >  1:54848/1‑150 (MQ=255)
                         ccTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAgg                                                                                                                         <  1:4958780/150‑1 (MQ=255)
                         ccTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAgg                                                                                                                         >  2:4371081/1‑150 (MQ=255)
                          cTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTACCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGc                                                                                                                        >  2:7875427/1‑150 (MQ=255)
                           tATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCTCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGACAAAGTCATCGCCACCGCCTCGATGAAGGCg                                                                                                                       <  2:1734250/150‑1 (MQ=255)
                            aTCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTTCATCGAGGCgg                                                                                                                      <  2:2650263/150‑1 (MQ=255)
                            aTCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCgg                                                                                                                      <  1:2713310/150‑1 (MQ=255)
                             tctcCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTTCATCGAGGCGGt                                                                                                                     <  2:5140356/150‑1 (MQ=255)
                             tctcCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGt                                                                                                                     <  1:7532497/150‑1 (MQ=255)
                              ctcCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTg                                                                                                                    >  1:7460609/1‑150 (MQ=255)
                               tcCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGa                                                                                                                   >  2:5944203/1‑150 (MQ=255)
                                cccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTATCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTTCATCGAGGCGGTggc                                                                                                                  <  2:4490165/150‑3 (MQ=255)
                                cccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTggc                                                                                                                  <  1:4802160/150‑3 (MQ=255)
                                 ccATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGaaa                                                                                                                 <  1:5909349/150‑1 (MQ=255)
                                  cATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTTCATCGAGGCGGTggcga                                                                                                                <  2:3273705/150‑5 (MQ=255)
                                  cATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAg                                                                                                                <  1:3932533/150‑1 (MQ=255)
                                      tCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCt                                                                                                            >  2:7181406/1‑150 (MQ=255)
                                      tCGAGAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCt                                                                                                            <  2:5771281/150‑1 (MQ=255)
                                        gagaGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGa                                                                                                          <  2:4755149/147‑1 (MQ=255)
                                        gATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGa                                                                                                          >  2:6118931/1‑150 (MQ=255)
                                          tAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGc                                                                                                        <  1:7732268/150‑1 (MQ=255)
                                              tAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGc                                                                                                    >  2:6542855/1‑150 (MQ=255)
                                               aaTGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGcc                                                                                                   >  2:2611788/1‑150 (MQ=255)
                                                aTGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCg                                                                                                  >  1:1613361/1‑150 (MQ=255)
                                                  gACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGtt                                                                                                <  2:5859725/150‑1 (MQ=255)
                                                    cTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGc                                                                                              <  2:7752230/150‑1 (MQ=255)
                                                    cTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGc                                                                                              <  2:5404654/150‑1 (MQ=255)
                                                       tACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCAc                                                                                           <  1:7316901/150‑1 (MQ=255)
                                                        aCCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCAct                                                                                          <  2:5713143/150‑1 (MQ=255)
                                                           ccATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCActgct                                                                                       >  2:808509/1‑150 (MQ=255)
                                                           ccATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCActgct                                                                                       >  2:2814180/1‑150 (MQ=255)
                                                            cATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCActgctg                                                                                      >  2:7590561/1‑150 (MQ=255)
                                                                cAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGa                                                                                  <  1:5607991/150‑1 (MQ=255)
                                                                cAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGa                                                                                  <  2:2921472/150‑1 (MQ=255)
                                                                cAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGa                                                                                  <  1:1989872/150‑1 (MQ=255)
                                                                 aGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAt                                                                                 >  1:5674243/1‑150 (MQ=255)
                                                                  ggCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtg                                                                                >  2:4612923/1‑150 (MQ=255)
                                                                   gCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtgg                                                                               >  2:6059977/1‑150 (MQ=255)
                                                                   gCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtgg                                                                               >  1:6169110/1‑150 (MQ=255)
                                                                     aaCGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAtggtg                                                                             <  2:5275364/150‑1 (MQ=255)
                                                                         tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTg                                                                         >  1:1094784/1‑150 (MQ=255)
                                                                         tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTg                                                                         >  1:4733398/1‑150 (MQ=255)
                                                                         tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTg                                                                         >  1:5314002/1‑150 (MQ=255)
                                                                           tAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGc                                                                       >  2:647639/1‑150 (MQ=255)
                                                                           tAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGc                                                                       >  1:6378806/1‑150 (MQ=255)
                                                                              aaTGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTAt                                                                    <  2:1475405/150‑1 (MQ=255)
                                                                                   ccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCAt                                                               <  2:4231894/150‑1 (MQ=255)
                                                                                    cGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCAtt                                                              <  2:2014649/150‑1 (MQ=255)
                                                                                    cGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCAtt                                                              >  1:81174/1‑150 (MQ=255)
                                                                                     gCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTg                                                             >  1:4443577/1‑150 (MQ=255)
                                                                                           tAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCt                                                       >  2:6432057/1‑150 (MQ=255)
                                                                                           tAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCt                                                       <  2:5497461/150‑1 (MQ=255)
                                                                                           tAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCt                                                       <  1:6478104/150‑1 (MQ=255)
                                                                                           tAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCt                                                       <  2:1020537/150‑1 (MQ=255)
                                                                                            aaCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTg                                                      <  1:1721621/150‑1 (MQ=255)
                                                                                            aaCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTg                                                      <  1:7782989/150‑1 (MQ=255)
                                                                                            aaCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTg                                                      <  1:4108663/150‑1 (MQ=255)
                                                                                               ggaggaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGAt                                                   <  2:310222/150‑1 (MQ=255)
                                                                                                   gaTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTg                                               <  1:3437005/150‑1 (MQ=255)
                                                                                                    aTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGa                                              >  2:2393326/1‑150 (MQ=255)
                                                                                                    aTCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGa                                              >  2:3462970/1‑150 (MQ=255)
                                                                                                       cTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGc                                           >  2:3668714/1‑150 (MQ=255)
                                                                                                        tGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCt                                          >  1:3707536/1‑150 (MQ=255)
                                                                                                        tGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCt                                          >  1:1736509/1‑150 (MQ=255)
                                                                                                         gCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtg                                         <  2:6394967/150‑1 (MQ=255)
                                                                                                           cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggt                                       >  2:3998771/1‑150 (MQ=255)
                                                                                                            aTGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtg                                      >  1:1868870/1‑150 (MQ=255)
                                                                                                             tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtgg                                     >  2:3365534/1‑150 (MQ=255)
                                                                                                                aTGCCTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTcc                                  >  2:7238451/1‑150 (MQ=255)
                                                                                                                   ccTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgca                               <  1:3873057/150‑1 (MQ=255)
                                                                                                                    cTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                              >  1:6936321/1‑150 (MQ=255)
                                                                                                                    cTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                              >  1:5143658/1‑150 (MQ=255)
                                                                                                                    cTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                              >  1:3169493/1‑150 (MQ=255)
                                                                                                                    cTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                              >  1:2919559/1‑150 (MQ=255)
                                                                                                                    cTATACCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                              >  1:1393523/1‑150 (MQ=255)
                                                                                                                    cTATACCTACCTCGCCATTGCCATCTGCCCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                              >  1:1185948/1‑150 (MQ=255)
                                                                                                                      ataCCTACCTCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcagca                            <  2:6815055/150‑1 (MQ=255)
                                                                                                                              ctCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTg                    <  2:53055/150‑1 (MQ=255)
                                                                                                                              ctCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTg                    <  2:230880/150‑1 (MQ=255)
                                                                                                                               tCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                   >  2:7417825/1‑150 (MQ=255)
                                                                                                                               tCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                   >  1:869235/1‑150 (MQ=255)
                                                                                                                               tCGCCATTGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                   >  1:7398583/1‑150 (MQ=255)
                                                                                                                                      tGCCATCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGc            <  2:3814183/150‑1 (MQ=255)
                                                                                                                                          aTCTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTAc        >  2:1566296/1‑150 (MQ=255)
                                                                                                                                            cTGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGc      >  2:6246662/1‑150 (MQ=255)
                                                                                                                                             tGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGcc     >  2:7125335/1‑150 (MQ=255)
                                                                                                                                             tGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGcc     >  2:492613/1‑150 (MQ=255)
                                                                                                                                             tGCGCCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGcc     <  2:7460609/150‑1 (MQ=255)
                                                                                                                                                gcCAAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATc  <  1:6613433/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAAGGTTCCACTGCTCGGCCATCAACCTATCTCCCATATCGATAGGTAATGACTTTACCCCATCCAGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATC  >  NC_002947/5610529‑5610822

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: