Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,610,742 +A coding (106/333 nt) emrE → small multidrug resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,610,7411.A98.7% 816.2 / 1.8 231coding (105/333 nt)emrEsmall multidrug resistance protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (116/112);  total (118/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

AGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCTATGGGCAGAAGCTG  >  NC_002947/5610594‑5610885
                                                                                                                                                    |                                                                                                                                                
aGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCCATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGaa                                                                                                                                                >  2:3597029/1‑149 (MQ=255)
    aaCGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGaatggt                                                                                                                                            >  2:4933317/1‑145 (MQ=255)
    aaCGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGaatggt                                                                                                                                            >  1:773554/1‑145 (MQ=255)
      cGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAAtggtg                                                                                                                                            <  2:208845/149‑1 (MQ=255)
       gTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATggtggt                                                                                                                                          >  1:4424402/1‑150 (MQ=255)
       gTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATggtggt                                                                                                                                          >  1:4827474/1‑150 (MQ=255)
        tCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGtt                                                                                                                                         <  2:773554/150‑1 (MQ=255)
          tAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGACAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTgg                                                                                                                                       >  2:3277582/1‑150 (MQ=255)
             aaTGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTa                                                                                                                                    >  2:2391620/1‑150 (MQ=255)
              aTGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTAt                                                                                                                                   <  2:4601007/150‑1 (MQ=255)
              aTGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTAt                                                                                                                                   <  1:323505/150‑1 (MQ=255)
                  ccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCa                                                                                                                               >  2:3676705/1‑150 (MQ=255)
                  ccGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCa                                                                                                                               >  1:3621927/1‑150 (MQ=255)
                   cGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCAt                                                                                                                              >  1:800090/1‑150 (MQ=255)
                    gCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGACAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCAtt                                                                                                                             <  2:1033726/150‑1 (MQ=255)
                    gCTCAGTAACGAAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCAtt                                                                                                                             <  1:2950016/150‑1 (MQ=255)
                      tCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGc                                                                                                                           >  1:1901032/1‑150 (MQ=255)
                        aGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGt                                                                                                                         <  2:3342769/150‑1 (MQ=255)
                          tAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTc                                                                                                                       <  1:42475/150‑1 (MQ=255)
                           aaCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCt                                                                                                                      <  2:436426/150‑1 (MQ=255)
                           aaCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCt                                                                                                                      <  2:2207121/150‑1 (MQ=255)
                           aaCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCt                                                                                                                      <  1:1090388/150‑1 (MQ=255)
                            aCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTg                                                                                                                     <  2:4436625/150‑1 (MQ=255)
                            aCGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTg                                                                                                                     >  1:259576/1‑150 (MQ=255)
                             cGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTgg                                                                                                                    <  1:3473376/150‑1 (MQ=255)
                                    tCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGa                                                                                                             >  1:4860766/1‑150 (MQ=255)
                                       tGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGc                                                                                                          >  2:2750237/1‑150 (MQ=255)
                                          cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtgg                                                                                                       <  2:2148777/150‑1 (MQ=255)
                                          cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtgg                                                                                                       <  1:3949217/150‑1 (MQ=255)
                                          cATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtgg                                                                                                       >  1:2164065/1‑150 (MQ=255)
                                           aTGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGGTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggt                                                                                                      >  2:589227/1‑150 (MQ=255)
                                           aTGAATGCATATACCTACCTCGCCATTGCCATCTGCTCCGAAGTCATAGCCACCTCCTCGATGAAGGCGGTGAAAGTCCTGAGCACGCAGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggt                                                                                                      >  2:1313715/1‑150 (MQ=255)
                                            tGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCtggtg                                                                                                     <  1:3557326/150‑1 (MQ=255)
                                               aTGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTc                                                                                                  <  1:2276077/150‑1 (MQ=255)
                                                 gCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCg                                                                                                >  1:1965268/1‑150 (MQ=255)
                                                   cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgca                                                                                              <  1:3458683/150‑1 (MQ=255)
                                                   cTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgca                                                                                              >  2:1603035/1‑150 (MQ=255)
                                                    tataCCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcag                                                                                             >  1:2376354/1‑150 (MQ=255)
                                                     ataCCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgcagc                                                                                            >  2:5210544/1‑150 (MQ=255)
                                                       acctacctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCAt                                                                                          >  2:1890483/1‑150 (MQ=255)
                                                           acctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCg                                                                                      <  1:4565401/150‑1 (MQ=255)
                                                           acctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCg                                                                                      >  2:515471/1‑150 (MQ=255)
                                                            cctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCgg                                                                                     <  2:4124737/150‑1 (MQ=255)
                                                             ctCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGt                                                                                    >  1:2485356/1‑150 (MQ=255)
                                                              tCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTg                                                                                   >  2:3733930/1‑150 (MQ=255)
                                                              gcGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTg                                                                                   <  2:4757169/149‑1 (MQ=255)
                                                               cGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgg                                                                                  <  1:2307006/150‑1 (MQ=255)
                                                                gCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTggg                                                                                 >  1:3375768/1‑150 (MQ=255)
                                                                gCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTggg                                                                                 >  2:2474066/1‑150 (MQ=255)
                                                                 ccATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTgggg                                                                                <  1:2675945/150‑1 (MQ=255)
                                                                        cATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCt                                                                         >  2:3584100/1‑150 (MQ=255)
                                                                         aTCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTa                                                                        <  2:4250131/150‑1 (MQ=255)
                                                                                ccGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATc                                                                 >  1:4817036/1‑150 (MQ=255)
                                                                                 cGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCt                                                                                                         <  2:1388237/109‑1 (MQ=255)
                                                                                 cGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCt                                                                                                         >  1:1388237/1‑109 (MQ=255)
                                                                                  gAAGTCATCGCCACCGCCTCGATGAAGGCGGTGACAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTg                                                               >  2:4192330/1‑150 (MQ=255)
                                                                                  gAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTg                                                               >  1:4274741/1‑150 (MQ=255)
                                                                                      tCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGt                                                                                    <  2:3183669/125‑1 (MQ=255)
                                                                                      tCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGt                                                                                    >  1:3183669/1‑125 (MQ=255)
                                                                                        aTCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCa                                                                   <  2:296572/140‑1 (MQ=255)
                                                                                        aTCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCggg                                                         >  2:1094351/1‑150 (MQ=255)
                                                                                        aTCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCggg                                                         >  2:1818899/1‑150 (MQ=255)
                                                                                         tCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCgggg                                                        <  2:1706378/150‑1 (MQ=255)
                                                                                          tttcACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGc                                                       >  1:1095508/4‑150 (MQ=255)
                                                                                           gCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGaatg                                                                                                                                              >  2:3924394/1‑58 (MQ=255)
                                                                                            ccACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTg                                                     >  2:288581/1‑150 (MQ=255)
                                                                                             cACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTgg                                                    >  1:3955875/1‑150 (MQ=255)
                                                                                             cACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTgg                                                    <  2:2573968/150‑1 (MQ=255)
                                                                                              aCCGCCTTCATCGAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTggg                                                   >  1:635962/1‑150 (MQ=255)
                                                                                              aCCGCCTTCATCGAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTggg                                                   >  2:5103727/1‑150 (MQ=255)
                                                                                              aCCGCCTTCATCGAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTggg                                                   >  2:1415383/1‑150 (MQ=255)
                                                                                              aCCGCCTTCATCGAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTggg                                                   >  2:850410/1‑150 (MQ=255)
                                                                                              aCCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTggg                                                   <  1:511121/150‑1 (MQ=255)
                                                                                              aCCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTggg                                                   >  1:4677400/1‑150 (MQ=255)
                                                                                                cGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGa                                                 <  1:4543560/150‑1 (MQ=255)
                                                                                                  ccTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGagt                                               >  1:3790667/1‑148 (MQ=255)
                                                                                                   cTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATtg                                              <  2:2324449/150‑1 (MQ=255)
                                                                                                    ttcatcgaGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATtgt                                             <  1:850410/143‑1 (MQ=255)
                                                                                                    tCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATtgt                                             <  1:5153809/150‑1 (MQ=255)
                                                                                                      gATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTCCGCCATCTGGTCGGGGCTGGGGATTGTGc                                           <  2:3465216/150‑1 (MQ=255)
                                                                                                      gATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGc                                           <  1:3676705/150‑1 (MQ=255)
                                                                                                         gAAGGCGGTGAAAGGCCTGAGCACGCCGTTGACACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGa                                        >  2:2131998/1‑150 (MQ=255)
                                                                                                             gCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAg                                    >  2:1702883/1‑150 (MQ=255)
                                                                                                              cGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAgt                                   <  2:1974739/150‑1 (MQ=255)
                                                                                                                gTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAgtgt                                 <  2:806956/150‑1 (MQ=255)
                                                                                                                   aaaGGCCTGAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTggc                              >  1:1349371/1‑150 (MQ=255)
                                                                                                                           gAGCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTgg                      <  1:4966311/150‑1 (MQ=255)
                                                                                                                             gCACGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTg                    <  2:3982522/150‑1 (MQ=255)
                                                                                                                               aCGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgca                                                                                              >  1:676639/1‑74 (MQ=255)
                                                                                                                               aCGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCgca                                                                                              <  2:676639/74‑1 (MQ=255)
                                                                                                                               aCGCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGAt                  >  1:3732295/1‑150 (MQ=255)
                                                                                                                                 gCCGTTGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCt                >  1:4360897/1‑150 (MQ=255)
                                                                                                                                     ttGCCACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCTATgg            <  1:1313715/150‑1 (MQ=255)
                                                                                                                                        ccACTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCTATGGGCa         >  1:3364278/1‑150 (MQ=255)
                                                                                                                                          aCTGCTGCTGAATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCTATGGGCAGa       >  1:141671/1‑150 (MQ=255)
                                                                                                                                          aCTGCTGCTGAAAGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGACAACTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCTATGGGCAGa       >  2:1922631/1‑150 (MQ=255)
                                                                                                                                            tgctgctgaATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCTATg             <  1:1658832/142‑1 (MQ=255)
                                                                                                                                              ctgctgaATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCTATGGGCAGAAGCt   <  1:2806880/150‑1 (MQ=255)
                                                                                                                                               tgctgaATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCTATGGGCAGAAGCTg  <  2:4828000/150‑1 (MQ=255)
                                                                                                                                               tgctgaATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCTATGGGCAGAAGCTg  >  1:1595699/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGGCAACGTCTAGAATGCCCGCTCAGTAACGGAGGATCCTGCCATGAATGCCTATACCTACCTCGCCATTGCCATCTGCGCCGAAGTCATCGCCACCGCCTCGATGAAGGCGGTGAAAGGCCTGAGCACGCCGTTGCCACTGCTGCTGATGGTGGTTGGCTATGGCATTGCGTTCTGGATGCTGACGCTGGTGGTCCGCAGCATTCCGGTGGGGATTGCCTACGCCATCTGGTCGGGGCTGGGGATTGTGCTGATCAGTGTGGCGGCGCTGGTGATCTATGGGCAGAAGCTG  >  NC_002947/5610594‑5610885

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: